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ARO  PROJECT  No.  DAAG55-98-1-0454 
MICROSTRUCTURE  AND  DYNAMIC  BEHA  VIOR 
CORRELATION  IN  TWO-PHASE  TiB2-Al203  CERAMIC 
(P.I.s  Naresh  Thadhani  and  Min  Zhou,  Georgia  Tech,  Atlanta) 


L  INTRODUCTION 

The  dynamic,  high-strain-rate  deformation  behavior  of  materials  (and  consequently  their 
ballistic  performance)  is  dramatically  influenced  by  their  microstructural  characteristics, 
including  phase  size,  phase  morphology,  composition  and  texture.  Such  a  trend  has  been 
revealed  not  only  for  steels  used  as  heavy  armor,  but  also  for  ceramics  considered  for  lightweight 
armor  applications  [1].  In  addition  to  their  low  density,  ceramics  exhibit  superior  hardness  and 
high  compressive  strength  that  enables  erosion  and  “interface  defeat”  of  projectiles.  This 
mechanism  makes  ceramics  highly  desirable  for  use  as  armor  materials.  Over  the  last  thirty  years 
a  number  of  studies  on  ballistic  performance  and  dynamic  behavior  [1-19]  have  been  performed 
on  various  ceramics  (e.g.,  AIN,  AI2O3,  B4C,  SiC,  TiB2,  WC,  Zr02).  These  studies  have  suggested 
that  the  desirable  characteristics  of  ceramics  that  are  beneficial  for  the  defeat  of  projectiles 
include  the  combination  of  high  yield  strength  or  hardness,  high  tensile  spall  strength,  high 
fracture  toughness,  high  Poisson’s  ratio,  and  high  coefficient  of  friction  [15].  The  properties 
typically  measured  to  characterize  the  dynamic  behavior  of  ceramics  are  the  Hugoniot  Elastic 
Limit  (cthel)  and  spall  strength.  In  general,  the  spall  strength  is  lower  than  the  elastic  limit. 
Materials  such  as  TiB2  have  a  spall  strength  that  decreases  with  increasing  impact  stress  and 
becomes  negligible  at  the  ghel-  In  contrast,  AI2O3  maintains  its  original  spall  strength  even  at 
impact  stress  levels  that  exceed  the  ghel  [3]. 

It  has  been  demonstrated  that  SiC  and  TiB2  exhibit  the  most  desirable  ballistic  properties 
based  on  results  of  experiments  performed  to  determine  the  transition  between  interface  defeat 
and  penetration  involving  W-projectiles  [14-16].  The  impact  velocity  for  dwell/penetration 
transition  for  SiC  and  TiB2  is  higher  than  that  for  B4C,  in  spite  of  the  higher  yield  strength  of  the 
latter.  This  trend  is  consistent  with  the  predictions  from  the  dwell/penetration  transition  model 
[17]  which  accounts  for  damage  mechanism  based  on  the  extension  of  mode-I  wing  cracks.  The 
predictions  illustrate  that  SiC  and  TiB2  ceramics  have  the  ability  to  suppress  the  formation  of 
wing  cracks  due  to  plastic  relaxation  of  pre-existing  flaws,  while  B4C  shows  a  brittle  behavior 
dominated  by  growth  of  wing  cracks.  Brittle  fracture  is  also  expected  to  dominate  the  behavior  of 
AI2O3.  Measurements  of  shock  wave  profiles  have  also  illustrated  the  ability  of  SiC  and  TiB2 
ceramics  to  undergo  deformation-induced  hardening  [18],  which  may  give  rise  to  improved 
ballistic  properties. 

Most  dynamic  behavior  studies  performed  to-date,  have  focused  on  single-phase 
monolithic  ceramics,  although  ceramics  with  glassy  (impurity)  phases,  e.g.,  AD85  A1203,  have 
also  been  investigated.  In  such  ceramics,  the  presence  of  intergranular  oxide  glass  has  been 
shown  to  significantly  lower  the  tensile  (spall)  strength,  and  therefore,  the  fracture  resistance 
[19].  Dynamic  behavior  of  ductile  metal-matrix  composites,  e.g.,  Al-alloys  consisting  of 
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embedded  ceramic  particles,  has  been  studied  to  a  limited  extent  [20].  It  has  been  shown  that 
both,  the  dynamic  yield  strength  and  the  spall  strength  are  reduced,  in  comparison  to  strength 
increases  observed  in  particle-reinforced  composite  materials  under  quasi-static  loading 
conditions.  The  composite  structure  acts  as  a  mechanical  energy  trap  due  to  scattering  of  waves 
from  incoherent  boundaries  and  interfaces  between  the  matrix  and  reinforcement  phases  having 
dissimilar  shock  impedance.  The  high-strain-rate  mechanical  behavior  of  ceramic-ceramic 
composites  shows  a  strain  rate  dependence  significantly  different  from  that  observed  in  metal- 
matrix  composites. 

Dynamic  behavior  of  ceramic-composites  has  been  limited  only  to  the  study  of  two-phase 
Al2C>3+TiB2  ceramics  (70/30  mass  ratio).  Past  work  on  these  two-phase  ceramics  has  revealed  an 
80%  increase  in  compressive  strength  with  increasing  strain  rate  (3.5  GPa  at  10'4  s'1  to  5.8  GPa  at 
103  s'1)  as  shown  in  Fig  1(a)  [21],  In  addition,  as  illustrated  in  Figure  1  (b),  the  two-phase 
AI2O3+T1B2  ceramics  have  superior  static  and  dynamic  mechanical  properties  than  their 
monolithic  constituents  [21-23],  The  AbCh+TiE^  ceramics  have  also  shown  better  penetration 
resistance  than  monolithic  AI2O3,  and  the  system  in  which  Til?2  is  an  interconnected  phase 
surrounding  AI2O3  has  been  shown  to  exhibit  a  superior  ballistic  performance  compared  with  the 
system  in  which  the  two  phases  are  simply  uniformly  interdispersed  [21,24],  Micromechanical 
simulations  have  also  demonstrated  the  effect  of  microstructural  bias  on  failure  resistance 
[25,26].  However,  the  influence  of  micro-structural  bias  on  the  fundamental  dynamic  properties 
of  these  ceramics  has  not  been  fully  established. 

The  objective  of  the  present  work  is  to  quantitatively  characterize  the  microstructure  and 
the  high-strain-rate  deformation  and  damage  response  of  four  types  of  microstructurally-biased, 
two-phase  TiB2+Al203  ceramics.  The  microstructural  bias  of  the  four  ceramics  falls  into  two 
morphological  categories.  The  first  category  involves  a  continuous  (interconnected)  TiB2 
network  that  surrounds  AI2O3  (qualitatively  termed  ‘T@A’)  and  the  second  category  involves 
TiB2  and  AI2O3  phases  that  are  interdispersed  and  uniformly  inter-twined  with  each  other 
(qualitatively  termed  ‘TinA’)  [21,27].  Quantitative  stereology  methods  [28]  were  used  to 
determine  the  two-phase  biased  microstructure.  Normal  plate-impact  experiments,  as  well  as 
Hopkinson  bar  experiments  were  performed  to  measure  the  compressive  yield  strength,  the 
Hugoniot  Elastic  Limit  (chel)  and  the  tensile  spall  strengths  of  the  materials.  The  measured 
responses  were  then  correlated  with  the  microstructure  morphological  characteristics  and  the 
phase  size  scales.  Micromechanical  modeling  based  on  a  cohesive  surface  formulation  model 
that  provides  explicit  account  for  arbitrary  microstructures  and  arbitrary  fracture  patterns  was 
also  developed  to  study  the  deformation  and  damage  sub-structure  evolution  during  impact 
loading.  Deterministic  and  stochastic  characterization  of  fracture  resistance  has  also  been 
performed  and  implemented. 

The  main  findings  of  the  work  performed  in  this  investigation,  are  compiled  in  the 
following  sections.  Only  a  brief  description  is  provided  along  with  a  listing  of  thesis  and  papers 
describing  the  work.  Copies  of  the  abstracts  of  thesis  and  reprints  of  some  of  the  key  papers  are 
attached  in  the  appendix. 


2.  MICROSTRUCTURAL  CHARACTERIZATION  AND  QUANTIFICATION 

Four  types  of  microstructurally-biased  TiB2+Al203  ceramics  have  been  produced:  (a) 
SHS  with  TiB2  surrounding  A1203  (T@A),  (b)  SHS  with  intermixed  TiB2  and  A1203  (TinA),  (c) 
MM  with  TiB2  surrounding  A1203  (T@A),  and  (d)  MM  with  intermixed  TiB2  and  A1203  (TinA). 
An  automated  quantitative  metallography  technique  was  developed  to  determine  the 
connectivity  of  the  TiB2  phase,  based  on  measurement  of  the  average  integral  curvature  defined 
by:  kave  =  2n  N/e!  /Q  =  2x  NAne‘  /(. P ),  +  Pp),  where,  N/et  is  the  net  numbers  of  loops  of  phase  A 

per  unit  area,  H  is  the  mean  linear  intercept,  and  P,,  and  Pp  are  numbers  of  TiB2  particles  and 
A1203  holes,  respectively.  The  number  counts  and  perimeters  for  size  measurement  are 
performed  on  several  montages  of  micrographs  obtained  from  both  planar  and  cross-sectional 
surfaces  using  the  SCOPE-PRO  software  and  MATERIALS-PRO  software,  which  allow  reliable 
and  repeatable  measurements.  Quantitative  metallography  measurements  performed  on  the  four 
sets  of  samples  with  the  different  microstructure  illustrate  that  the  samples  do  not  reveal  a  100% 
microstructural  bias.  The  general  trend  reveals  Sample  ‘C’  having  the  most  interconnected  TiB2 
phase  surrounding  the  A1203  phase,  and  Sample  ‘D’  having  the  least  interconnected  TiB2. 

Published  Work: 

1.  Louis  Ferranti,  Processing  and  Quantitative  Microscopy  Characterization  of  Microstructurally-biased 
TiB2+Al20j  two-phase  ceramic,  M.S.  thesis,  Georgia  Institute  of  Technology,  October  2001. 

2.  G.  Kennedy,  L.  Ferranti,  R.  Russell,  M.  Zhou,  and  N.N.  Thadhani,  Dynamic  High-Strain-Rate  Behavior  of 

Microstrcuturally-biased  two-phase  TiB2+A1203  Ceramics,  Journal  of  Applied  Physics  vol91  no  4 
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3.  L.  Fen-anti  and  N.N.  Thadhani,  “Quantitative  Characterization  of  Microstructurally  Biased  Two-Phase 
TiB2+Al203  Ceramics  Using  Mean  Integral  Curvature,”  submitted  to  Metallurgical  Transactions  October 
2002. 


3.  CHARACTERIZATION  OF  HIGH-STRAIN-RATE  MECHANICAL  PROPERTIES 

Hugoniot  Elastic  Limit:  Measurements  of  the  Hugoniot  Elastic  Limit  (oHel)  were 
obtained  from  gas-gun  impact  experiments  conducted  at  a  nominal  750  m/s  impact  velocity. 
Input-  and  backer-wave  stress  profiles  were  recorded  using  the  PVDF  stress  gauges,  from  which 
the  Hugoniot  elastic  limit  (an el)  and  wave-speed  were  determined.  The  Hugoniot  Elastic  Limit 
(cthel),  a  measure  of  the  axial  stress  at  which  a  solid  loaded  under  uniaxial  compression  can  no 
longer  support  elastic  deformation  and  begins  to  flow  plastically  or  undergoes  fracture,  is 
observed  to  be  a  strong  function  of  the  microstructure  including  the  phase  size.  The  <jHel  is 
found  to  be  the  lowest  for  Sample-B  (4.4+1. 2  GPa),  due  to  its  ~4%  porosity  and  highest  for 
Sample-D  (8. 5+4. 5  GPa),  although  its  standard  deviation  is  also  large.  The  cthel  values  reported 
previous  by  Grady  for  TiB2  is  9-18  GPa,  and  ~6.7  for  A1203.  While  (ctHEl)  identifies  the  limit  of 
elastic  response  in  a  dynamic  (shock)  uniaxial  strain  loading,  it  is  also  common  to  assume  a  von 
Misses  condition  of  yield  which  asserts  that  yield  initiates  when  the  second  deviatoric  stress 
invariant  attains  a  critical  value.  Through  this  formalism,  the  yield  stress  in  simple  tension  can  be 
easily  related  to  (cthel)  through  ays  =  2  (Cs2  /  Cl2)  x  (cthel)-  The  yield  stress  in  simple  tension 
(aYs)  for  the  various  samples  follows  the  same  trend  as  the  Hugoniot  Elastic  Limit. 


Tensile  (spall)  Strength  Properties:  Tensile  (spall)  strength  experiments  were 
performed  (using  VISAR  interferometry)  on  samples  of  microstructure  ‘A’  at  impact  velocities 
corresponding  to  input  stresses  above  and  below  the  HEL,  and  on  samples  of  microstructures  CC5 
and  CD’  under  elastic  loading  conditions  to  ensure  that  the  compression-induced  damage  does  not 
influence  the  tensile  response.  The  spall  strength  for  sample-^’  is  0.320  MPa  at  input  stresses  of 
3.7  and  7.9  GPa.  However,  with  input  stress  increasing  to  11.8  GPa,  the  tensile  spall  strength 
decreases  to  0.160  GPa.  Tensile  spall  experiments  were  also  performed  on  samples  of 
microstructures  CC’  and  CD\  at  an  input  stress  of  -3.8  GPa.  A  spall  strength  of  0.311  GPa  for 
Sample-C  and  0.222  GPa  for  Sample-  D,  was  observed.  The  comparison  of  spall  test  results 
between  Sample-C  and  Sample-D  illustrate  that  the  latter  sample  with  dispersed  microstructure, 
has  a  lower  spall  strength,  while  ceramic  with  TiB2  interconnected  microstructure  has  a  higher 
spall  strength.  The  measured  high  value  of  the  spall  strength  of  the  two-phase  ceramic  is  similar 
to  the  published  spall  strength  of  TiB2  (0.33  GPa)  but  lower  than  that  of  A12C>3  (0.45  GPa). 

Dynamic  Compressive  Strength:  Experimental  measurement  of  dynamic  compressive 
strength  has  been  carried  out  at  strain  rates  of  103  s’1,  using  the  split  Hopkinson  compression  bar 
configuration.  Experiments  conducted  so  far  show  that  the  compressive  strength  of  the  four 
materials  is  between  4-5.2  GPa.  There  is  clearly  a  dependence  of  strength  on  microstructure.  It  is 
found  that  there  is  up  to  15%  difference  in  strength  between  the  four  microstructures,  and  order 
of  decreasing  strength  is  similar  to  that  predicted  by  numerical  simulations. 
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4.  MICROMECHANICAL  MODELING  OF  DYNAMIC  BEHAVIOR 

A  micromechanical  model  that  explicitly  describes  arbitrary  microstructures  and  arbitrary 
fracture  patterns  was  formulated  using  the  cohesive  finite  element  (CFEM)  method.  The 
approach  uses  both  a  constitutive  law  for  the  bulk  solid  constituents  and  a  constitutive  law  for 
fracture  surfaces.  The  model  based  on  a  phenomenological  characterization  for  atomic  forces  on 
potential  crack/microcrack  surfaces  was  used  to  analyze  dynamic  crack  propagation  in  the  four 
types  of  TiB2+Al203  samples.  The  analysis  considered  arbitrary  phase  distributions  in  the  actual 
microstructures  of  the  two-phase  ceramics.  Numerical  results  reveal  that  interfacial  bonding 
strength  significantly  influences  failure.  At  the  same  load  intensity,  weak  interfacial  bonding 
strengths  give  rise  to  microcrack  nucleation  and  coalescence.  In  the  presence  of  strong  interfacial 
bonding  strength,  the  extension  and  branching  of  a  main  crack  is  found  to  be  the  primary  failure 
mode.  Crack  speeds  are  found  to  be  <60%  of  the  Rayleigh  wave  speed  of  the  more  compliant 
constituent  in  the  microstructure.  Analyses  of  the  effects  of  microstructure  on  dynamic  fracture 
resistance  showed  that  (a)  fine  phase  size  is  beneficial  to  enhancement  of  fracture  resistance  and 
significant  gains  (up  to  60%)  in  fracture  resistance  can  be  realized  though  manipulation  of  phase 
morphology  and  phase  distribution. 
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6.  SUMMARY  OF  RESULTS 
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evolution  during  impact  loading.  The  Ph.D.  and  M.S.  thesis,  as  well  as  publications  and 
presentations  resulting  from  this  work,  are  listed  below. 

7.  GRADUATE  THESIS  (COMPLETED  and  IN  PROGRESS) 

1.  Jun  Zhai,  “Micromechanical  Modeling  of  Dynamic  Fracture  in  Heterogeneous  Materials,” 
Ph.D.  dissertation,  Georgia  Institute  of  Technology,  January  2000.  {currently  Research 
Scientist  at  Motorola) 

2.  A.  Keller,  “An  Experimental  Analysis  of  the  Dynamic  Failure  Resistance  of  TiB2/Al203 
Composites,”  Georgia  Tech  M.S.  Thesis,  October  2000  (currently  at  Ford  Motor  Company) 

3.  Louis  Ferranti,  Processing  and  Quantitative  Microscopy  Characterization  of 
Microstructurally-biased  TiB2+Al203  two-phase  ceramic,  M.S.  thesis,  October  2001 
(currently  Ph.D.  student  at  Georgia  Tech  working  on  graduate  research  internship  through 
EGLIN  AFB,  Florida). 

4.  Greg  Kennedy,  “Effect  of  Microstructural  Bias  on  the  High  Strain  Rate  Properties  of 
TiB2+Al203  two-phase  ceramic,”  M.S.  thesis  (in  final  stage  of  thesis  writing;  currently  on  a 
3-month  international  research  experience  program  funded  through  an  award  from 
Kumamoto  University  in  Japan). 

5.  Vikas  Tomar,  “Stochastic  analysis  of  failure  resistance  of  two-phase  ceramics”,  PhD  thesis  in 
progress. 

8.  LIST  OF  PUBLICATIONS 

1.  J.  Zhai  and  M.  Zhou,  “Finite  Element  Analysis  of  Micromechanical  Failure  Modes  in 
Heterogeneous  Brittle  Solids”,  International  Journal  of_ Fracture,  special  issue  on  Failure 
Mode  Transition  in  Solids,  R.  C.  Batra,  Y.  D.  S.  Rajapakse,  and  A.  J.  Rosakis,  eds.,  101,  pp. 
161-180,2000; 

2.  M.  Zhou  and  J.  Zhai,  Micromechanical  Characterization  of  the  Fracture  Resistance  of  A 
TiB2/Al203  Ceramic  Composite  System,  American  Ceramic  Society  102nd  Annual  Meeting 
&  Exposition,  St.  Louis,  Missouri,  April  30  -  May  3,  2000; 

3.  G.  Kennedy,  J.  Zhai,  R.  Russell,  M.  Zhou,  and  N.  Thadhani,  High-Strain-Rate  Behavior  of 
Microstructurally-biased  Two-Phase  TiB2+Al203  Ceramics,  Proceedings^  of  the  24th 
Annual  Cocoa  Beach  Conference  of  The  American  Ceramic  Society,  January  23-28,  2000; 

4.  Greg  Kennedy,  Rod  Russell,  Louis  Ferranti,  and  Naresh  Thadhani,  “Fabrication  of 

Microstructurally-biased  Two-phase  TiB2+Al203  Ceramics  and  their  Dynamic  Properties,”  in 


Processing  and  Fabrication  of  Advanced  Materials,  ASM  Conf.  Proc.,  eds.  T.S.  Srivatsan  and 
R.A.  Varin. 

5.  Greg  Kennedy,  Andrew  Keller,  Rod  Russell,  Louis  Ferranti,  Jun  Zhai,  Min  Zhou,  and  Naresh 

Thadhani;  “Dynamic  Mechanical  Properties  of  Microstructurally-biased  Two-phase 
TiB2+Al203  Ceramics,”  in  Proc.  of  EXPLOMET  2000  Conference,  eds.  K.P.  Staudhammer, 
L.E.  Murr,  and  M.A.  Meyers. 

6.  G.  Kennedy,  L.  Ferranti,  R.  Russell,  M.  Zhou,  and  N.  Thadhani,  Dynamic  High-Strain-Rate 

Mechanical  Behavior  of  Microstructurally-biased  Two-Phase  TiB2+Al203  Ceramics,  J. 
Applied  Physics,  91(4),  1921-1927,  2002;  ’  ~ 

7.  A.  R.  Keller  and  M.  Zhou,  Effect  of  Microstructure  on  Dynamic  Failure  Resistance  of 
TiB2/Al203  Ceramics,  submitted  to  Journal  of  theAmerican  Ceramic  Society,  2001; 

8.  J.  Zhai,  V.  Tomar  and  M.  Zhou,  (2002)  Explicit  Analysis  of  Microscopic  Fracture  Processes 
Using  a  Cohesive  Finite  Element  Method,  Part  I-Mesh  Dependence  and  Convergence, 
manuscript  in  preparation. 

9.  J.  Zhai,  V.  Tomar  and  M.  Zhou,  (2002)  Explicit  Analysis  of  Microscopic  Fracture  Processes 
Using  a  Cohesive  Finite  Element  Method,  Part  II-Numerical  Results,  manuscript  in 
preparation. 

10.  A.  Keller  and  M.  Zhou,  Microscopic  and  Nanoscopic  Fracture  in  Heterogeneous  in  a 
Heterogeneous  Material  System,  Proceedings  of  2002  SEM  Annual  Conf.  and  Expo,  on  Exp. 
&  Appl.  Mechanics,  June  10-12, 2002; 

11.  A.  Keller,  G.  Kennedy,  L.  Ferranti,  M.  Zhou,  and  N.  Thadhani,  Correlation  of  Dynamic 
Behavior  with  Microstructural-bias  in  Two-phase  TiB2+A1203  Ceramics,  Impact  Engineering 
and  Application,  Proceedings  of  the  Fourth  International  Symposium  on  Impact  Engineering 
(ISIE/4),  pp.  321-326,  July  16-18,  Kumamoto,  Japan,  Elsevier,  2001; 

12.  M.  Zhou,  Numerical  Simulation  of  Micro  and  Nano  Scale  Fracture,  Proceedings  of  The  10th 
International  Congress  of  Fracture,  December  2-7, 2001,  Honolulu,  Hawaii; 

13.  A.  R.  Keller  and  M.  Zhou,  Experimental  characterization  of  the  dynamic  failure  of 
TiB2/A1203  composites,  Shock  compression  of  condensed  matter-2001,  ed.  M.  D.  Furnish, 
N.  N.  Thadhani,  and  Y.  Horie,  American  Institute  of  Physics,  795-798, 2001 ; 

14.  G.  Kennedy,  L.  Ferranti,  R.  Russell,  M.  Zhou,  and  N.  Thadhani,  Influence  of  microstructural 
bias  on  the  hugoniot  elastic  limit  and  spall  strength  of  two-phase  TiB2+A1203  ceramics, 
Shock  compression  of  condensed  matter-2001,  ed.  M.  D.  Furnish,  N.  N.  Thadhani,  and  Y. 
Horie,  American  Institute  of  Physics,  755-758, 2001; 

15.  G.  Kennedy  and  N.N.  Thadhani,  Investigation  of  the  spall  strength  of  microstructurally- 
biased  two-phase  TiB2+A1203  ceramics,  2001  (in  preparation); 

9.  PRESENTATIONS:  (♦invited  talk! 

1.  Greg  Kennedy,  Micromechanical  Characterization  of  the  Fracture  Resistance  of  A 
TiB2/Al203  Ceramic  Composite  System,  102nd  Annual  Meeting  &  Exposition  of  the 
American  Ceramic  Society,  St.  Louis,  MO,  April  30  -  May  3, 2000; 

2.  Min  Zhou,  Effect  of  Microstructure  on  Dynamic  Fracture  Resistance  of  Ceramic  Composites, 
The  20' ’  Southeastern  Conference  on  Theoretical  and  Applied  Mechanics,  April  16-18, 2000, 
Callaway  Gardens,  GA; 

3.  Computational  Characterization  of  the  effect  of  Phase  Morphology  and  Interphase  bonding 
on  Fracture  Toughness  of  Ceramic  Composites,  ASME IMECE,  Orlando,  FL,  Nov.  2000; 


4.  Naresh  Thadhani  (♦),  “Fabrication  of  Microstructurally-biased  Two-phase  TiB2+Al203 
Ceramics  and  their  Dynamic  Properties,”  in  Ninth  Inti.  Conf.  on  Processing  and  Fabrication 
of  Advanced  Materials,  ASM  Meeting,  St.  Louis,  October  2000. 

5.  Naresh  Thadhani;  “Dynamic  Mechanical  Properties  of  Microstructurally-biased  Two-phase 
TiB2+Al203  Ceramics,”  EXPLOMET  2000  Conference,  Albuquerque,  June  2000. 

6.  Naresh  Thadhani  “Influence  Of  Microstructural  Bias  On  The  Hugoniot  Elastic  Limit  And 
Spall  Strength  Of  Two-Phase  TiB2+Al203  Ceramics,”  at  the  APS  Topical  Group  Meeting  in 
Shock-compression  of  Condensed  Matter  -  2001,  June  24-29,  2001. 

7.  M.  Zhou  (♦),  Microscopic  and  Nanoscopic  Fracture  in  Heterogeneous  in  a  Heterogeneous 
Material  System,  2002  SEM  Annual  Conf.  and  Expo,  on  Exp.  &  Appl.  Mechanics,  June  10- 
12,  2002; 

8.  M.  Zhou  (♦),  Computational  and  Experimental  Characterization  of  Micro  and  Nano  Scale 
Deformation  and  Fracture,  NASA  Langley  Research  Center,  Hampton,  August  2001 

9.  M.  Zhou,  Correlation  of  Dynamic  Behavior  with  Microstructural-bias  in  Two-phase 
TiB2+Al203  Ceramics,  The  Fourth  International  Symposium  on  Impact  Engineering  (ISIE/4), 
July  16-18,  2001,  Kumamoto,  Japan; 

10.  Numerical  Simulation  of  Micro  and  Nano  Scale  Fracture,  The  10th  International  Congress  of 
Fracture,  December  2-7,  2001,  Honolulu,  Hawaii; 

11.  Experimental  Characterization  of  the  Dynamic  Failure  Resistance  of  Effect  TiB2/Al203 
Composites,  APS  SCCM  meeting,  June  27-29,  2001,  Atlanta,  GA; 

12.  A.  Keller  and  M.  Zhou,  Microscopic  and  Nanoscopic  Fracture  in  Heterogeneous  in  a 
Heterogeneous  Material  System,  2002  SEM  Annual  Conf.  and  Expo,  on  Exp.  &  Appl. 
Mechanics,  June  10-12,  2002; 

10.  LIST  OF  PERSONNEL  INVOLVED 

Faculty:  Prof.  Naresh  N.  Thadhani  and  Prof.  Min  Zhou 

Grad  Students:  Jun  (Charlie)  Zhai  (Ph.D./ME),  Andrew  Keller  (M.S./ME), 

Greg  Kennedy  (M.S./MSE)  and  Louis  Ferranti  (M.S./MSE),  Vikas  Tomar  (ME) 
UG  Students:  Andy  Gapin  (UG/MSE),  Thomas  Frosell  (ME),  and  Nicholas  Liggett  (ME/MSE) 
Summer  Undergraduate  Research  Assistants:  Reinhard  Powell  (Georgia  Tech,  1999),  Joshua 
Clarke  (Georgia  Tech,  2000),  Celeste  Davis  (California  State  University,  Pamona,  2001),  Adam 
Rowe  (Colorado  School  of  Mines,  2002) 

High  School  Teacher  Participants:  Marie  Johnson  (Chamblee  HS,  Atlanta,  2000  and  2001) 

11.  LIST  OF  HONORS  AND  AWARDS 

•  Prof.  Min  Zhou  was  awarded  the  NSF  Career  Award 

•  Prof.  Naresh  Thadhani  was  elected  FELLOW  of  ASM  International 

•  Greg  Kennedy,  awarded  ASM  Atlanta  Chapter  outstanding  graduate  student  award 

•  Greg  Kennedy,  International  Research  Fellowship  Award  for  research  at  Kumamoto 
University  in  Japan. 

12.  INTERACTION  WITH  ARL  SCIENTISTS 

Interactions  have  been  maintained  with  Ernie  Chin,  Jim  McCauley,  Datta  Dandekar,  Jane 
Adams,  and  Gary  Guildy  at  ARL.  Through  the  financial  support  provided  by  ARL,  Greg 
Kennedy  was  able  to  spend  two  2-month  internships  at  ARL.  Greg  worked  on  microstructural 
characterization  of  samples  hot-pressed  at  ARL. 


REFERENCES 

1.  D.J.  Viechnicki,  M.L.  Slavin,  and  M.I.  Kliman,  Ceram.  Bull.,  70  (1991),  p.  1035. 

2.  D.A.  Shockey,  A.H.  Marchand,  S.R.  Skaggs,  G.E.  Cort,  M.W.  Burkett,  and  R.  Parker,  Int.  J. 
of  Impact  Eng.,  Vol.  9(3),  (1990)  p.  263. 

3.  D.E.  Grady,  “Dynamic  Material  Properties  of  Armor  Ceramics,”  Sandia  National 
Laboratories  Reports,  SAND91-0147  (March  1991). 

4.  D.E.  Grady  and  J.L.  Wise,  Sandia  National  Laboratories  Reports,  SAND93-0610  (Sept  ‘93). 

5.  D.P.  Dandekar  and  D.C.  Benfanti,  J.  Appl.  Phys.,  73,  (1993)  pp.  673-679. 

6.  D.E.  Grady,  Sandia  National  Laboratories  Reports,  SAND94-3266  (February  1995). 

7.  D.L.  Orphal,  R.R.  Franzen,  A.C.  Charters,  T.L.  Menna,  and  A.J.  Piekutowski,  Int.  J.  Impavct 
Eng.,  Vol.  19  (1997)  p.  15. 

8.  R.  Feng,  Y.M.  Gupta,  and  M.K.  Wong,  J.  Appl.  Phys.,  82  (1997)  pp.  2845-2854. 

9.  J.  Lankford,  W.W.  Predebon,  J.M.  Staehler,  G.  Subhash,  B.J.  Plekta,  and  C.E.  Anderson, 
Mech.  Matls.,  29  (1998)  p.  205. 

10.  G.R.  Johnson  and  T.J.  Holmquist,  J.  Appl.  Phys.,  85(12)  (1999)  p.  8060. 

11.  D.P.  Dandekar  and  P.I.  Bartkowski,  “Spall  Strength  of  Silicon  Carbides  Under  Shock 
Loading,”  in  Fundamental  Issues  and  Applications  of  Shock  Wave  and  High  Strain  rate 
Phenomena,  eds.  K.P.  Staudhammer,  L.E.  Murr,  and  M.A.  Meyers,  Elsevier,  (2001),  71-77. 

12.  G.E.  Hauver,  P.H.  Netherwood,  R.F.  Blenck,  and  L.K.  Kecskes,  19th  Army  Science 
Conference,  Florida,  June  20-24, 1994. 

13.  H.D.  Espinosa,  N.S.  Brar,  G.  Yuan,  Y.  Xu,  and  V.  Arrietta,  Int.  J.  Solids  Structures,  37  [36], 
(2000),  p.  4893. 

14.  F.I.  Grace,  Interface  Defeat  of  impacting  rods  against  ceramic  targets,”  in  Fundamental 
Issues  and  Applications  of  Shock  Wave  and  High  Strain  rate  Phenomena,  eds.  K.P. 
Staudhammer,  L.E.  Murr,  and  M.A.  Meyers,  Elsevier  Science,  (2001),  pp.  421-428. 

15.  J.C.  LaSalvia,  E.J.  Horwarth,  E.  Rapacki,  C.  James  Shih,  and  M.A.  Meyers,  “Microstructural 
and  micromechanical  aspects  of  ceramic/long-rod  projectile  interactions:  dwell/penetration 
transitions,”  in  Fundamental  Issues  and  Applications  of  Shock  Wave  and  High  Strain  rate 
Phenomena,  eds.  K.P.  Staudhammer,  et  al.,  Elsevier,  (2001),  437-446. 

16.  P.  Lundberg,  R.  Renstrom,  and  B.  Lundberg,  Int.  J.  Impact  Eng.,  24  (2000)  259. 

17.  H.  Horii  and  S.  Nemat-Nasser,  Phil.  Trans.  R.  Soc.  London,  A,  319  (1986),  337. 

18.  D.E.  Grady,  Mechanics  of  Materials,  Vol.  29,  (1998)  pp.  181-203. 

19.  S.  Sundaram  and  R.J.  Clifton,  Mechanics  of  Materials,  Vol.  29,  (1997),  pp.  233-251. 


20.  J.N.  Johnson,  R.S.  Hixon,  and  G.T.  Gray,  J.  Appl.  Phys.,  76  (10)  (1994)  5706-5718. 

21.  K.V.  Logan,  Ph.D.  dissertation,  Georgia  Institute  of  Technology,  1993. 

22.  A.  Keller,  “An  Experimental  Analysis  of  the  Dynamic  Failure  Resistance  of  T^/ALCE 
Composites,”  Georgia  Tech  M.S.  Thesis,  October  2000. 

23.  Andrew  Keller,  Greg  Kennedy,  Louis  Ferranti,  Min  Zhou,  and  Naresh  Thadhani, 
“Correlation  of  Dynamic  Behavior  With  Microstructural-Bias  In  Two-Phase  TiB2+AL03 
Ceramics,”  to  be  published  in  Proc.  of  Fourth  International  Symposium  on  Impact 
Engineering,  Kumamoto,  Japan,  July  2001 . 

24.  G.  Gilde,  J.W.  Adams,  M.  Burkins,  M.  Motyka,  P.J.  Patel,  E.  Chin  M.  Sutaria,  M.  Rigali,  and 
L.  Prokurat  Franks,  “Processing  of  Aluminum  oxide  and  Titanium  Diboride  Composites  for 
Penetration  Resistance,”  in  Proceedings  of  the  25th  Annual  Cocoa  Beach  Conference  of  The 
American  Ceramic  Society,  January  23-28,  2001. 

25.  J.  Zhai  and  M.  Zhou,  “Finite  Element  Analysis  of  Micromechanical  Failure  Modes  in 
Heterogeneous  Brittle  Solids”,  International  Journal  of  Fracture,  special  issue  on  Failure 
Mode  Transition  in  Solids,  R.  C.  Batra,  Y.  D.  S.  Rajapakse,  and  A.  J.  Rosakis,  eds.,  101,  pp. 
161-180,2000. 

26.  M.  Zhou  and  J.  Zhai,  Micromechanical  Characterization  of  the  Fracture  Resistance  of  A 
TiB2/Al203  Ceramic  Composite  System,  American  Ceramic  Society  102nd  Annual  Meeting 
&  Exposition,  St.  Louis,  Missouri,  April  30  -  May  3,  2000. 

27.  K.V.  Logan,  “Shaped  Refractory  Products  and  the  Method  of  Making  the  Same,”  U.S. 
Patent,  #  5,141,900,  August  25,  1992 

28.  E.E.  Underwood,  Quantitative  Stereology,  Addison  Wesley  Publishing  Company,  1970. 

29.  G.  Kennedy,  J.  Zhai,  R.  Russel,  M.  Zhou,  and  N.N.  Thadhani,  “Dynamic  High  Strain  Rate 
Mechanical  Behavior  of  Microstructurally-Biased  Two-phase  T^+ALCE  Ceramic,”  24th 
Annual  Int.  Conf.  On  Advanced  Ceramics  and  Composites,  Cocoa  Beach,  Jan  23-28,  2001. 


APPENDIX 

RE-  and  PRE-PRINTS 
of 

PAPERS  PUBLISHED 
Based  on 

ARO-SUPPORTED  RESEARCH 


Processing  and  Characterization  of  Microstructually  Biased  Two-Phase 
Titanium  Diboride/Alumina  Ceramic  (TiB2+Al203) 


A  Thesis 
Presented  to 
the  Academic  Faculty 

by 

Louis  Ferranti  Jr. 


In  Partial  Fulfillment 
of  the  Requirements  for  the  Degree 
Master  of  Science  in  Materials  Science  and  Engineering 


Georgia  Institute  of  Technology 
October  2001 


SUMMARY 


The  development  of  microstructurally  biased  two-phase  titanium  diboride/alumina 
(TiB2+Al203)  ceramic  is  influenced  by  several  factors  during  processing.  A  standard  method 
for  quantitatively  characterizing  the  different  types  of  microstructures  produced  is  necessary  to 
comprehend  the  mechanism(s)  that  control  the  microstructural  bias  (e.g.  phase  distribution 
and  morphology)  generated  during  processing,  and  their  subsequent  effect  on  mechanical 
properties.  In  this  study,  two-phased  TiB2+Al203  ceramic  powders  were  fabricated  using  self- 
propagating  high  temperature  synthesis  (SHS)  reactions,  and  by  manually  mixing  (MM) 
commercially  purchased  pure  TiB2  and  A1203  powders.  Dense  ceramic  disks  were  formed  by 
hot-pressing  MM  and  SHS  formed  powders  after  various  ball-milling  times.  The 
microstructures  developed  as  a  function  of  powder  processing  technique  (SHS  reacdon  or 
MM),  ball-milling  time,  and  varying  hot-pressing  parameters  including  varying  degrees  of 
connectivity  of  the  TiB2  minor  phase  surrounding  A1203  and  different  phase  sizes  and  size 
distributions.  These  features  were  characterized  by  stereological  measurements  of  two- 
dimensional  cross-sections  of  the  samples,  based  on  phase  size  and  morphology  in  terms  of 
phase  connectivity  and  dispersion  of  the  minor  phase  (TiR).  Mean  integral  curvature  was 
used  to  quantitatively  describe  the  connectivity  of  the  TiB2  phase  around  an  A1203  matrix.  It 
was  observed  that  the  TiB2  phase  connectivity  is  established  direedy  during  the  SHS  synthesis 
process,  and  extended  periods  of  ball-milling  only  reduce  the  inherent  connectivity  while 
decreasing  the  phase  size  and  altering  the  phase  distribution.  Thus,  an  optimal  hall-milling 
time  of  about  4  hours  for  the  SHS  synthesized  powders  and  lesser  for  the  manually  mixed 


powders  was  needed  to  achieve  the  maximum  degree  of  TiB2  phase  connectivity  and 
phase/particle  size  distribution  for  maximum  final  densification  during  hot-pressing. 
Microhardness  and  fracture  toughness  measurements  performed  on  the  dense  compacts 
showed  that  the  presence  of  TiB2  dramatically  improves  the  fracture  toughness  of  the  two- 
phase  ceramic  in  contrast  to  monolithic  alumina. 
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Summary 


Recently,  alumina/titanium  diboride  (A^O/TiBz)  two  phase  ceramic  composites  having 
a  wide  range  of  micro  and  nano  phase  sizes  and  phase  morphologies  have  been  developed, 
variations  in  processing  conditions  hav4e  allowed  formation  of  these  ceramics  with  biased 
microstructures  which  are  expected  to  influence  the  dynamic  failure  behavior.  In  order  to 
understand  failure  mechanisms  in  this  material  system  and  the  influence  of  phase  morphologies 
and  phase  size  on  fracture  resistance,  a  numerical  framework  is  needed  to  explicitly  account  for 
arbitrary  microstructures  and  fracture  patterns. 

Micromechanical  modeling  and  simulation  provide  an  important  approach  for  analyzing 
the  effects  of  material  inhomogeneity  and  anisotropy  over  a  range  of  microscopic  length  scales.  A 
framework  is  proposed  in  this  research  for  explicit  modeling  and  simulation  of  microscopic 
damage/fracture/failure  processes.  The  model  and  approach  account  for  the  real  arbitrary 
microstructural  morphologies.  A  cohesive  finite  element  method  (CFEM)  based  on  cohesive 
surface  theory  is  used.  A  fully  dynamic  kinetic  framework  and  finite  deformation  kinematic 
formulation  are  used.  Mesh  independence  of  solution  is  studied  and  verified. 

Idealized  microstructures  containing  circular  and  elliptical  particles  and  real 
microstructures  with  arbitrary  morphologies  are  used  to  investigate  the  effects  of  phase 
morphologies,  phase  size  and  phase  anisotropy  on  fracture  of  this  ceramic  composite  system. 
Numerical  results  show  that  microstructural  variations  give  rise  to  a  range  of  fracture  resistance. 
It  is  found  that  small  inter-phase  spacing,  small  phase  size,  uniform  distribution  of  reinforcements 
and  microstructures  with  enhanced  bridging  effects  will  contribute  to  increasing  fracture 
resistance  of  materials. 

The  failure  mode  is  found  to  be  significantly  influenced  by  the  interfacial  bonding 
strength  between  the  phases.  Two  distinct  failure  modes  are  observed  for  strong  and  weak 
bonding.  While  the  primary  failure  mode  is  the  main  crack  propagation  and  branching  in 
microstructures  with  strong  interfaces,  the  primary  failure  mode  is  microcrack  nucleation  and 
coalescence  in  microstructures  with  weak  interfaces. 

The  loading  rate  effects  are  also  investigated.  It  is  found  that  energy  release  rate  increases 
with  loading  rates,  due  to  the  extensive  damage  development  at  higher  loading  rates. 


This  research  contributes  to  the  elucidation  of  microscopic  deformation  and  failure 
mechanisms  of  microscopically  inhomogeneous  materials  over  a  range  of  size  scales.  The 
research  is  important  for  the  development  of  more  failure  resistant  ceramic  composites  through 
microstructura!  design  and  revision. 
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under  quasi-static  conditions.  Clearly,  there  is  a  strong  dependence  of  failure  stress  on 
the  strain  rate.  The  strength  at  higher  rates  is  approximately  21%  greater  than  that 
observed  at  quasi-static  loading  rates. 

The  use  of  a  soft  recovery  mechanism  in  the  experiment  allowed  the  specimen 
fragments  to  be  recovered  for  postmortem  analysis.  Scanning  electron  microscopy 
(SEM)  and  Energy  Dispersive  Spectrometry  (EDS)  were  used  to  analyze  the  post 
mortem  fracture  surfaces.  Inspection  of  the  fracture  surfaces  indicated  that  failure 
associated  with  the  AI2O3  phase  of  the  composite  is  characterized  by  transgranular 
cleavage  in  all  four  materials.  Depending  upon  the  microstructure,  failure  associated 
with  the  TiB2  phase  was  characterized  by  a  combination  of  transgranular  cleavage  and 
intergranular  pull  out. 

In  addition  to  the  SEM  analysis  of  the  post  mortem  fragments,  a  quantification 
of  the  fragment  size  distribution  was  conducted  for  each  of  the  four  materials. 
Materials  with  larger  grain  sizes  and  higher  strength  tended  to  yield  larger  average 
fragment  sizes.  The  distribution  of  the  fragment  sizes  were  compared  to  theoretical 
models  suggested  by  Lienau  (1936),  Mott  and  Linfoot  (1943),  and  Grady  and  Kipp 
(1985,  1990).  It  was  found  that  the  models  of  Grady  and  Kipp  provided  the  best  fit  to 
the  distribution  of  the  fragments  in  these  experiments. 
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The  dynamic  high-strain-rate  behavior  of  two-phase  TiB2-FAl203  ceramics  with  biased 
microstructures  was  investigated  in  this  study.  The  microstructural  bias  includes  differences  in  phase 
(grain)  size  and  phase  distribution  such  that  in  one  case  a  continuous  (interconnected)  TiB2  network 
surrounds  the  A1203  phase  (qualitatively  termed  “T@A”)  and  in  another  case  the  TiB2  and  A1203 
phases  are  interdispersed  and  uniformly  intertwined  with  each  other  (qualitatively  termed  “TinA”). 

Quantitative  microscopy  was  used  to  characterize  the  phase  size  and  the  integral  curvature  which  is 
taken  as  a  measure  of  TiB2  phase  connectivity  around  A1203 .  Dynamic  compression  and  tension 
(spall)  properties  were  measured  using  plate  impact  experiments.  The  measurements  used 
piezoelectric  polyvinyldine  fluoride  stress  gauges  to  obtain  the  loading  profile  and  to  determine  the 
Hugoniot  elastic  limit.  In  addition,  velocity  interferometry  system  for  any  reflector  interferometry 
was  used  to  obtain  the  spall  signal  and  the  tensile  dynamic  strengths  of  the  materials.  Experimental 
results  reveal  that  while  the  crHEL  and  the  compressive  strengths  of  TiB2  +  Al203  ceramics  are 
dependent  on  the  average  grain  (phase)  size,  the  tensile  (spall)  strength  scales  with  the  TiB2 -phase 
connectivity.  This  result  suggests  that  the  interconnected  TiB2/Al203  microstructural  morphology 
provides  a  stronger  impediment  to  failure  in  tension  compared  with  the  morphology  with  simply 
interdispersed  phases.  ©  2002  American  Institute  of  Physics.  [DOI:  10.1063/1.1429770] 


I.  INTRODUCTION 

The  dynamic,  high-strain-rate  deformation  behavior  of 
materials  (and  consequently  their  ballistic  performance)  is 
dramatically  influenced  by  their  microstructural  characteris¬ 
tics,  including  phase  size,  phase  morphology,  composition 
and  texture.  Such  a  trend  has  been  revealed  not  only  for 
steels  used  as  heavy  armor,  but  also  for  ceramics  considered 
for  lightweight  armor  applications.1  In  addition  to  their  low 
density,  ceramics  exhibit  superior  hardness  and  high  com¬ 
pressive  strength  that  enables  erosion  and  “interface  defeat” 
of  projectiles.  This  mechanism  makes  ceramics  highly  desir¬ 
able  for  use  as  armor  materials.  Over  the  last  30  years  a 
number  of  studies  on  ballistic  performance  and  dynamic 
behavior1-19  have  been  performed  on  various  ceramics  (e.g., 
AIN,  AI203,  B4C,  SiC,  TiB2,  WC,  Zr02).  These  studies 
have  suggested  that  the  desirable  characteristics  of  ceramics 
that  are  beneficial  for  the  defeat  of  projectiles  include  the 
combination  of  high  yield  strength  or  hardness,  high  tensile 
spall  strength,  high  fracture  toughness,  high  Poisson’s  ratio, 
and  high  coefficient  of  friction.15  The  properties  typically 
measured  to  characterize  the  dynamic  behavior  of  ceramics 
are  the  Hugoniot  elastic  limit  (rxHEL)  an^  spall  strength.  In 
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general,  the  spall  strength  is  lower  than  the  elastic  limit.  Ma¬ 
terials  such  as  TiB2  have  a  spall  strength  that  decreases  with 
increasing  impact  stress  and  becomes  negligible  at  the  crHEL . 
In  contrast,  A1203  maintains  its  original  spall  strength  even 
at  impact  stress  levels  that  exceed  the  crHEL.3 

It  has  been  demonstrated  that  SiC  and  TiB2  exhibit  the 
most  desirable  ballistic  properties  based  on  results  of  experi¬ 
ments  performed  to  determine  the  transition  between  inter¬ 
face  defeat  and  penetration  involving  W  projectiles.14-16  The 
impact  velocity  for  dwell/penetration  transition  for  SiC  and 
TiB2  is  higher  than  that  for  B4C,  in  spite  of  the  higher  yield 
strength  of  the  latter.  This  trend  is  consistent  with  the  predic¬ 
tions  from  the  dwell/penetration  transition  model17  which  ac¬ 
counts  for  damage  mechanism  based  on  the  extension  of 
mode-I  wing  cracks.  The  predictions  illustrate  that  SiC  and 
TiB2  ceramics  have  the  ability  to  suppress  the  formation  of 
wing  cracks  due  to  plastic  relaxation  of  preexisting  flaws, 
while  B4C  shows  a  brittle  behavior  dominated  by  growth  of 
wing  cracks.  Brittle  fracture  is  also  expected  to  dominate  the 
behavior  of  A1203 .  Measurements  of  shock  wave  profiles 
have  also  illustrated  the  ability  of  SiC  and  TiB2  ceramics  to 
undergo  deformation-induced  hardening,18  which  may  give 
rise  to  improved  ballistic  properties. 

Most  dynamic  behavior  studies  performed  to  date  have 
focused  on  single-phase  monolithic  ceramics,  although  ce- 
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ramies  with  glassy  (impurity)  phases,  e.g.,  AD85  A1203, 
have  also  been  investigated.  In  such  ceramics,  the  presence 
of  intergranular  oxide  glass  has  been  shown  to  significantly 
lower  the  tensile  (spall)  strength,  and,  therefore,  the  fracture 
resistance.19  Dynamic  behavior  of  ductile  metal-matrix  com¬ 
posites,  e.g.,  Al  alloys  consisting  of  embedded  ceramic  par¬ 
ticles,  has  been  studied  to  a  limited  extent.20  It  has  been 
shown  that  both  the  dynamic  yield  strength  and  the  spall 
strength  are  reduced,  in  comparison  to  strength  increases  ob¬ 
served  in  particle-reinforced  composite  materials  under  qua¬ 
sistatic  loading  conditions.  The  composite  structure  acts  as  a 
mechanical  energy  trap  due  to  scattering  of  waves  from  in¬ 
coherent  boundaries  and  interfaces  between  the  matrix  and 
reinforcement  phases  having  dissimilar  shock  impedance. 
The  high-strain-rate  mechanical  behavior  of  ceramic- 
ceramic  composites  shows  a  strain  rate  dependence  signifi¬ 
cantly  different  from  that  observed  in  metal-matrix  compos¬ 
ites.  Dynamic  behavior  of  ceramic  composites  has  been 
limited  only  to  the  study  of  two-phase  Al203  +  TiB2  ceram¬ 
ics  (70/30  mass  ratio).  Past  work  on  these  two-phase  ceram¬ 
ics  has  revealed  an  80%  increase  in  compressive  strength 
with  increasing  strain  rate  (3.5  GPa  at  10~4  s_l-5.8  GPa  at 
103  s'3)  as  shown  in  Fig.  1(a).21  In  addition,  as  illustrated  in 
Fig*  Kb),  the  two-phase  Al203  +  TiB2  ceramics  have  static 
and  dynamic  mechanical  properties  superior  to  their  mono¬ 
lithic  constituents.21"23  The  Al203  +  TiB2  ceramics  have  also 
shown  better  penetration  resistance  than  monolithic  A1203, 
and  the  system  in  which  TiB2  is  an  interconnected  phase 
surrounding  A!203  has  been  shown  to  exhibit  a  superior  bal¬ 
listic  performance  compared  with  the  system  in  which  the 
two  phases  are  simply  uniformly  interdispersed.21,24  Micro¬ 
mechanical  simulations  have  also  demonstrated  the  effect  of 
microstructural  bias  on  failure  resistance.25,26  However,  the 
influence  of  microstructural  bias  on  the  fundamental  dy¬ 
namic  properties  of  these  ceramics  has  not  been  fully  estab¬ 
lished. 

The  objective  of  the  present  work  is  to  characterize  the 
high-strain-rate  deformation  and  damage  response  of  four 
types  of  microstructural Iy  biased,  two-phase  TiB2  +  A!203 
ceramics.  The  microstructural  bias  of  the  four  ceramics  falls 
into  two  morphological  categories.  The  first  category  in¬ 
volves  a  continuous  (interconnected)  TiB2  network  that  sur¬ 
rounds  A1203  (qualitatively  termed  “T@A”)  and  the  second 
category  involves  TiB2  and  A1203  phases  that  are  interdis- 
persed  and  uniformly  intertwined  with  each  other  (qualita¬ 
tively  termed  “Tin A”)-21  Normal  plate  impact  experiments 
were  used  to  measure  the  <rHEL  and  the  spall  strength  of  the 
materials.  The  measured  responses  were  then  correlated  with 
the  microstructure  morphological  characteristics  and  the 
phase  size  scales. 


11.  PREPARATION  AND  CHARACTERIZATION  OF  THE 
TWO-PHASE  CERAMICS 

The  two-phase  TiB2  +  A1203  ceramics  were  fabricated  as 
—80  mm  diameter  by  20-mm-thick  disks  by  hot-pressing 
powders.  The  powders  were  produced  through  either  self- 
propagating  high-temperature  synthesis  (SHS)  reactions  be¬ 
tween  powder  precursors  (3Ti02+ 3B203+  10Al=3TiB2 
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FIG.  1.  (a)  Strain-rate  dependence  of  compressive  strength  for  A1203 
+  TiB2  (solid  line)  in  contrast  to  that  for  Al203  (dashed  line)  (see  Ref.  23). 
(b)  A  comparison  of  the  dynamic  compressive  strengths  of  two-phase  TiB2 
+  A12Oj  with  those  of  the  monolithic  constituents  (sec  Ref.  23). 
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-FA1203;  33:67  weight  ratio),  or  through  manual  mixing 
(MM)  of  the  constituents  (nominally,  30  wt%  TiB2 
+  70wt%  A1203)  which  are  obtained  using  conventional 
powder  processing  techniques.  Following  the  SHS  synthesis 
or  manual  mixing,  the  powders  were  also  ball  milled  to  gen¬ 
erate  the  microstructural  bias.  Details  of  the  processing  ap¬ 
proaches  used  for  fabricating  these  two-phase  TiB2  +  Al203 
ceramics  are  described  elsewhere.21,27 

Quantitative  metallography28  was  performed  to  deter¬ 
mine  the  phase  size  (as  measured  by  the  linear  intercept 
length)  and  phase  connectivity  (based  on  measurement  of 
integral  curvature).  The  connectivity  of  a  phase  in  a  multi¬ 
phase  material  can  be  related  to  a  quantitative  measurement 
of  curvature,  k,  defined  as  the  angle,  dO ,  of  a  finite  segment 
of  an  arc  divided  by  the  length  of  the  arc  d\,  or 

k=d0fd\.  (]) 

Through  a  division  of  line  integrals,  the  average  curvature, 
*ave  can  be  calculated  by  dividing  the  angle  of  arc  fff  (in 
phase  A)  by  total  length  of  arc,  LA  .  Hence, 

^avc“  =  {  0A  —  0a}ILa  ,  (2) 
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FIG.  2.  Schematic  illustrations  of  the  setups  used  for  measurements  of  (a) 
stress  profiles  with  PVDF  gauges  and  (b)  frcc-surfacc  velocity  traces  with 
VISAR  interferometry. 


where,  0A  is  angle  for  a  convex  arc  and  0A  is  angle  for  a 
concave  arc  in  phase  A.  If  the  phases  are  considered  as 
closed  (as  in  the  case  of  our  two-phase  ceramic),  then  the 
angle  of  the  arc  is  2tt.  Hence  the  above  equation  can  be 
written  as 

kmc=  (2  ir)Nf/LA  =  2  tt{N+a  -  N~}ILa  ,  (3) 


TABLE  I.  Measured  average  values  of  integral  curvature,  and  sizes  of  con¬ 
stituent  phases. 


Microstructure 

type 

Average  integral 
curvature  (/xm-1) 

TiB2  phase 
size  {/xm) 

A1203  phase 
size  (/xm) 

A 

-  0.3 16±  0.022 

7.0 

10.4 

B 

-0.476±  0.046 

6.2 

9.1 

C 

-0.074±  0.028 

8.7 

25.1 

D 

-0.375  ±0.031 

7.9 

12.3 

where,  NAet,  NA  ,  and  NA  are  the  net,  positive,  and  negative 
numbers  of  loops  of  phase  A.  The  length  of  the  arc,  LA  ,  of  a 
closed  loop  is  the  perimeter  and  can  be  related  to  the  mean 
linear  intercept  length  L  using  La  —  {ttI2)(  ML).  The  average 
integral  curvature  of  phase  A  can  then  be  expressed  as 

k^=4'NnfL.  (4) 

In  the  present  work,  the  average  integral  curvature  val¬ 
ues  were  determined  using  Eq.  (4).  The  mean  linear  intercept 
length  L  and  the  number  of  net  loops  of  “A”  (or  TiB2)  phase 
around  the  “B”  (or  A1203,)  phase  were  determined  using 
IMAGEPRO-PLUS,  a  commercial  software  package  (product  of 
Media  Cybernetics).  To  account  for  anisotropy  effects,  the 
analysis  was  performed  on  at  least  ten  micrographs  of  both 
planar  and  cross-sectional  surfaces.  However,  the  differences 
due  to  anisotropy  were  found  to  be  in  the  range  of  the  stan¬ 
dard  deviation.29 

The  impact  experiments  used  an  80-mm-diam,  single- 
stage  gas  gun.  Measurements  of  the  <rHEL  and  the  shock 
wave  speeds  under  dynamic  compression  were  obtained 
from  stress  profiles  recorded  using  polyvinyldine  fioride 
(PVDF)  stress  gauges.  As  shown  in  Fig.  2(a),  a  TiB2  flyer 
plate  backed  by  an  air  gap  and  mounted  at  the  head  of  an 
aluminum  projectile,  was  used  to  impact  the  target  assembly 
consisting  of  — 3-mm-thick  ceramic  sample  backed  by  a 
TiB2  backer  plate.  Impact  experiments  were  conducted  at  a 
nominal  velocity  of —0.750  km/s  (—15  GPa  nominal  impact 
stress)  to  ensure  shock  loading  of  samples  under  similar  con¬ 
ditions.  The  impact  velocity  was  measured  using  arrival  time 
pins.  One  PVDF  gauge  package  was  placed  at  the  impact 
surface  and  another  was  placed  between  the  sample  and  TiB2 


Microstructure 


FIG.  3.  Optical  micrographs  of  two-phase  ceramics  (TiB2  white  and  A1203 
dark  phase),  (a)  Sample  A,  SHS,  T@A  microstructure,  (b)  sample  B,  SHS, 
TinA  microstructure,  (c)  sample  C,  MM,  T@A  microstructure,  and  (d) 
sample  D,  MM,  TinA  microstructure. 


FIG.  4.  Plot  comparing  the  average  integral  curvature  values  of  the  TiB2 
+  A1203  ceramic  samples  of  four  different  types  of  microstructures.  Note 
that  the  TiB2  phase  connectivity  is  greater  in  sample  A  than  in  sample  B,  and 
likewise  in  sample  C  than  in  sample  D. 
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TABLE  II.  Measured  values  of  density,  elastic  wave  velocity,  and  Young’s  and  shear  moduli  of  the  four  types 
of  microstructurally  biased  ceramic  samples. 

Sample 

type 

Density 
(g/cm3;  %TD) 

Long,  vcl 
(km/s) 

Shear  vcl. 
(km/s) 

Young’s 
mod  (GPa) 

Shear  mod 
(GPa) 

A 

4.095;  99.4% 

11.114 

6.477 

505.8 

171.8 

B 

3.939;  95.6% 

10.336 

6.145 

420.8 

148.7 

C 

4.108;  99.7% 

10.796 

6.527 

478.8 

175.0 

D 

4.083;  99.1% 

10.773 

6.523 

473.9 

173.7 

TiB2a 

4.509 

10.790 

7.43 

524.9 

248.9 

Al2<V 

3.55 

9.280 

5.47 

306.2 

106.4 

8 Values  of  TiB2  and  A!203  arc  from  D.  E.  Grady,  Sandia  Report  No.  SAND-91 -01 4,  UC-704,  1991. 


backer  plate.  The  experimental  setup  was  designed  such  that 
a  planar-parallel  shock  wave  propagates  through  the  target 
and  the  input  and  propagated  stress  profiles  are  measured 
with  little  interference  from  radial  reflected  waves.  The 
PVDF  gauges  were  connected  to  current  viewing  resistors  to 
allow  current  versus  time  data  acquisition  with  a  1  GHz  fre¬ 
quency  oscilloscope.  The  current  versus  time  data  were  then 
numerically  integrated  to  yield  stress  versus  time  profiles  for 
the  impact  and  propagated  wave  gauges.  The  arrival  time  at 
the  respective  gauges  was  taken  as  the  travel  time  through 
the  sample  thickness,  from  which  the  shock  wave  speed  was 
determined. 

Figure  2(b)  shows  the  experimental  setup,  illustrating  the 
target  sample  (without  any  backer)  being  impacted  by  a  pro¬ 
jectile  consisting  of  a  4140  steel  (or  silicon  carbide)  flyer 
plate  and  an  aluminum  sabot.  An  air  gap  exists  between  the 
flyer  plate  and  the  sabot  to  allow  full  unloading  required  for 
the  spall  strength  measurements.  All  sample  surfaces  were 
lapped  for  flatness  and  parallelism.  The  ceramic  targets  were 
polished  with  5  diamond  paste  to  ensure  reflectivity  re¬ 
quired  by  the  velocity  interferometer  system  for  any  reflector 
(VISAR)  beam.  The  time-resolved  longitudinal  motion  of  the 
sample  free  surface  was  measured  in  the  form  of  interference 
fringes  with  a  Valyn  VISAR  and  recorded  on  a  digital  oscil¬ 
loscope.  The  interference  fringes  were  then  converted  to 
time-resolved  history  of  particle  velocity. 

III.  RESULTS  AND  DISCUSSIONS 
A.  Microstructural  characterization 

Optical  micrographs  of  the  four  types  of  microstructur- 
ally  biased  samples  studied  are  shown  in  Figs.  3(a)~3(d). 
The  micrographs  correspond  to  sample  A  made  by  the  SHS 
process  with  TiB2  surrounding  A1203  (T@A);  sample  B  also 
made  by  SHS  with  intermixed  TiB2  and  A1203  (TinA) 
phases;  sample  C  made  by  MM  with  TiB2  surrounding 
A1203  (T@A);  and  sample  D  made  by  MM  with  intermixed 
TiB2  and  A1203  (TinA)  phases.  It  should  be  noted  that  the 
samples  do  not  reveal  100%  microstructural  bias,  i.e.,  while 
the  micrographs  of  samples  A  and  C  generally  illustrate  a 
nearly  continuous  (interconnected)  TiB2  phase  surrounding 
A1203 ,  regions  where  the  converse  is  true  are  also  present. 
Likewise,  samples  B  and  D  show  not  only  an  intermixed 
structure  but  also  the  continuous  phase  microstructures. 

Table  I  lists  the  values  of  average  integral  curvature  and 
the  average  size  of  TiB2  and  A1203  phases.  In  general,  the 


average  sizes  of  TiB2  and  AI203  phases  in  both  SHS  samples 
(A  and  B)  are  smaller  than  those  in  the  manually  mixed 
samples  (C  and  D),  with  sample  C  showing  the  largest  size 
for  both  constituents.  Sample  C  also  represents  the  micro¬ 
structure  with  the  highest  value  of  connectivity  for  TiB2 ,  as 
illustrated  in  Fig.  4.  All  samples,  in  fact,  show  a  negative 
curvature,  since,  within  a  given  area  of  measurement,  the 
number  of  loops  around  smaller  particles  is  greater  than  that 
surrounding  the  larger  particles.  This  skews  the  average  in¬ 
tegral  curvature  towards  negative  values.  Nevertheless,  the 
measured  values  provide  a  quantitative  measure  of  the  mi¬ 
crostructural  bias  on  the  basis  of  phase  connectivity,  and  il¬ 
lustrate  that  the  connectivity  of  TiB2  phase  is  greater  in 
sample  A  than  in  sample  B,  and  likewise  in  sample  C  than  in 
sample  D.  The  grain  size  of  the  respective  phases  also  shows 
a  similar  trend,  with  sample  A  being  coarser  than  sample  B, 
and  sample  C  being  coarser  than  sample  D. 

B.  Characterization  of  elastic  properties 

The  elastic  properties  of  the  ceramics  were  characterized 
using  an  Ultran  Laboratories  ultrasonic  test  apparatus.  The 
setup  allows  determination  of  the  longitudinal  and  shear 
wave  velocities  from  which  the  elastic  moduli  are  obtained 
using  densities  measured  on  individual  samples  by  the 
Archemedean  method.  Table  II  lists  the  results  of  elastic 
wave  velocity  measurements  and  corresponding  Young’s  and 
shear  moduli.  It  can  be  seen  that  the  moduli  of  the  various 
ceramic  samples  are  relatively  similar  except  in  the  case  of 
sample  B,  which  has  ~4.4%  porosity. 


FIG.  5.  Typical  input  and  backer  gauges  stress  profiles  (sample  B)  obtained 
from  Expt.  9916. 
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TABLE  III.  Summary  of  results  of  impact  experiments  under  dynamic  compression. 


Sample 
type  and 
expt.  No. 

Target 

thickness 

(mm) 

Density 

(%TMD) 

Wave 

speed 

(km/s) 

Hugoniot 
elastic  limit 
crHEL  (GPa) 

Tensile 
yield  str. 
o’ys  (GPa) 

A-9914 

3.03 

99.4% 

8.24±.83 

6.2±3.I 

4.2  ±2.1 

B-9916 

3.42 

95.6% 

9.67±  1.0 

4.4±  1.2 

Td- 

oo 

O 

+1 

ro 

C-9920 

3.36 

99.7% 

9.08  ±.74 

5.5±2.3 

4.02  ±1.7 

D-9921 

3.06 

99.1% 

8.31  ±.78 

8.5±4.5 

6.23  ±3.3 

C.  Measurements  of  dynamic  high-strain-rate 
mechanical  behavior 

Stress  wave  profiles  obtained  using  PVDF  stress  gauges 
were  used  to  determine  the  mechanical  behavior  under  dy¬ 
namic  compression.  VISAR  interferometry  was  used  to  ob¬ 
tain  free-surface  velocity  profiles  to  determine  the  tensile 
(spall)  strength. 

1.  Dynamic  compression  behavior 

Measurements  of  the  Hugoniot  elastic  limit  (o*HEL)  were 
obtained  from  normal  planar  impact  experiments  conducted 
at  an  impact  velocity  of  750  m/s.  Figure  5(b)  shows  ex¬ 
amples  of  stress  histories  recorded  by  the  “input”  and 
“backer”  PVDF  stress  gauges  from  experiment  9916,  for  the 
sample  B  ceramic.  Table  III  lists  the  sample  thickness,  den¬ 
sity,  impact  velocity  measured  using  shorting  pins,  wave 
speed  measured  by  considering  the  times  of  travel  through 
the  sample  thickness  as  recorded  by  input  and  backer  gauges, 
the  <thel,  and  the  yield  stress  in  simple  tension  (<jy s)  cal- 
culated  from  the  <rHEL. 

The  wave  speed  in  the  material,  which  is  a  function  of 
sample  density  and  micro  structure  in  addition  to  loading 
conditions,  is  measured  to  be  similar  for  all  four  samples 
(within  the  range  of  experimental  scatter).  The  Hugoniot 
elastic  limit  (trHEL),  the  axial  stress  at  which  a  solid  loaded 
under  conditions  of  uniaxial  strain  begins  to  exhibit  plastic 
deformation,  is  observed  to  be  a  strong  function  of  the  mi¬ 
crostructure  including  phase  size.  The  cr HEL  was  determined 
by  considering  one  half  the  value  between  the  peak  and 
trough  of  the  elastic  precursor  wave  [Fig.  5(b)],  and  one  half 
the  difference  was  considered  as  the  range.  The  <rHEL 
found  to  be  the  lowest  for  sample  B  (4.4±  1.2  GPa),  due  to 
its  high  level  of  porosity  (—4%)  and  highest  for  sample  D 
(8.5±4.5GPa)  (with  a  large  standard  deviation).  The  cr HEL 
values  reported  previously  by  Grady4,6  are  9-18  GPa  for 
TiB2 ,  and  —6.7  for  A1203 .  While  crHEL  identifies  the  limit  of 
elastic  response  under  dynamic  (shock)  uniaxial  strain  load¬ 
ing,  it  is  also  common  to  assume  a  Von  Misses-type  yield 
condition  which  asserts  that  plastic  flow  initiates  when  the 
second  deviatoric  stress  invariant  attains  a  critical  value. 
Through  this  formalism,  the  yield  stress  in  simple  tension 
has  been  shown  to  be  related  to  (crHEL)  through  crYs 
=  2(C^/C^)X (crHEL).6  The  yield  stress  in  simple  tension, 
cr  ys  for  the  various  micro  structurally  biased  samples  follows 
the  same  trend  as  the  Hugoniot  elastic  limit,  with  sample  D 
showing  the  highest  and  sample  B  the  lowest  yield  strength. 


2.  Tensiie  (spall)  strength 

Spall  experiments  were  perfonned  (with  VISAR  interfer¬ 
ometry)  on  samples  of  microstructure  A  at  impact  velocities 
corresponding  to  input  stresses  above  and  below  the  crHEL, 
and  on  samples  of  microstructures  ‘C’  and  ‘D’  under  elastic 
loading  conditions  to  ensure  that  the  compression-induced 
damage  does  not  influence  the  tensile  response.  Figure  6 
shows  the  free-surface  velocity  traces  obtained  on  samples  of 
microstructure  A  at  237,  495,  and  758  m/s,  using  a  2.66-mm- 


b)  Time(ps) 


FIG.  6.  Free-surface  velocity  traces  obtained  on  samples  of  Microstructure 
A  at  237(a),  495(b),  and  758(c)  m/s,  using  a  2.66-mm-thick  AISI  4140  steel 
flyer  plate.  Arrows  indicate  spall  signal. 


Downloaded  30  Oct  2002  to  128.61.137.209.  Redistribution  subject  to  AIP  license  or  copyright,  see  http://ojps.aip.org/japo/japcr.jsp 


1926  J.  Appl.  Phys.,  Vol.  91,  No.  4,  15  February  2002 


Kennedy  et  at. 


TABLE  IV.  Loading  conditions  and  summary  of  results  of  spall  experiments. 


Sample 

Flyer 

Flyer  & 

Target 

Impact 

Input 

Spall 

type  and 

plate 

target  th. 

density 

velocity 

stress 

strength 

expt.  No. 

material 

(mm) 

(%TMD) 

(m/s) 

(GPa) 

(GPa) 

A-0005 

4140  steel 

2.66/16.22 

98.8% 

495 

7.9 

0.320 

A-0007 

4140  steel 

2.66/16.20 

98.8% 

237 

3.7 

0.320 

A-000S 

4140  steel 

2.66/7.99 

98.8% 

758 

11.8 

0.160 

C-9925 

SiC 

4.30/7.54 

99.7% 

244 

3.8 

0.311 

D-9926 

SiC 

4.77/7.10 

99.1% 

239 

3.8 

0.222 

thick  AISI  4140  steel  flyer  plate.  The  velocity  traces  reveal  a 
spall  signal  (velocity  decrease,  AC7fs)  which  was  used  to 
compute  the  spall  strength  (<rspal,=  l/2p0C0At/fs.)  As 
shown  in  Table  IV,  the  spall  strength  for  sample  A  is  0.320 
MPa  at  input  stresses  of  3.7  and  7.9  GPa.  However,  with 
input  stress  increasing  to  11.8  GPa,  the  spall  strength  de¬ 
creases  to  0.160  GPa.  It  can  be  seen  that  while  the  spall 
strength  decreases  with  increasing  input  stress,  the  two-phase 
TiB2-fAl203  ceramic  maintains  non-ncgligible  spall 
strength  even  at  input  stresses  exceeding  the  Hugoniot  elastic 
limit  (6.2±  3.4  GPa). 

Tensile  spall  experiments  were  also  performed  on 
samples  of  microstructures  C  and  D,  at  an  input  stress  of 
~~3.8  GPa.  Figure  7  shows  the  free  velocity  traces  for  these 
samples.  A  spall  strength  of  0.311  GPa  for  sample  C  and 
0.222  GPa  for  sample  D  was  observed.  A  comparison  of 
spall  test  results  between  sample  C  and  sample  D  shows  that 
the  latter  sample  with  dispersed  microstructurc  has  a  lower 
spall  strength,  while  the  microstructure  with  interconnected 
TiB2  has  a  higher  spall  strength,  even  though  sample  D  has  a 
smaller  phase  size.  The  measured  high  value  of  spall  strength 
of  the  two-phase  ceramic  is  similar  to  the  published  spall 
strength  of  TiB2  (0.33  GPa)  but  lower  than  that  of  A!203 
(0.45  GPa). 

IV.  DISCUSSION  AND  SUMMARY  OF  RESULTS 

The  two-phase  TiB2  +  Al203  ceramics,  made  either  by 
the  SHS  or  mechanical  milling  methods,  reveal  differences 
in  microstructure  which  qualitatively  show  TiB2  as  a  con¬ 
tinuous  (interconnected)  phase  surrounding  A1203  (T@A), 
or  TiB2  and  A1203  intermixed  with  each  other  (TinA).  Quan- 


FIG.  7.  Free- surface  velocity  traces  obtained  for  samples  C  and  D,  at  similar 
loading  conditions. 


titative  microscopy  analysis  based  on  the  measurement  of  the 
integral  curvature  showed  that  the  samples  investigated  do 
not  exhibit  100%  microstructura!  bias.  However,  the  overall 
trend  of  the  influence  of  microstructural  bias  emerging  from 
the  results  of  experiments  performed  to  date  illustrates  that 
the  dynamic  yield  strength  and  the  crHEL  are  more  domi¬ 
nantly  dependent  on  the  phase  size.  Sample  C  prepared  by 
manual  mixing  and  having  the  largest  phase  (grain)  size 
shows  the  lowest  values  in  contrast  to  the  other  samples  of 
similar  (~99%)  density.  In  contrast,  the  tensile  spall  strength 
appears  to  scale  with  the  continuity  of  the  TiB2  phase. 
Sample  C,  which  has  the  most  interconnected  TiB2  phase, 
has  the  highest  value  of  the  tensile  spall  strength. 

The  results  therefore,  illustrate  that  while  the  Hugoniot 
elastic  limit  and  the  dynamic  compressive  yield  strength  of 
Al203H“TiB2  are  dependent  on  the  average  grain  (phase) 
size,  the  tensile  spall  strength  scales  with  the  TiB2-phase 
connectivity,  and  less  so  with  the  average  phase  size.  It  is 
possible  that  the  interconnected  phase  morphology  is  more 
effective  in  impeding  initiation  and  progression  of  fracture 
under  tensile  conditions.  If  so,  TiB2  as  an  interconnected 
phase  has  the  ability  to  suppress  and  relax  the  cracks  formed 
in  A1203 .  Alternatively,  it  is  also  possible  that  in  the  two- 
phase  ceramic  with  the  microstructure  containing  dispersed 
A1203  and  TiB2  phases,  the  mechanical  energy  trapping  due 
to  scattering  of  waves  from  incoherent  boundaries  and  inter¬ 
faces  results  in  the  lowering  of  the  tensile  (spall)  strength. 
Further  work  is  currently  in  progress  to  more  clearly  delin¬ 
eate  the  effects  of  interconnected  versus  dispersed  TiB2 
phase,  and  to  eventually  fabricate  mo-phase  Al203  +  TiB2 
ceramics  with  the  microstructural  bias  that  yields  the  most 
optimal  dynamic  properties. 
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Abstract 

Microstructurally  biased  two-phase  TiB2+Al203  ceramic  with  an  interconnected  or  dispersed  TiB2 
(minor)  phase  can  be  produced  by  variations  in  processing  parameters.  A  standard  method  for 
quantitatively  characterizing  the  different  types  of  microstructures  is  necessary  to  comprehend 
mechamsm(s)  that  control  the  generation  of  the  microstructural  bias  during  processing  and  consequently  the 
resultant  properties.  In  this  work,  techniques  derived  from  stereology  were  used  to  quantitatively 
characterize  the  microstructure  based  on  phase  size,  and  phase  morphology  in  terms  of  phase  connectivity 
and  dispersion  of  the  minor  phase  (TiB^.  Mean  integral  curvature  based  on  area  particle  count  and  area 
tangent  count  methods  were  used  to  quantitatively  describe  the  connectivity  of  the  TiB2  phase  around  the 
Al203  major  phase.  The  results  illustrate  that  in  spite  of  partial  and  mixed  bias,  integral  curvature 
measurements  (particularly  those  based  on  area  tangent  count  method)  provide  an  accurate  means  for 
quantitative  characterization  of  the  microstructure  and  its  correlation  with  measured  properties. 
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1.  Introduction 

Most  traditional  or  advanced  ceramics  contain  more  than  a  single  constituent  phase.  The  second 
phase  may  either  be  intentionally  added  or  get  picked  up  during  processing.  For  example,  silicon  carbide 
[1],  titanium  carbide  [2],  and  boron  carbide  [3]  have  previously  been  added  to  an  alumina  matrix  to  inhibit 
the  growth  of  cracks  by  providing  a  plastic  “energy-absorbing”  barrier  and  thereby  enhancing  material 
properties.  The  addition  of  a  second  phase  can  also  assist  in  ceramic  fabrication.  For  example,  sintering 
aids  such  as  B4C  assist  in  the  pressure-less  sintering  of  TiB2  ceramics  resulting  in  enhanced  density  and  grain 
growth  suppression  [4].  On  the  other  hand,  impurities  such  as  Fe  (~  0.5 wt%)  cause  abnormal  grain  growth 
and  also  yield  a  low  sintered  density.  Likewise,  presence  of  an  unintentionally  formed  glassy  second  phase 
can  weaken  the  ceramic  and  deteriorate  its  properties. 

Efforts  to  design  ceramics  with  unique  microstructures  that  provide  improved  material  properties 
are  receiving  great  interest.  Multiphase  ceramics  are  being  developed  to  subside  the  inherent  brittle  nature 
of  these  materials  and  attain  a  combination  of  high  strength  and  toughness,  and  wear  and  impact  resistance. 
For  these  materials,  while  it  may  be  required  to  have  a  certain  desired  distribution  of  the  respective  phases, 
the  control  in  forming  the  complex  microstructure  is  often  limited.  The  properties  of  two-phase  ceramics 
are  influenced  by  microstructural  characteristics,  including  type  of  phases  present,  their  size,  morphology, 
and  chemical  composition.  Mechanical  properties  establish  material  limitations  for  specific  engineering 
applications.  Many  of  these  properties  are  dictated  not  only  by  the  intrinsic  properties  of  the  constituent 
phases,  but  also  the  interaction  between  the  phases.  Therefore,  it  is  important  to  correlate  processing 
conditions  and  constituent  phase  properties,  with  final  microstructure  and  the  bulk  properties  of  engineered 
ceramics. 

In  recent  years,  several  studies  have  been  conducted  to  produce  two-phase  TiB2+ Al203  ceramics 
(nominally  30%  TiB2  and  70%  Al203  by  weight)  for  a  variety  of  engineering  applications  [5-10],  These  have 
been  fabricated  by  manual  mixing  (MM)  of  commercially  acquired  TiB2  and  Al203  powders  and  also  by 
using  thermite-type  self-propagating  high-temperature  synthesis  (SHS)  reactions  in  precursor  powders  of 
TiOz,  B203,  and  Al  [19,20].  Logan  [5]  and  Carney  [6]  have  synthesized  microstructurally  biased  two-phase 
TiB2+Al203  ceramics  with  an  interconnected  or  dispersed  TiB2  (minor)  phase  by  simple  variations  in 
processing  parameters.  Micrographs  of  samples  shown  in  figure  l(a,b)  were  produced  by  SHS  formed 
powder  after  ball-milling  and  finally  hot-pressing.  Micrographs  of  samples  shown  in  figure  l(c,d)  were 
produced  by  manually  mixing  commercially-acquired  TiB2  and  Al203  powders  and  subsequent  hot-pressing 
to  form  the  dense  ceramic  disks.  Microstructures  labeled  “A”  (in  figure  1(a))  and  “C”  (in  figure  1(c)) 
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comprising  of  a  relatively  continuous  interconnected  TiB2  network  surrounding  an  Al203  phase  were 
qualitatively  termed  “T@A”.  Microstructures  labeled  “B”  (in  figure  1(b))  and  “D”  (in  figure  1(d)),  with  the 
TiB2  phase  interdispersed  within  a  continuous  network  of  Al203  were  qualitatively  termed  “TinA”. 
Microstructure  B  reveals  a  TiB2  phase  that  is  qualitatively  more  dispersed  (“TinA”)  in  Al203,  although, 
regions  with  “T@A”  are  also  observed.  Furthermore,  the  size  of  the  two  phases  is  much  smaller  in 
microstructure  B  than  in  microstructure  A .  It  has  been  shown  that  these  two-phase  ceramics  exhibit 
mechanical  properties  that  are  not  necessarily  based  on  the  rule  of  mixtures  [12],  with  their  compressive 
strengths  being  higher  than  either  of  the  constituents,  both  at  low  and  high  strain  rates.  It  has  also  been 
observed  that  while  the  compressive  yield  strength  at  high  strain  rates,  scales  with  the  grain  size  (increasing 
with  decreasing  grain  size),  the  tensile  spall  strength  increases  with  increase  in  connectivity  of  the  TiB2  phase 
[13,14].  The  dissimilar  influence  of  phase  size  and  microstructural  bias  on  the  static  and  dynamic 
mechanical  properties  [11-14]  of  these  two-phase  TiB2+Al203  ceramics  is  further  complicated  by  the  fact 
that  the  microstructural  bias  revealed  in  a  given  sample  is  not  100%,  particularly  in  samples  made  by  the 
SHS  process.  It  has  also  not  been  established  as  to  which  stage(s)  of  processing  most  dominantly  influences 
the  generation  of  microstructural  bias  in  terms  of  phase  connectivity  or  dispersion,  starting  from  pre-mixing 
of  constituent  powders  of  specific  morphology,  SHS  reaction  process,  ball-milling  of  reaction  products,  or 
hot-pressing  parameters  used  during  powder  densification. 

The  objective  of  the  present  work  was  to  quantitatively  characterize  the  biased  microstructure  of 
two-phase  TiB2+ Al203  ceramics  to  better  understand  the  evolution  of  microstructure  and  control  of 
processing  mechanism(s).  The  approach  involved  quantitative  characterization  of  the  distribution  and 
morphology  of  the  respective  TiB2  and  Al203  phases,  followed  by  the  correlation  of  these  characteristics 
with  the  various  stages  of  processing  of  the  two-phase  ceramic.  In  the  following  section,  the  experimental 
procedure  employed  for  processing  the  ceramics,  microstructure  analysis,  and  quantitative  characterization 
will  be  described  first,  followed  by  the  presentation  of  observed  results  and  their  discussion. 


2.  Experimental  Procedure 

2.1  Processing  of  Two-Phase  TiB2 + A1203  Ceramics 

Processing  of  two-phase  TiB2 + A1203  ceramics  (30/70  nominal  weight  percent  ratio)  was  performed 
either  by  manually  mixing  TiB2  and  Al2Q3  powders  (referred  to  as  MM  powders),  or  by  using  SHS  reaction 
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(referred  to  as  SHS  powders)  based  on  the  following  reaction:  3Ti02  +  3B203  +  10A1  — >  3TiB2  +  5Al203. 
The  MM  as  well  as  SHS  synthesized  powders  were  ball-milled  for  various  times  and  subsequently  hot- 
pressed  into  98%  dense  disks,  approximately  7.62mm  diameter  by  1.27mm  thick  (3”  x  V2”),  and  containing 
different  types  of  biased  microstructures,  similar  to  those  described  in  previous  studies  [5-7].  Variations  in 
ball-milling  times,  and  the  pressure  and  temperature  during  hot-pressing  were  used  to  produce  the 
microstructural  bias  generated  as  a  result  of  the  complex  interactions  between  the  various  constituent 
particles  during  sintering  under  pressure. 

The  hot-pressed  TiB2+ Al203  ceramics  exhibit  a  microstructural  bias,  which  shows  TiB2  either  as  an 
interconnected  phase  surrounding  the  Al203  phase,  or  the  two  phases  being  simply  dispersed,  as  illustrated 
in  figure  l(a-d).  Quantitative  characterization  using  stereological  methods  was  employed  to  quantitatively 
describe  the  relevant  features  of  the  microstructures.  Features  of  interest  included  phase  size,  interface 
curvature,  and  phase  connectivity.  Microstructural  characterization  was  performed  on  polished  two- 
dimensional  sections  (surfaces  perpendicular  and  parallel  to  pressing  direction)  of  samples  cut  from  the  bulk 
disks  of  hot-pressed  SHS  and  MM  compacts  made  for  this  study.  Characterization  was  also  performed  on 
hot-pressed  samples  with  microstructures  A  through  D  (surfaces  perpendicular  to  pressing  direction  only) 
prepared  in  previous  work  by  Carney  [6].  Micrographs  from  planar  and.  cross-sectioned  surfaces  were 
obtained  by  imaging  the  samples  using  optical  microscopy.  Quantitative  stereology  techniques  were  used  to 
determine  the  mean  integral  curvature  (using  both  area  particle  and  area  tangent  count  methods)  and  phase 
size.  The  measurements  were  performed  on  a  total  of  four  sections,  two  planar  top  and  bottom  surfaces 
and  two  cross-sectional  surfaces  of  opposite  orientations,  cut  from  the  hot-pressed  disk  shaped  samples. 
Measurements  were  made  on  flat  sections  of  the  material  and  then  transformed  from  two-dimensional  to 
three-dimensional  approximations  of  the  geometric  relations  between  the  different  features. 


2.2  Quantitative  Stereological  Characterization 

An  automated  quantitative  characterization  approach  was  devised  for  the  measurement  of  mean 
integral  curvature  using  both  area  tangent  and  area  particle  count  methods  [15-18].  Both  techniques 
measure  the  degree  of  connectivity  of  the  TiB2  (minor)  phase  present  in  the  two-phase  ceramic,  and  are 
described  in  more  detail  in  the  Appendix.  The  measure  of  average  curvature  according  to  the  particle  count 
method  is  given  by: 
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where  NA  +  is  the  number  of  positive  loops  per  unit  area  (or  TiB2  surrounding  A1203),  NA_  is  the  number  of 
negative  loops  per  unit  area  (or  TiB2  dispersed  in  Al203  matrix),  and  the  total  length  of  the  loops  per  unit 
area,  LA  is  actually  the  perimeter  of  the  closed  loops  formed  by  the  phase  boundary,  which  is  easily  related 

to  the  mean  linear  intercept  length,  l .  The  average  curvature  according  to  the  area  tangent  count  method  is 
given  by: 


(2) 


where,  TA+  and  TA.  are  tangents  of  convex  and  concave  surfaces  of  Al203  holes  and  NL  is  the  average 
number  of  intersections  of  phase  boundaries  with  a  test  line  per  unit  length. 

Physical  observations  and  quantitative  characterization  of  the  microstructures  were  performed  on  a 
Leica  DM  IRM  (Bannockburn,  IL)  optical  microscope  operated  in  a  bright  field  mode  with  a  20x  objective 
lens.  Digital  images  were  obtained  using  a  CoolSNAP  CCD  video  camera  (Media  Cybernetics;  Silver 
Springs,  MD),  capable  of  capturing  high-resolution  (1390  x  1040  pixels,  4.65  micron  pixels)  digital  images. 
Images  were  imported  into  Image-Pro  Plus  (Media  Cybernetics;  Silver  Springs,  MD)  analytical  imaging 
software  where  the  number  of  features,  volume  fraction,  and  average  phase  size,  etc.,  were  measured. 

Calculations  of  mean  integral  curvature,  phase  size,  and  phase  distributions  were  performed  using 
spatially  calibrated  images.  Direct  measurements  and  “real”  values  could  be  taken  directly  from  the  images. 
Captured  images  were  1390  x  1040  pixels,  which  correspond  to  a  0.2325pm/pixel  resolution  for  a  20x 
objective  lens  with  the  smallest  resolvable  measured  feature  greater  than  or  equal  to  0.70pm.  Measurements 
were  taken  on  five  randomly  selected  regions  for  each  of  the  planar  and  cross-sectional  specimen  sections. 

Curvature  and  phase  size  calculations  were  accomplished  by  making  “basic”  measurements  on  each 
of  the  captured  images.  Basic  measurements  included  number  of  features  or  phases,  average  feature  size, 
mean  liner  intercept,  tangent  count,  and  phase  volume  fraction.  Light  contrast  features  (TiB^  were  first 
separately  counted  and  then  dark  contrast  features  (Al203)  were  counted  for  the  same  image. 

Volume  fraction  and  mean  linear  intercept  [18]  were  automatically  calculated  using  Materials-Pro 
plug-in  (Media  Cybernetics;  Silver  Springs,  MD)  software  for  Image-Pro  Plus.  Volume  fraction  was 
calculated  according  to  pixel  area  fraction  for  each  phase  of  the  material.  Automatic  thresholding  was 
performed  for  each  image  before  any  measurements  were  taken.  Volume  fraction  results  for  respective 
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phases  were  used  to  check  the  automatic  thresholding  feature.  Typically,  automatic  thresholding  results 
were  close  to  the  nominal  30%  TiB2  and  70%  Al203  phase  distribution.  Mean  linear  intercept  values  were 
measured  using  127  horizontal  test  lines.  The  software  reported  the  total  line  length  and  number  of  phase 
intercept  counts. 

The  area  tangent  count,  used  for  calculating  curvature,  was  measured  using  The  Image  Processing 
Tool  Kit  3.0  software  developed  by  J.  C.  Russ  (Reindeer  Games,  Inc.;  Asheville,  NC).  The  software 
measures  the  number  of  positive  and  negative  tangents  formed  with  the  phase  boundary  and  reports  the 
two  values  separately. 


3.  Results  of  Quantitative  Microstructural  Characterization 

The  ability  to  quantitatively  describe  two-phase  ceramics  requires  the  consideration  of  geometrical 
and  topological  features  such  as  number,  size,  and  shape  of  phases  and  the  continuity  or  connectivity  of  the 
constituent  phases.  Topological  properties  are  unchanged  when  the  structure  is  stretched  or  distorted, 
which  has  to  do  with  the  connectivity  of  a  structure  [15].  A  multi-phase  structure  becomes  more 
interconnected  as  the  numbers  of  separate  features  (from  one  phase)  decrease,  and  impinge  upon  each  other 
forming  a  system  of  closed  continuous  loops.  The  total  connectivity  of  a  structure  can  be  shown  to  equal 
the  number  of  features  it  contains  and  the  value  of  the  net  angle  subtended  by  the  closed  curve  forming 
around  the  feature  (equal  to  2ti  radians  regardless  of  shape).  Thus,  as  described  in  more  detail  in  Appendix 
-A,  the  net  subtended  angle  is  a  “topological  invariant”  for  a  closed  curve  and  is  related  to  the  measure  of 
integral  curvature  [15]. 

Microstructures  A  through  D  fabricated  from  previous  work  [6]  and  shown  earlier  in  figure  1  (a)-(d), 
were  examined  first.  Microstructures  A  and  C  qualitatively  show  the  most  connective  TiB2  phase  for  the 
SHS  and  MM  samples  respectively  as  shown  in  figure  l(a,c).  Microstructure  D  has  qualitatively  the  most 
dispersed  TiB2  phase,  while  Microstructure  B  has  a  combination  of  features  from  microstructures  A  and  D. 
Sample  A  has  several  very  small  TiB2  particles  (<  1pm)  and  also  areas  where  the  TiB2  phase  is  very 
connective.  Quantitative  measurements  of  mean  integral  curvature  indicated  the  same  trend.  Figure  2(a) 
shows  the  values  obtained  using  both  curvature  calculation  methods  (area  particle  count  and  area  tangent 
count),  calculated  according  to  equations  (1)  and  (2),  respectively.  The  structure  becomes  more  connective 
as  the  curvature  value  approaches  zero.  Typically,  area  particle  count  and  area  tangent  count  revealed  the 
same  basic  trend  with  differences  in  magnitude  ranging  from  approximately  10  to  20%  between  the  two 
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methods  (for  measurements  performed  on  five  metallographic  locations  for  each  sample).  The  area  tangent 
count  values  were  however,  typically  slightly  greater  than  the  area  particle  count  values  as  depicted  in  figure 
2(a).  Microstructure  C  has  the  most  connective  phase,  indicated  by  area  tangent  count  curvature 

approaching  zero  ( k  =  -0.12pm1)  while  microstructure  D  has  the  most  dispersed  TiB2  phase 
(k  =  -0.46pm'1).  Although  SHS  formed  microstructures  exhibit  both  connective  and  dispersed  TiB2 
regions,  microstructure  A  revealed  a  more  connective  TiB2  phase  (it  =  -0.40pm'1)  in  contrast  to 
microstructure  B  (fc  =  -0.64pm'1). 

Compacts  were  also  made  in  the  present  work  to  duplicate  the  microstructures  of  samples  A,  B,  C, 
and  D,  produced  in  prior  work.  Samples  A  and  B  were  reproduced  by  ball-milling  SHS  formed  material  for 
4  and  30  hours,  respectively,  while  samples  C  and  D  were  reproduced  by  using  MM  technique  with  no  ball¬ 
milling  and  2  hour  ball-milling  times,  respectively.  Figure  2(b)  shows  a  plot  of  the  curvature  values  (based 
on  area  tangent  count  only)  measured  for  samples  yl  through  D  made  in  prior  work  and  representative 
samples  made  from  the  present  work.  Good  agreement  of  the  overall  trends  among  the  curvature  values 
was  observed  between  the  two  groups  of  samples. 

Additional  measurements  were  performed  on  SHS  formed  samples  that  were  subjected  to  various 
ball-milling  periods  (from  1  to  30  hours)  before  hot-pressing.  Optical  micrographs  in  figure  3(a-g)  show 
differences  in  the  microstructural  bias  following  ball-milling  for  1,  2,  4,  8,  15,  23,  and  30  hours.  Each 
micrograph  exhibits  regions  that  have  both  interconnected  and  dispersed  TiB2  phases.  It  can  be  seen  that 
the  microstructures  of  samples  ball-milled  for  1, 2,  and  4  hours  (figure  3(a-c))  appear  to  reveal  a  qualitatively 
greater  degree  of  TiB2  phase  connectivity  surrounding  Al203.  However,  as  the  ball-miliing  time  increases  (8 
to  30  hours),  the  bias  of  TiB2  around  Al203  gradually  shifts  to  a  more  dispersed  microstructure  as  shown  in 
figure  3(d-g).  Mean  integral  curvature  measured  for  all  SHS  prepared  samples  ball-milled  for  1  to  30  hours 
is  shown  in  figure  4.  Samples  ball-milled  for  less  than  4  hours,  have  the  most  connective  TiB2  phase  around 
A1203,  while  the  30  hour  milled  samples  had  a  curvature  value  corresponding  to  the  least  interconnected  or 
uniformly  dispersed  TiB2  in  Al203  phase. 

The  effect  of  ball-milling  on  MM  prepared  samples  was  also  determined.  Figure  5(a,b)  shows 
micrographs  of  MM  prepared  hot-pressed  samples  with  zero  and  2  hours  of  ball-milling.  The  zero  hour 
ball-milled  sample  shows  the  “T@A”  microstructure,  with  continuously  interconnected  bright-contrast  TiB2 
surrounding  the  large  Al203  phase,  while  the  2-hour  ball-milled  sample  shows  the  “TinA”  microstructure 
with  TiB2  interdispersed  in  Al203.  Hence,  ball-milling  not  only  decreases  the  size  of  the  constituent  phases, 
but  also  drastically  influences  the  resulting  microstructural  bias,  more  so  in  MM  samples  than  in  SHS 
samples.  Quantitative  analysis  also  measured  a  2  hour  milled  SHS  sample  having  a  mean  integral  curvature 
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of  -0.24pm1  and  a  MM  sample,  following  the  same  ball-milling  time,  showed  a  curvature  of  -0.32pm1. 
Microstructures  of  corresponding  SHS  and  MM  samples  on  which  the  curvature  and  phase  size  values  were 
measured  are  shown  in  figures  6(a)  and  (b),  respectively. 

The  phase  size  (calculated  using  the  mean  linear  intercept  method  described  in  Appendix  A)  for 
samples:  A,  B,  C,  and  D  from  prior  work  (figure  7)  and  all  of  the  SHS  hot-pressed  samples  made  in  the 
present  study  shown  in  figure  8.  Typically,  the  TiB2  phase  for  all  samples  was  smaller  than  that  of  the  Al203 
phase.  As  illustrated  in  figure  7,  the  average  size  of  respective  phases  in  SHS  samples  A  and  B  is  not  very 
different  although,  the  Al203  size  distribution  in  sample  A  is  greater  then  in  sample  B.  In  contrast,  for  the 
MM  sample,  the  size  and  size  distribution  of  both  phases  (TiB2  and  Al203)  in  sample  C  is  greater  then  in 
sample  D.  The  influence  of  powder  ball-milling  time  on  the  final  phase  size  of  the  hot-pressed  samples 
shows  that  the  size  of  both  TiB2  and  Al203  phases  decreases  as  the  powder  is  milled  for  longer  periods  of 
time  as  illustrated  in  figure  8.  The  1  hour  milled  sample  has  a  TiB2  and  Al203  phase  size  of  6.3  and  12.7pm, 
respectively,  while  the  30  hour  milled  sample  shows  a  TiB2  phase  size  of  4.5pm  and  9.3pm  for  Al203.  Also, 
typically  a  narrower  phase  size  distribution  is  achieved  with  increasing  ball-milling  time.  Sampling  different 
orientations  did  not  appear  to  have  a  significant  effect  on  the  measured  phase  size  values.  As  an  example, 
the  values  of  average  phase  sizes  for  TiB2  and  Al203  from  each  sample  section  of  an  SHS  sample  range 
between  5.6  pm  and  6.3  pm  for  TiB2  and  11.6  pm  and  12.5  pm  for  Al203. 

To  determine  the  effects  of  hot-pressing  variables,  including  applied  pressure,  temperature,  and  soak 
time  at  temperature,  each  was  separately  evaluated  by  measuring  phase  size  and  mean  integral  curvature  for 
the  SHS  samples.  Typically,  the  density  of  each  sample  regardless  of  the  hot-pressing  method  resulted  in  a 
variation  less  than  1%. 

The  phase  size  measurements  were  examined  for  samples  to  see  the  effect  of  soak  time  during  hot 
pressing.  Sample  SHS-025  and  SHS-021  was  soaked  for  150  and  240  minutes,  respectively,  during  hot- 
pressing  at  a  temperature  of  1620°C.  Correspondingly,  the  TiB2  and  Al203  phase  sizes  were  5.8  and  9.3pm 
for  SHS-025  and  4.8  and  10.1pm  for  SHS-021,  illustrating  a  minor  effect  of  soak  time  on  phase  size. 
Pressure  application  did  not  appear  to  affect  the  final  phase  size  for  each  constituent  phase.  Samples  SHS- 
033  and  SHS-040  were  virtually  identical  for  constant  application  of  pressure  versus  applying  the  pressure 
when  the  powder  compact  reaches  sintering  temperature.  Sample  SHS-033  (34.5MPa  constant  application) 
and  SHS-040  (34.5MPa  applied  once  reaching  1620°C)  had  a  phase  size  of  5.1  and  5.5  pm  for  the  TiB2 
phase  respectively.  Similarly,  the  Al203  phase  had  sizes  of  10.4  and  10.6  pm  for  the  two  samples.  Changing 
the  sintering  temperature  also  did  not  show  any  effect  on  the  overall  phase  size.  Samples  SHS-028  and 
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SHS-023  were  sintered  at  1620°  and  1700°C,  respectively,  and  differences  between  the  two  samples 
measured  phase  sizes  were  5.0%  for  TiB2  and  7.7%  for  Al203. 

In  addition  to  evaluating  hot-pressing  conditions  affecting  the  constituent  phase  sizes,  the  final 
microstructure  morphologies  (bias)  were  also  considered.  Sample  SHS-025  was  ball-milled  for  30  hours  and 
had  an  average  powder  particle  size  of  6.0  pm.  This  sample  was  held  at  temperature  for  150  minutes  during 
hot-pressing.  Sample  SHS-021  was  also  ball-milled  for  30  hours  and  had  a  similar  average  powder  particle 
size  of  5.8  pm,  and  was  held  at  temperature  for  a  longer  time  period  of  240  minutes.  Mean  integral 
curvature  for  SHS-021  and  SHS-025  samples  were  -0.50  and  -0.38pm"1  respectively.  The  sample  held  at 
temperature  for  a  longer  time  period  (SHS-021)  had  a  more  dispersed  TiB2  phase  based  on  mean  integral 
curvature  values.  All  other  parameters  for  processing  the  two  samples  were  similar. 

Additionally,  no  detectable  differences  were  noticed  between  the  samples  that  were  formed  using 
alternate  pressure  applications  and  sintering  temperatures  during  processing.  Samples  SHS-033  (constant 
pressure  of  34.5MPa)  and  SHS-040  (pressure  of  34.5MPa  once  temperature  is  reached)  had  mean  integral 
curvature  values  of  -0.40  and  -0.38pm 1  respectively  and  the  pressure  application  method  did  not  appear  to 
affect  the  curvature  value.  Similarly,  SHS-037  and  SHS-023  were  sintered  at  1620°  and  1700°C  and  mean 
integral  curvature  values  ranged  from  -0.44  and  -0.47pm"1  respectively.  Thus,  the  variation  of  hot-pressing 
parameters  measured  by  curvature  appears  to  vary  more  with  the  starting  powder  particle  size  as  a  function 
of  ball-milling  time  than  the  effects  of  sintering  temperature  and  pressure  application. 


4.  Discussion  of  Results 


The  main  objectives  of  this  research  were  to  quantitatively  characterize  the  biased  microstructure  of 
two-phase  TiB2+Al203  ceramics  and  understand  the  effects  of  processing  variable  that  control  the  phase 
development  and  microstructural  bias.  Quantitative  characterization  of  phase  connectivity  in  each 
microstructure  was  achieved  by  measuring  the  mean  integral  using  both  area  particle  count  and  area  tangent 
count  methods.  The  measure  of  mean  integral  curvature  simply  involves  counting  positive  and  negative 
loops  when  the  microstructure  is  composed  of  isolated,  convex  features  as  seen  in  microstructure  D  (figure 
1(d)).  However,  it  becomes  more  challenging  to  measure  a  microstructure  similar  to  microstructure  A 
(figure  1  (a)),  which  comprises  of  complex  interactions  between  the  two  phases.  Phases  intersecting  the  edge 
of  the  micrograph  are  not  counted  and  provoke  measurement  errors.  However,  the  area  tangent  count 
method  does  not  suffer  from  this  edge  effect  problem  because  a  test  line  is  swept  across  the  entire  field  of 
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view  and  all  features  are  counted  in  the  measurement.  Also,  this  method  looks  at  phase  intersection  points 
that  are  either  within  or  outside  the  micrograph  boundary,  hence,  it  is  a  more  accurate  method.  In  general, 
area  particle  counting  can  be  used  for  structures  that  have  isolated  features  of  interest  within  a  matrix 
material,  while  area  tangent  count  is  used  for  more  complicated  microstructures  where  neither  phase  is  fully 
isolated,  similar  to  those  produced  by  the  SHS  method. 

The  phase  size  for  each  microstructure  is  also  expected  to  influence  the  mean  integral  curvature 
values  since  curvature  is  essentially  the  inverse  of  the  radius  of  particles  within  the  structure.  This  can  be 
illustrated  by  looking  at  the  average  phase  size  for  the  TiB2  phase.  Essentially,  if  the  phase  size  is  thought  of 
as  a  particle  size,  then  the  average  curvature  can  be  calculated  by  taking  the  inverse  of  the  particles  radius. 
From  this  notion  and  using  the  TiB2  phase  size  values  plotted  in  figure  7,  the  most  connective  phase  should 
be  sample  C,  which  is  similar  to  that  revealed  by  the  mean  integral  curvature  values  indicated  using  the  area 
tangent  count  method.  Thus,  the  phase  size  is  influencing  the  curvature,  which  is  not  unexpected.  In  terms 
of  the  effects  of  processing  variables,  ball-milling  of  the  powders  prior  to  hot-pressing  was  observed  to  have 
the  greatest  influence  on  phase  size  and  the  microstructural  bias.  Hot-pressing  temperature  and  pressure 
did  not  appear  to  affect  the  microstructural  bias  or  phase  size  significantly.  However,  the  increased  hot- 
pressing  soak  time  did  increase  the  TiB2  phase  connectivity  while  only  slightly  changing  the  grain  (phase) 
size. 


5.  Conclusion 

Microstructurally  biased  two-phase  TiB2 + Al203  ceramics  were  produced  by  hot-pressing  powders 
prepared  using  SHS  and  MM  processing  techniques.  Stereological  measurements  were  used  to 
quantitatively  distinguish  the  connectivity  and  dispersion  of  TiB2  in  the  TiB2+Al203  ceramics. 
Measurements  of  a  number  of  features  and  the  mean  liner  intercept  were  used  to  quantify  the  mean  integral 
curvature  as  a  measure  of  TiB2  phase  connectivity.  The  area  tangent  count  method  was  found  to  be  more 
accurate,  particularly  because  the  complicated  microstructure  of  the  SHS  synthesized  samples  in  which 
neither  phase  is  fully  isolated.  Amongst  the  various  processing  variables,  ball-milling  time  prior  to  hot- 
pressing  of  the  powders  was  observed  to  have  the  greatest  influence  on  phase  size  and  breakdown  of  the 
microstructural  bias. 
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APPENDIX 


A.l  Stereology  of  Mean  Integral  Curvature 

Curvature  has  been  found  to  be  particularly  useful  in  quantifying  phase  connectivity  in  two-phase 
materials.  It  can  be  shown  that  the  connectivity  of  a  two-dimensional  area  is  related  to  the  curvature  of  its 
boundary  [15].  The  curvature,  which  illustrates  how  sharply  a  curve  bends  in  2D  or  3D  space,  can  be 
obtained  for  any  point  along  a  curve  and  has  a  local  value  for  a  unique  point  on  that  curve.  Mathematically, 
curvature  is  defined  as  the  angle,  dOxhax  is  subtended  on  a  finite  arc  segment  divided  by  the  length  of  the 
arc  dA,  as  illustrated  in  figure  Al(a)  for  any  two-phase  material.  The  angle  dO expresses  how  much  the 
curve  has  changed  direction  for  the  segment  dA.  So,  the  average  rate  of  change  is  given  by 


k 


=  d0_ 
dX  ' 


(1) 


The  local  curvature  value  is  comparable  to  the  reciprocal  of  the  radius  of  a  unique  circle  that  exhibits  the 
shape  of  the  curve  at  a  distinctive  point.  For  complex  shapes,  such  as  phases  comprised  of  TiB2  in  an  Al203 
matrix,  a  mean  curvature  can  be  calculated  for  all  points  along  the  curves  formed  from  their  phase 
boundaries.  The  average  curvature  is  then  given  by 


\k{X)dX 

jdX 


The  mean  integral  curvature  can  be  found  by  simply  taking  the  local  mean  curvature  and  integrating  it  over 
the  entire  surface  of  the  curve.  Substituting  equation  (1)  into  (2)  and  evaluating  the  integral,  results  in 


k  = 


A9 
&k  ’ 


(3) 


where,  A9 is  net  change  of  the  angle  subtended  by  all  the  arc  segments  of  the  curve.  If  the  particle  sections 
are  simply  connected,  the  integral  mean  curvature  can  be  estimated  from  the  area  particle  count  method 
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described  by  DeHoff  and  Rhines  [15].  For  a  closed  curve  that  encloses  the  reference  matrix  phase  (TiB2 
surrounding  Al203),  A0=  +2x.  If  the  reference  phase  encloses  the  second  phase  (TiB2  dispersed  in  Al203) 
then,  A0  =  -2n.  The  measure  of  average  curvature  is  then  given  by 


La 


(4) 


In  the  above  equation,  NA+  is  the  number  of  positive  loops  per  unit  area  (or  TiB2  surrounding  Al203  within 
the  TiB2  phase)  and  NA.  is  the  number  of  negative  loops  per  unit  area  (or  TiB2  dispersed  in  Al203  matrix). 
This  can  be  simplified  further  by  recognizing  that  the  total  length  of  the  loops  per  unit  area,  LA  is  actually 
the  perimeter  of  the  closed  loops  formed  by  the  phase  boundary,  which  is  easily  related  to  the  mean  linear 
intercept  length,  I . 


7t 

Y 


(5) 


In  equation  (5),  NL  is  the  average  number  of  intersections  of  phase  boundaries  with  a  test  line  per  unit 

length.  Substituting  for  LA  in  equation  (4)  results  in  the  average  curvature  of  isolated  TiB2  phase  areas  or 
“particles”  and  is  given  by 


k=4NA(J. 


(6) 


For  surfaces  that  enclose  a  second  phase,  the  quantity  is  related  to  the  number  of  positive  and  negative 
loops  divided  by  the  total  perimeter  of  both.  By  convention,  the  net  angle  subtended  by  the  reference 
phase,  in  this  case,  A1203  has  a  positive  value  if  enclosed  by  the  minor  phase  (TiB^  and  has  a  negative  value 
if  the  inverse  is  true.  For  this  reason,  any  time  a  TiB2  “particle”  is  counted,  it  is  considered  a  negative  loop. 
The  same  is  true  for  counting  dispersed  Al203,  which  can  be  thought  of  as  “holes”  in  the  TiB2  phase. 
Determination  of  curvature  then  becomes  nothing  more  than  an  exercise  in  counting  the  numbers  of  TiB2 
particles  and  Al203  holes  and  the  number  of  intersections  that  test  lines  make  with  the  phase 
boundaries  of  the  two  materials.  The  mean  integral  curvature  based  on  loop  count  is,  thus,  determined  by 
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equation  (6).  There  is  no  restriction  on  the  shape  of  connectedness  of  the  “particles”  and  those  of  the 
boundaries  between  the  reference  “particles”  and  the  matrix  area  [16]. 

DeHoff  [15,17]  suggests  another  way  of  calculating  mean  integral  curvature  using  a  technique 
devised  by  Rhines  [15]  for  very  complex  microstructures.  The  definition  of  mean  integral  curvature  is  the 
same  as  in  equation  (1)  but  the  evaluation  of  the  angle  subtended  by  the  curves  40  is  modified.  The 
sampling  procedure  involves  sweeping  a  test  line  across  a  two-dimensional  section  of  the  specimen  and 
counting  where  test  lines  form  tangents  with  the  reference  phase.  This  is  illustrated  in  figure  A  1(b)  for  a 
two-phase  ceramic.  The  total  number  of  tangents  divided  by  the  total  area  evaluated  is  designated  Tj.  The 
total  area  that  is  sampled  by  all  the  test  lines  is  equal  to  NL2,  where  N  is  the  number  of  test  lines  and  L  is  the 
length  of  each  test  line.  The  fraction  of  test  lines  that  form  a  tangent  with  the  phase  boundary  arc  of  length 
dl  is  equal  to  ddh r.  Then  the  number  of  tangents  formed  with  the  arc  length  dl  per  unit  area  of  the  sample 
is  simply  d6/nA.  For  a  curve  of  finite  length,  the  total  number  of  tangents  formed  is 


(7) 


where  0a  is  the  total  angle  subtended  per  unit  area. 

As  described  earlier  for  area  particle  count,  it  is  necessary  to  choose  a  convention  for  defining  the 
reference  phase.  Following  the  same  sign  convention  as  before,  the  Al203  phase  is  again  chosen  as  the 
reference  phase  and  a  positive  sign  denotes  a  convex  arc  segment  while  a  concave  arc  segment  is  denoted 
with  a  negative  sign.  A  convex  arc  segment  has  a  normal  vector  pointing  into  the  reference  phase  (Al203) 
while  the  concave  arc  segments  normal  vector  is  pointing  away  from  the  reference  phase.  What  follows 
then  is  that  for  all  convex  arcs, 


rp  _  @a+ 

1  a+  ~ 


K 


and  for  all  concave  arcs, 


(8) 


(9) 
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Combining  equations  (8)  and  (9)  and  substituting  into  equation  (3)  where,  AA  =  LA  =  ^=>  gives  the 
expression  for  mean  integral  curvature  using  area  tangent  count, 

k  ~  2TA(nJ .  (10) 


A.2  Stereology  of  Phase  Size 

Phase  size  was  measured  for  each  of  the  microstructures  discussed  in  this  work  using  mean  intercept 
length.  The  geometrical  features  were  approximated  and  related  to  an  equivalent  diameter.  To  understand 
how  this  approximation  works,  imagine  an  irregularly  shaped  convex  particle  with  an  array  of  test  lines 
passing  through  it.  The  length  of  each  line  within  the  convex  body  is  measured  and  varies  from  zero  to 
some  maximum  value.  The  mean  value  of  the  intercept  lengths  will  fall  between  these  two  extremes. 
Relating  the  geometry  of  a  sphere  to  mean  intercept  length,  which  is  equal  to  two-thirds  the  diameter,  one 
can  approximate  the  geometrical  features  in  terms  of  an  equivalent  average  diameter. 

The  mean  intercept  length  for  a  single  particle  is  defined  as 


I  N 

i  =—YJii 


(11) 


where  /,•  are  the  intercept  lengths  from  ^intersections  of  test  lines  with  phase  boundaries.  For  a  given  test 
line  length  L, 


(12) 


where  Li  is  the  fractional  length  of  a  test  line  through  a  phase  and  Nl  is  the  number  of  intercepts  per  unit 
length  of  a  test  line.  From  basic  stereological  relations  [15,18]  and  considering  a  multiphase  material  gives 
the  expression 
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(13) 


/ 


=  4 


(Vy)a 

(Sr).' 


where  Vy  is  the  volume  fraction  and  Sy  is  the  total  particle  interface  area  per  unit  volume.  Considering  now 
a  single  spherical  particle  again,  and  substituting  the  volume  of  a  sphere  for  Vy  and  the  area  of  a  sphere  for 
Sy  from  equation  (13),  gives 


J  =  —R  ~—D,  (14) 

3  3 

where  R  is  the  particle  radius  and  D  is  the  particle  diameter.  For  a  two-phase  material,  D  is  actually  the 
phase  size  of  the  a-phase  (reference  phase).  Substituting  equation  (14)  into  equation  (13)  gives 

D=3{Vv)J,  (15) 


which  is  used  to  calculate  the  average  phase  size  of  a  two-phase  material.  This  gives  approximate  phase  sizes 
for  a  complex  two-phase  material  that  is  otherwise  difficult  to  characterize  because  of  irregular  shaped 
features. 
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Figure  1:  Two-phase  TiB2+Al203  ceramics  from  Carney  [6]  showing  interconnected  TiB2  (T@A) 
microstructures  A  and  C  (a,c),  and  dispersed  TiB2  (TinA)  microstructures  B  and  D  (b,d). 

Figure  2:  a)  Mean  integral  curvature  measured  for  microstructures  A,  B,C,  and  D  using  area  particle 
count  and  area  tangent  count  methods,  b)  Comparison  of  mean  integral  curvature 
measurements  using  area  tangent  count  method  for  samples  fabricated  by  Carney  [6]  (prior 
work)  and  Ferranti  [7]  (present  work). 

Figure  3:  Optical  micrographs  of  hot-pressed  sample  cross-sections  fabricated  from  SHS  powders  ball- 
milled  a)  1,  b)  2,  c)  4,  d)  8,  e)  15,  f)  23,  and  g)  30  hours. 

Figure  4:  Effect  of  ball-milling  on  measured  mean  integral  curvature  values  for  SHS  formed  ceramics. 

Figure  5:  Micrographs  of  manually  mixed  samples  with  a)  no  ball-milling  (MM-017)  and  b)  ball-milling 
for  2  hours  (MM-042). 

Figure  6:  Micrographs  of  a)  SHS-039  and  b)  MM-042  samples  both  ball-milled  for  2  hours. 

Figure  7 :  Effect  of  ball-milling  time  on  the  average  measured  phase  size  for  individual  TiB2  and  Al203 
phases  of  SHS  formed  TiB2+Al203  ceramics. 

Figure  8:  Average  measured  phase  sizes  of  TiB2  and  Al203  for  microstructures  A,  B,  C,  and  D made  from 
previous  work  [6]. 

Figure  Al:  Cross-section  of  a  two-phase  material,  (a)  Arc  segments  located  at  the  phase  boundary  subtend 
some  angle  d0,  which  is  used  in  calculating  mean  integral  curvature,  (b)  A  horizontal  test  line 
“sweeping”  across  a  two-dimension  cross-section  forms  positive  and  negative  tangent  points 
illustrating  area  tangent  count  method  for  measuring  mean  integral  curvature. 
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Figure  1:  Two-phase  TiB2+ Al203  ceramics  from  Carney  [6]  showing  interconnected  TiB2  (T@A) 
microstructures  A  and  C  (a,c),  and  dispersed  TiB2  (TinA)  microstructures  B  and  D  (b,d). 
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Figure  2:  a)  Mean  integral  curvature  measured  for  microstructures  A,  B,  C,  and  D  using  area  particle 
count  and  area  tangent  count  methods,  b)  Comparison  of  mean  integral  curvature 
measurements  using  area  tangent  count  method  for  samples  fabricated  by  Carney  [6]  (prior 
work)  and  Ferranti  [7]  (present  work). 


Figure  5:  Micrographs  of  manually  mixed  samples  with  a)  no  ball-milling  (MM-017)  and 
b)  ball-milling  for  2  hours  (MM-042). 
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Figure  6:  Micrographs  of  a)  SHS-039  and  b)  MM-042  samples  both  ball-milled  for  2  hours. 
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Figure  7:  Average  measured  phase  sizes  of  TiB2  and  Al203  for  microstructures  A>  B ,  C,  and 
D  made  from  previous  work  [6]. 
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Figure  Al:  Cross-section  of  a  two-phase  material,  (a)  Arc  segments  located  at  the  phase  boundary 
subtend  some  angle  dd,  which  is  used  in  calculating  mean  integral  curvature,  (b)  A 
horizontal  test  line  “sweeping”  across  a  two-dimension  cross-section  forms  positive  and 
negative  tangent  points  illustrating  area  tangent  count  method  for  measuring  mean 
integral  curvature. 
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Abstract.  A  micromechanical  model  that  provides  explicit  accounts  for  arbitrary  microstructures  and  arbitrary 
fracture  patterns  is  developed  and  used.  The  approach  uses  both  a  constitutive  law  for  the  bulk  solid  constituents 
and  a  constitutive  law  for  fracture  surfaces.  The  model  is  based  on  a  cohesive  surface  formulation  of  Xu  and 
Needleman  and  represents  a  phenomenological  characterization  for  atomic  forces  on  potential  crack/microcrack 
surfaces.  This  framework  of  analysis  does  not  require  the  use  of  continuum  fracture  criteria  which  assume,  for  ex¬ 
ample.  the  existence  of  AT -fields.  Numerical  analyses  carried  out  concern  failure  in  the  forms  of  crack  propagation 
and  microcrack  formation.  Actual  microstructures  of  brittle  alumina/titanium  diboride  (Al2C>3/TiB2)  composites 
are  used.  The  results  demonstrate  the  effects  of  microstructure  and  material  inhomogeneities  on  the  selection  of 
failure  modes  in  this  material  system.  For  example,  the  strength  of  interfaces  between  the  phases  is  found  to 
significantly  influence  the  failure  characteristics.  When  weak  interfacial  strength  exists,  interfacial  debonding  and 
microcrack  initiation  and  growth  are  the  principal  mode  of  failure.  When  strong  interfacial  strength  is  derived  from 
material  processing,  advancement  of  a  dominant  crack  and  crack  branching  are  observed. 

Key  words:  Micromechanical  modeling,  cohesive  force,  fracture,  heterogeneous  materials,  elasticity,  fracture 
modes,  failure  modes,  numerical  simulation,  crack  propagation,  ceramic  composites. 


1*  Introduction 

Microscopically  inhomogeneous  materials  derive  significantly  higher  strength  and  toughness 
from  microscopic  reinforcements  such  as  fibers  and  particles.  Composite  materials  with  fail¬ 
ure  resistance  superior  to  those  of  their  individual  constituents  have  been  developed.  The 
enhancement  is  due  to  the  higher  toughness  and  strength  of  the  additional  phases  and  de¬ 
formation  mechanisms  that  do  not  come  into  play  for  monolithic  materials.  In  addition  to 
different  length  scales  associated  with  material  inhomogeneities,  different  time  scales  are 
also  introduced  by  composite  microstructures  under  transient  loading  due  to  different  stress 
wave  speeds  or  the  inertia  effect.  These  spatial  and  temporal  effects  provide  opportunities 
for  material  property  enhancement.  Material  heterogeneities  also  give  rise  to  multiple  failure 
mechanisms.  For  example,  fracture  can  occur  in  different  phases  and  along  phase  boundaries 
in  heterogeneous  materials. 

Ceramic  composites  with  microstructural  reinforcements  over  a  range  of  size  scales  are 
emerging  materials  with  increasing  applications  in  heat  engines,  gas  turbine  blade  coatings, 
cutting  tools,  drill  bits,  wear  parts,  sensors,  magnetic  recording  media,  structural  components, 
electronic  components,  and  biomedical  devices  (e.g.  prosthetic  articulate  joints,  orthopedic 
load-bearing  hip  implants,  spinal  surgery  implants,  dental  crowns  and  bridges).  These  appli- 
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cations  take  advantage  of  outstanding  mechanical  and  thcrmomechanical  properties  at  high 
temperatures  (Lange,  1979;  Zeng  et  al.,  1992:  Matsui  et  ah,  1994).  good  wear  resistance, 
high  elastic  moduli  and  excellent  chemical  stability  (Jones,  1996)  of  the  materials.  For  in¬ 
stance,  ceramic  nanocomposites  make  high  performance  cutting  tools  for  high-speed  ma¬ 
chining  (Komanduri.  1989).  Ceramic  composites  are  also  excellent  biomaterials  since  they 
are  biocompatible  (fully  oxidized  and  chemically  stable)  and  do  not  cause  adverse  effects 
within  a  physiological  environment.  Conversely,  their  reliability  and  performance  are  unaf¬ 
fected  by  the  biological  environment  in  a  body.  An  important  issue  for  ceramic  materials 
in  these  applications,  industrial  or  biomedical,  is  their  failure  resistance,  including  strength 
and  fracture  toughness  (Komanduri.  1989;  Messer,  1995).  The  fracture  toughness  of  these 
materials  is  at  least  an  order  of  magnitude  lower  than  those  of  metals  and  polymers.  Low  fail¬ 
ure  resistance  is  the  most  important  impediment  to  their  applications,  e.g.  Willmann  (1996a, 
b).  Progress  has  been  made  in  developing  advanced  ceramic  materials  using  the  fact  that  the 
materials  derive  significantly  higher  toughness  from  microscopic  or  nanosized  reinforcements. 
(Niihara  et  al.,  1990,  1991;  Lehman,  1989;  Zhao  et  al.,  1993).  For  example.  Niihara  et  al. 
(1993)  reported  that  a  5  percent  population  of  SiC  nanoparticles  increases  the  tensile  strength 
of  Si3N_i  from  350  MPa  to  1  GPa  and  improves  its  fracture  toughness  from  3.25  MPa^m 
to  4.7MPav'm.  Recently,  alumina/titanium  diboride  (ALCtyTC^)  composites  with  a  wide 
range  of  microstructural  morphologies  demonstrate  a  range  of  failure  resistance  and  a  strong 
dependence  of  fracture  toughness  on  microstructure  in  experiments.  (Logan,  1996).  These 
materials  are  composed  of  titanium  diboride  (TiB2)  reinforcements  embedded  in  a  matrix 
of  alumina  (AL03).  The  principal  mode  of  failure  observed  in  experiments  is  transgranular 
fracture  of  AL03  when  strong  interfacial  strength  between  the  phases  is  obtained  through 
processing.  In  contrast,  the  principal  mode  of  failure  is  interfacial  fracture  through  microc¬ 
rack  formation  when  a  weak  inteifacial  bonding  between  the  constituents  is  found.  Clearly, 
interfacial  strength  significantly  influences  the  overall  behavior  of  the  composites.  Although 
microstructure-induced,  size-dependent  toughening  mechanisms  at  the  micro  and  nano  levels 
are  demonstrated  approaches  for  property  enhancement,  the  physics  for  such  effects  has  not 
been  well  quantified.  In  order  to  develop  more  advanced  materials,  it  is  necessary'  to  character¬ 
ize  the  influences  of  phase  morphology,  phase  length  scale,  and  interfacial  behavior  on  failure 
behavior  and  fracture  toughness  of  these  materials. 

The  dynamic  failure  of  brittle  materials  has  been  extensively  analyzed  by,  e.g.  Shockey  et 
al.  (1974),  Grady  and  Kipp  (1979),  Lankford  (1989),  Shockey  et  al.  (1985),  Brockenbrough 
et  al.  (1988).  Longy  and  Cagnoux  (1989),  Kishi  et  al.  (1990),  Curtin  (1991),  Shockey  et  al. 
(1990),  Suresh  et  al.  (1990),  Yang  et  al.  (1990),  Evans  (1991),  Kishi  (1991),  Kobayashi  (1991), 
Espinosa  et  al.  (1992),  Ahrens  and  Rubin  (1983),  Vekinis  et  al.  (1993).  Lankford  (1994), 
Woodward  et  al.  (1 994),  and  Zhou  and  Curtin  ( 1 995).  Most  available  models  for  the  failure  are, 
for  the  most  part,  continuum  damage  theories  in  which  the  net  effect  of  fracture  is  idealized 
as  a  degradation  of  the  elasticity  modulus,  see  e.g.  Seaman  et  al.  (1985),  Curran  et  al.  (1987, 
1993),  Rajendran  (1994),  Johnson  and  Hohnquist  (1992),  Walter  (1992),  Espinosa  et  al.  (1992, 
1995),  Ravichandran  et  al.  (1995),  and  Gao  et  al.  (1997).  While  capturing  the  macroscopic 
or  effective  response,  these  models  do  not  explicitly  consider  the  discrete  nature  of  fracture 
through  crack  growth  and  coalescence.  Thus,  the  models  lack  the  ability  to  account  for  the 
interaction  between  cracks  and  resolve  specific  failure  modes  and  failure  patterns.  In  addition, 
the  effects  of  microstructural  inhomogeneities  of  different  size  scales  such  as  inclusions,  fibers 
and  grains  on  crack  path  and  fracture  toughness  cannot  be  explicitly  analyzed. 
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The  complex  morphologies  of  material  microstructures  preclude  the  application  of  analyti¬ 
cal  methods.  Explicit  micromechanical  modeling  and  simulation  represent  a  unique  and  attrac¬ 
tive  means  for  analyzing  micro  and  meso  failure  mechanisms  and  for  elucidating  scaling  laws. 
Through  the  consideration  of  representative  samples  of  actual  microstructures,  the  effects 
of  various  fracture  mechanisms  can  be  delineated.  The  required  features  of  this  framework 
should  include 

(1)  explicit  account  of  real,  arbitrary  material  microstructures, 

(2)  explicit  resolution  of  fracture  in  a  non-constrained  (arbitrary  crack  paths  or  microcrack 
patterns)  manner,  and 

(3)  freedom  from  limitations  of  fracture  criteria  applicable  only  over  certain  length  scales 
(e.g.  continuum  criteria  which  assume  the  existence  of  X-fields). 

Combined  use  of  the  cohesive  surface  approach  of  Xu  and  Needleman  (1994,  1997)  and  bulk 
constitutive  laws  is  a  good  candidate  for  providing  such  a  framework.  A  similar  technique 
has  also  been  used  by  Camacho  and  Ortiz  (1996)  and  Ortiz  ( 1996)  in  the  analysis  of  dynamic 
failure  of  materials. 

A  cohesive  finite  element  method  (CFEM)  for  explicit  micromechanical  fracture  analysis 
is  developed  and  used  here.  This  approach  involves  the  combined  use  of  a  cohesive  surface 
characterization  for  crack  surfaces  and  bulk  constitutive  laws  for  solid  constituents.  In  addition 
to  traditional  finite  elements,  all  boundaries  between  the  finite  elements  are  cohesive  surfaces 
serving  as  potential  crack  paths.  Like  in  Xu  and  Needleman  (1994),  the  crack  surfaces  are 
regarded  as  cohesive  surfaces  exhibiting  traction  forces  which  are  functions  of  interfacial 
separations.  The  concept  of  cohesive  crack  faces  can  be  traced  back  to  the  pioneering  work  of 
Dugdale  (1960)  and  Barrenblatt  (1962).  The  additional  interfacial  discretization  in  the  CFEM 
allows  the  cohesive  surfaces  to  permeate  the  whole  microstructure  as  an  intrinsic  part  of 
the  material  characterization.  Consequently,  fracture  is  an  inherent  attribute  of  the  discrete 
model.  Explicit  account  is  taken  of  arbitrary  crack  and  microcrack  patterns  as  w'ell  as  of 
arbitrary  microstructures.  Intergranular  and  transgranular  fractures  evolve  as  natural  outcomes 
of  the  cohesive  responses  within  each  constituent,  the  cohesion  between  the  phases,  and  the 
bulk  constitutive  behavior  of  the  phases.  The  cohesive  relations  represent  phenomenological 
characterizations  of  atomistic  attraction  and  repulsion  forces  which  vary  w'ith  inter-atomic 
displacements.  It  is  surmised  that  characterizations  of  the  inter-atomic  interactions  are  possible 
on  the  nano,  micro,  meso  and  continuum  levels.  When  combined  with  bulk  constitutive  laws 
appropriate  at  the  corresponding  length  scales,  these  cohesive  models  can  provide  a  unified 
framework  for  explicit  account  and  simulation  of  fracture  processes  over  multiple  length 
scales.  This  novel  approach  does  not  require  any  a  priori  crack  initiation  or  propagation  cri¬ 
teria  required  by  traditional  fracture  mechanics  approaches,  such  as  the  attainment  of  critical 
values  of  stress  intensity  factors  or  maintenance  of  constant  energy  release  rate.  The  difference 
is  significant  because  the  concept  or  validity  of  singular  crack  tip  fields  ceases  to  exist  when 
one  approaches  the  grain,  subgrain,  or  atomic  scales.  Since  the  CFEM  model  has  the  inherent 
properties  of  deformation  and  fracture  upon  loading,  it  possesses  a  predictive  power  similar  to 
that  of  discrete  molecular  dynamics  models  concerning  fracture  initiation,  fracture  path,  crack 
speed  and  microcrack  patterning. 

Analyses  in  this  paper  focus  on  the  TiEL/ALCLs  material  system  developed  by  Logan 
(1996)  because  of  its  application  potential  and  because  of  the  fact  that  processing  has  been 
demonstrated  to  be  an  effective  means  to  alter  the  mechanical  properties  of  the  materials 
through  microstructural  modification.  This  selection  is  also  motivated  by  the  opportunity  for 
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a  joint  program  combining  micromcchanica)  modeling,  mechanical  testing  and  materials  syn¬ 
thesis  aimed  at  improving  the  fracture  resistance  of  a  class  of  ceramic  composites.  While 
a  specific  material  system  is  considered  here,  the  approach  for  explicit  fracture  modeling 
can  be  directly  applied  to  other  material  systems.  Understanding  of  microstructure-induced 
toughening  mechanisms  is  also  directly  relevant  for  other  composite  materials. 

2.  Problem  formulation 


To  account  lor  finite  strain  involved  in  crack  tip  regions,  a  Lagrangian  finite  deformation 
formulation  is  used.  The  independent  variables  are  the  position  of  a  material  point  in  the 
reference  configuration  x,  and  time  t.  Relative  to  a  fixed  Cartesian  frame  {£'),  a  material  point 
initially  at  x  occupies  position  x  in  the  current  configuration.  The  displacement  vector  and  the 
deformation  gradient  are  defined  as  u  =  x  -  x  and  F  =  3x/3x.  respectively.  The  principle  of 
virtual  work  includes  a  contribution  from  the  cohesive  surfaces  and  is  written  as 

[s:8TdV-[  T  -5A  dS  =  f  T -SudS  -  [  •  SudV,  (1) 

Jv  Jsim  Js*\i  Jv  d< 

where  s  :  5F  =  s'J8Fji.  s  is  the  nonsymmetric  first  Piola-Kirchhoff  stress;  A  is  the  displace¬ 
ment  jump  across  a  pair  of  cohesive  surfaces;  V,  S„,  and  5ln,  are  the  volume,  external  surface 
area  and  internal  cohesive  surface  area,  respectively,  of  the  body  in  the  reference  configuration. 
The  density  of  the  material  in  the  reference  configuration  is  p.  Also,  5F,  8  A,  and  5u  denote 
admissible  variations  in  F.  A  and  u  respectively.  The  traction  vector  T  and  the  surface  normal 
in  the  reference  configuration  n  are  related  through  T  =  n  s.  The  volumetric  constitutive  law 
is  hyperelastic  so  that 


where  S  =  s  -  F  7  is  the  second  Piola-Kirchhoff  stress.  The  strain  energy  density  IV  is  taken 
to  be 

W  =  iE  :  L  :  E  n\ 


being  the  tensor  of  isotropic  elastic  moduli.  E  and  v  are  the  Young’s  modulus  and  Poisson’s 
ratio,  respectively.  E  is  the  Lagrangian  strain  given  by 

E  =  |(Fr  •  F  -  I).  (5) 

Also  in  the  above  formulas,  II  is  the  fourth  order  identity  tensor,  I  is  the  second  order  identity 
tensor,  I  ®  I  denotes  the  tensor  product  of  two  second  order  tensors,  and  ( )T  and  ( )~T  denote 
inverse  and  inverse  transpose,  respectively.  The  resulting  stress-strain  relation  specifies  a  linear 
relation  between  S  and  E.  i.e., 

S  =  — — E  + - — - (trE)I. 

\  +  v  ( I  +  v)(l  -  2v) K  ’ 


(6) 
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Under  the  conditions  of  infinitesimal  strains,  this  relation  approximates  the  isotropic,  linear 
elastic  behavior  A  discussion  of  linear  stress-strain  relations  for  finite  deformations  is  given 
by  Batra  (1999)* 

The  constitutive  law  for  cohesive  surfaces  relates  the  traction  and  displacement  jumps 
across  crack  surfaces  and  is  also  taken  to  be  hyperelastic  so  that  any  dissipation  associated 
with  separation  is  neglected.  Assuming  the  surface  potential  energy  is  < p ,  the  traction  on  the 
cohesive  surfaces  can  be  derived  through 


(7) 


In  two  dimensions,  the  specific  form  of  0  is  given  by  Xu  and  Needleman  (1994)  as 


0(A)  =  0()  —  0Q 


(8) 


where  A„  =  n  •  A  and  A,  =  t  *  A  are  the  normal  and  tangential  displacement  jumps,  with  n 
and  t  denoting  unit  vectors  along  the  surface  normal  and  tangent  in  the  reference  configuration, 
respectively.  Potential  0  is  written  such  that  <p{  A)|iAm)  =  o  and  therefore  4>0  =  0(A)[|A,_^oc  is  the 
work  of  separation.  Two  special  paths  of  decohesion  are  considered  to  illustrate  the  variations 
of  the  cohesive  traction  components.  The  first  path  is  pure  normal  separation  with  tangential 
traction  component  Tt  —  t  •  T  =  0  and  tangential  separation  A,  —  0.  The  second  path  is  pure 
tangential  separation  with  normal  traction  component  Tn  =  t  T  =  0  and  normal  separation 
A„  —  0.  The  works  of  normal  and  tangential  separations  along  these  two  paths  are  equal  to 
each  other  and  are  taken  to  be  0O  =  eoimxSn  =  y|r maxS,,  with  e  =  exp(l)  —  2.718281829.  8n 
and  <5,/a/2  are  normal  and  shear  displacement  jump  values  at  which  maximum  normal  stress 
amnx  and  maximum  shear  stress  rmax  occur,  respectively. 

The  cohesive  surface  traction  components  obtained  from  (7)  and  (8)  are 


These  relations  are  illustrated  in  Figure  1.  Clearly,  there  are  cross  dependencies  of  normal 
traction  (Tn)  on  shear  displacement  jump  (A,)  and  shear  traction  (Tr)  on  normal  displacement 
jump  (A/7).  In  general,  the  traction  components  first  increase  with  the  displacement  jumps  and 
decrease  with  further  increase  in  separations  after  displacement  jumps  exceed  certain  values. 
Both  traction  components  approach  zero  as  |  A  |  — >  oo.  Figure  1(a)  shows  the  dependence  of 
normal  traction  T„  on  normal  separation  An  and  shear  separation  A,.  When  normal  separation 
is  attempted  while  A,  is  held  to  be  zero  (A„  >  0  and  A,  =  0),  Tn  first  increases  with  A„ 
for  0  <  A/?  <  (5;j;  it  reaches  maximum  stress  crmax  at  A„  =  8n ,  decreases  with  An  thereafter 
and  approaches  zero  as  A„  goes  to  oo.  A  nonzero  shear  separation  A,  always  reduces  the 
value  of  T„  required  for  causing  the  same  amount  of  normal  separation  A„.  The  normal 
traction  component  increases  rapidly  with  negative  normal  separations  (A*  <  0,  representing 
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(a)  (b) 

Figure  1.  An  illustration  of  (he  constitutive  behavior  of  cohesive  surfaces,  (a)  normal  traction,  (b)  shear  traction. 


interpenetration  of  the  crack  surfaces).  This  part  of  the  relation  simulates  the  strong  atomistic 
repulsion  during  compression.  Figure  1(b)  shows  the  variation  of  T,  with  A,  and  A„.  When 
shear  separation  is  attempted  while  A„  is  held  to  be  zero  (A,  ^  0  and  A„  =  0).  T,  first 
increases  with  |A,|  forO  <  |A,|  <  it  reaches  maximum  stress  rinax  at  |A,|  =  5, /a/2. 

decreases  with  |A,|  thereafter  and  approaches  zero  as  |A,|  goes  to  oo.  A  nonzero  normal 
separation  A„  always  reduces  the  value  of  T,  required  for  causing  the  same  amount  of  shear 
separation  A,. 

The  balance  of  energy  requires  that  the  total  external  work  X  done  to  a  body  to  be  equal  to 
the  sum  of  the  strain  energy  W  stored  in  the  bulk  material,  the  kinetic  energy  X  in  the  body 
and  the  cohesive  surface  energy  <t>  in  crack  surfaces,  i.e. 

P  =  [  f  T~dSdr  =  X  +  W  +  d>.  nil 

Jo  Jscx,  dr  v  ' 

where 


„  r  ,  9u  3u 

-W  =  fy  w  dV. 

">  =  /5,'W 


(12) 


Cohesive  energy  is  the  amount  of  energy  spent  on  generating  new  crack  surfaces,  there¬ 
fore  represents  a  measure  of  the  energy  consumption  on  fracture.  Under  the  same  loading 
conditions,  the  higher  the  <J>  required  for  generating  the  same  amount  of  crack  surfaces,  the 
higher  the  fracture  resistance  of  the  material.  In  the  analyses  carried  out  there,  this  energy  is 
used  to  compare  the  fracture  resistance  of  specimens  with  different  microstructures  under  the 
same  overall  remote  loading. 
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Figure  2.  Problem  analyzed:  a  centered-crncked  specimen  under  tensile  loading. 


3.  Problem  analyzed 

Computations  are  earned  out  for  a  center-cracked  specimen,  as  illustrated  in  Figure  2.  The 
specimen  has  an  initial  height  of  2 H  —  1.6  mm  and  an  initial  width  of  2 L  —  1.6  mm.  An 
initial  crack  of  length  2 a{  =  0.4  mm  exists  along  the  £!  axis.  Only  one  half  of  the  specimen 
(£!  >0)  is  discretized  and  modeled  in  the  simulations  because  of  the  symmetry  with  respect 
to  the  £2  axis.  Conditions  of  plain  strain  are  assumed  to  prevail.  The  small  region  in  front  of  the 
crack  tip  contains  microstructures  digitized  from  micrographs  of  actual  composite  materials, 
see  Figure  3.  Inside  this  region,  the  material  inhomogeneities  and  material  distribution  are 
explicitly  modeled.  Digitized  microstructures  of  actual  AFO^/TiB?  composites  are  used  in  the 
analyses.  The  particular  microstructural  morphologies  analyzed  are  shown  in  Figure  4.  The 
microstructure  in  Figure  4(a)  has  a  TiB2  grain  volume  fraction  of  approximately  32  percent. 
The  average  grain  size  is  10  //m.  This  microstructure  is  referred  to  as  the  base  microstructure 
throughout  this  paper.  In  addition  to  the  base  microstructure,  two  variations  are  also  con¬ 
sidered.  The  first  variation,  shown  in  Figure  4(b),  is  obtained  from  the  base  microstructure 
by  reducing  the  volume  fraction  of  the  TiB2  grains  from  32  to  15  percent.  This  reduction  is 
achieved  by  successively  removing  layers  of  TiB2  from  the  grains  and  converting  them  into  the 
matrix  material,  A12C>3.  The  second  variation  is  shown  in  Figure  4(c)  and  is  obtained  from  the 
base  microstructure  by  shrinking  uniformly  the  constituents  in  the  microstructure.  The  ratio  of 
reduction  is  50  percent  in  each  direction  therefore  the  average  grain  size  in  Figure  4(c)  is  one 
half  of  that  in  Figure  4(a).  The  volume  fractions  of  the  phases,  nevertheless,  remain  the  same 
as  those  in  the  base  microstructure.  Each  of  the  phases  behaves  hyperelastically  according 
to  (2)  and  (3).  Depending  on  material  properties,  boundary  conditions  and  loading,  fracture 
is  possible  inside  each  of  the  phases  and  along  the  interfaces  between  the  phases.  The  bulk 
properties  of  each  finite  element  are  those  for  either  the  grains  or  those  for  the  matrix.  The 
properties  of  each  segment  of  potential  fracture  surface  are  specified  according  to  its  location 
as  those  belonging  to  the  A1203  matrix,  the  TiB2  grains  or  the  matrix/grain  interfaces. 

Materials  outside  the  crack-tip  region  are  assumed  to  be  homogeneous  and  are  assigned 
a  set  of  effective  parameters  representative  of  those  for  the  Al203/TiB2  ceramic  composite. 
Both  regions  are  discretized  in  the  same  manner,  using  both  the  bulk  and  the  cohesive  surface 
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Figure  4 .  Microstructural  morphologies  analyzed:  (a)  base  microstructure,  (b)  lower  volume  fraction,  (c)  smaller 


constitutive  descriptions.  The  specimen  is  stress-free  and  at  rest  initially.  Tensile  loading  is 
applied  by  imposing  symmetric  velocity  boundary  conditions  along  the  upper  and  lower  edges 
of  the  specimen.  For  the  results  discussed  here,  the  imposed  boundary  velocity  is  Vq  —  50  m/s 
for  each  edge  with  a  linear  ramp  from  zero  to  this  maximum  velocity  in  the  first  0.01  /zs  of 
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Table  I.  Bulk  material  properties 


Compound 

Density 

(Kg/m3) 

*1 C 

(MPa^/m) 

E 

(GPa) 

V 

Cl 

(m/s) 

C5 

(m/s) 

CR 

(m/s) 

ai2o3 

3990 

4.0 

340 

0.23 

9939 

5885 

5379 

TiB? 

4520 

7.2 

500 

0.12 

10694 

7027 

6267 

Composite 

4120 

3.6 

415 

0.15 

10300 

6609 

5936 

Table  2.  Constitutive  parameters  for  cohesive  surfaces 


Cohesive  surface  pair  rrmax  Tinax  ^,.5/  00 

(GPa)  (GPa)  (//  m)  (J/nr) 


AI2O3 

34 

TiB2 

50 

Composite 

41.5 

A1203  /TiB2 

34 

(Strong  Interface) 

AI2D3  /TiB2 

1.7 

(Weak  Interface) 

78.2 

0.0005 

46.2 

1 1.5 

0.001 

135.9 

95.5 

0.0003 

33.8 

78.2 

0.001 

92.4 

3.91 

0.001 

4.6 

loading.  All  other  specimen  surfaces  have  traction-frec  boundary  conditions.  Specifically,  the 
loading  conditions  are 


—  V„.  i  <  0.01  (is. 

—  L  <  £  <  E; 


±Vo,  i  >  0.0 1  ns. 

=  0.  —L  <£'</.; 

T](±L^7.t)  =  T2(±L.f-,t)  =  Q.  -H  <f-  <  H. 


03) 

(14) 

(15) 


The  material  and  model  parameters  are  listed  in  Table  1  (bulk  properties)  and  Table  2 
(cohesive  surface  constitutive  properties).  For  comparison  and  analysis  purposes,  the  speeds 
for  the  longitudinal  stress  waves  (q_),  the  shear  stress  waves  (cv)  and  the  Rayleigh  surface 
waves  ( cr )  are  also  listed  in  Table  1 .  The  choice  of  the  cohesive  law  parameters  assumes  that 

=  E/10,  with  E  being  the  Young’s  modulus,  and  <p0  =  --l^'’£>A|C,  with  KiC  being  the 
mode-I  fracture  toughness  of  the  materials  in  question. 


4.  Finite  element  method 

Finite  element  discretization  is  based  on  linear-displacement  triangular  elements  arranged  in  a 
crossed-triangle  quadrilateral  pattern.  Neighboring  elements  are  connected  through  cohesive 
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Figure  5.  Distributions  of  L1(  four  different  times  for  a  specimen  of  pure  ABO;,.  Vq  —  50  m/s. 

surfaces.  Hence,  for  the  uniform  mesh  region  in  front  of  the  crack  tip  in  Figure  3.  the  cohesive 
surfaces  are  initially  oriented  along  four  directions,  horizontal  (0C)?  vertical  (90°).  positive 
and  negative  45  degrees  (±45°).  Since  a  very  fine  mesh  is  used  (the  element  size  is  2jum). 
arbitrary  fracture  paths  or  patterns  can  be  resolved.  When  the  finite  element  discretization 
of  the  displacement  field  is  substituted  into  the  principle  of  virtual  work  (1).  the  discretized 
equations  of  motion  take  the  form 


where  U  is  the  vector  of  nodal  displacements,  M  is  the  nodal  mass  matrix  and  R  is  the  nodal 
force  vector  consisting  of  contributions  from  the  bulk  elements  and  the  cohesive  surfaces. 
A  lumped  mass  matrix  is  used  in  (16)  instead  of  the  consistent  mass  matrix  for  reasons 
of  efficiency  and  accuracy  during  explicit  time-integration.  Krieg  and  Key  (1973).  The  ex¬ 
plicit  time-integration  scheme  based  on  the  Newmark  /^-method  with  fi  —  0  and  y  —  0.5 
(Belytschko  et  ah.  1976)  is  employed  to  integrate  (16). 
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Figtnc  6.  Dishibmions  of  022  M  four  different  times  for  11  specimen  with  microstructore  in  Ficure  4(a)  and  a 
strong  interfacial  bonding  between  the  grains  and  the  matrix.  r0  =  M)  m/s. 


5.  Results 

For  comparison  purposes,  a  calculation  is  first  carried  out  for  a  specimen  made  of  pure  AFO, 
matrix  without  reinforcement.  The  material  for  the  whole  specimen  has  the  bulk  and  cohesive 
sui face  piopenies  for  AFO3  listed  in  Tables  1  and  2.  This  calculation  has  an  applied  boundary 
velocity  of  V'0  =  50  m/s.  The  distributions  of  the  vertical  component  of  the  Cauchy  stress  o22 
at  four  dilferent  times  are  shown  in  Figure  5.  The  location  of  the  tip  of  the  initial  crack  is  at 
0.2  mm  in  the  horizontal  direction.  Crack  propagation  initiates  at  approximately  0.065  /is  after 
the  beginning  of  loading.  The  crack  first  propagates  horizontally  briefly  and  then  divides  into 
two  propagating  branches  which  are  oriented  at  approximately  ±45°  from  the  initial  direction 
of  propagation.  This  behavior  has  been  described  by  Xu  and  Needleman  (1994).  The  average 
crack  speed  in  this  case  is  3200  m/s  for  the  period  shown.  Note  that  the  Rayleieh  wave  speed 
in  this  material  is  5379  m/s. 

The  tesults  of  two  calculations  using  the  base  microstructure  in  Figure  4(a)  are  shown 
in  Figures  6  and  7.  These  two  sets  of  results  are  obtained  for  the  same  loading  condition  of 
Vb  =  50  m/s.  the  same  that  as  in  Figure  5.  In  Figure  6.  (he  interfaces  between  the  grains  and  the 
matrix  are  assumed  to  have  a  strong  bonding,  as  characterized  by  a  maximum  normal  stress  of 
=  34  GPa  and  n  maximum  shear  stress  of  rlim  =  78.2  GPa"  These  values  correspond  to  a 
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Figure  7.  Distributions  of  022  at  four  different  times  for  a  specimen  with  microstructure  in  Figure  4(a)  and  a  weak 
interfacial  bonding  between  the  grains  and  the  matrix.  V'q  =  50m/.s. 


work  of  separation  of  <p0  =  92.4  J/m2.  In  Figure  7.  the  interfaces  between  the  two  constituent 
phases  are  assumed  to  have  a  weak  bonding  with  anm  =  1.7GPa,  rirmx  —  3.91  GPa  and 
<f>o  ~  4.6J/nr. 

Figure  6  shows  that  when  strong  interfacial  strength  is  assumed  fracture  occurs  through  the 
extension  of  the  pre-existing  crack.  The  crack  propagation  exhibits  a  zigzag  pattern  because  of 
the  materia]  inhomogeneities.  The  crack  eventually  splits  into  two  branches,  in  a  way  similar 
to  what  is  seen  in  Figure  5  for  fracture  in  uniform  AFO3.  However,  the  crack  extension  before 
branching  is  longer  than  that  in  Figure  5.  The  crack  propagation  occurs  mainly  in  the  matrix 
and  along  grain-matrix  interfaces,  with  occasional  fracture  of  grains.  In  Figure  7,  the  mode 
of  failure  is  initiation,  growth  and  coalescence  of  microcracks.  Clearly,  the  relatively  weak 
interfacial  bonding  allows  microcracks  to  initiate  at  multiple  interfacial  sites  upon  tensile 
loading.  This  process  can  occur  independently  of  the  main  pre-crack.  Coalesced  microcracks 
eventually  link  up  with  the  main  crack,  resulting  in  the  failure  of  the  material. 

The  evolutions  of  the  crack  lengths  for  the  three  cases  discussed  so  far  are  compared  in 
Figure  8.  The  curves  represent  the  histories  of  the  cumulative  length  of  cracks/microcracks  in 
each  constituent  or  along  the  phase  boundaries,  regardless  of  the  number  of  segments  involved. 
The  history  of  combined  total  crack  length  in  all  fracture  sites  is  also  shown.  The  cumulative 
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Figure  8.  A  comparison  of  crack  length  histories  for  the  cases  in  Figures  5-7. 

crack  length  in  the  grains  is  quite  small  and  is,  therefore,  not  plotted.  The  specimen  made 
of  pure  Ah03  shows  the  smallest  crack  length  at  all  stages  of  deformation.  Clearly,  material 
inhomogeneities  allowed  significantly  more  crack  surfaces  to  be  generated  under  the  same 
remote  deformation.  The  total  amount  of  crack  surface  area  generated  in  Figure  6  with  the 
strong  interfacial  bonding  is  approximately  three  times  (~  300 /vm  at  0.08  /vs)  that  in  Figure  5 
("'  'DO /nn  at  0.08 /vs).  The  total  crack  surface  area  for  the  weak  interface  in  Figure  7  is 
24  times  (~  2400 /vm  at  0.08 /vs)  that  in  Figure  5.  There  is  a  difference  in  the  distributions 
of  ft  act ut e  sites  in  the  two  cases  as  well.  Jn  Figure  6.  the  strong  interfaces  cause  cracks  to 
propagate  predominantly  in  the  matrix.  While  in  Figure  7.  the  weak  interphase  bonding  al¬ 
lows  microcracks  to  form  mainly  through  grain-matrix  interfacial  debonding.  The  results  here 
are  consistent  with  experimental  observations  that  materials  with  strong  interphase  bonding 
show  transgranular  fracture  and  materials  with  weak  interphase  bonding  show  intergranular 
fracture.  Logan  (1996). 

The  work  of  separation  for  the  weak  interface  is  only  1/20  that  for  the  strong  interface 
(sec  Table  2),  corresponding  to  significantly  less  energy  consumption  per  unit  crack  surface 
area  generated.  However,  the  higher  amount  of  crack  surface  seen  in  Figure  8  more  than 
compensates  foi  this  lower  energy  value  per  unit  area.  Figure  9  shows  the  evolutions  of 
the  strain  energy  stored  in  the  specimen  materials,  the  kinetic  energy  in  the  specimen  and 
the  cohesive  energy  or  surface  energy  along  cracks  generated.  It  can  be  seen  that  the  rate 
of  energy  dissipation  through  crack  surface  generation  is  actually'  higher  initially  for  weak 
interfaces.  This  is  due  to  the  simultaneous  occurrence  of  microcracks  in  multiple  locations. 
Consequently,  more  energy  is  delivered  into  the  crack  surfaces  during  this  time  period  than  in 
the  case  of  strong  interfaces.  As  crack  growth  continues,  the  cohesive  energy  (<t>)  for  the  strong 
interface  material  eventually  exceeds  that  for  the  weak  interface  case.  The  mode  of  failure 
through  miciocracks  also  corresponds  to  less  strain  energy'  (3V)  jn  the  specimen  and  a  higher 
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Figure  9.  A  comparison  of  energy  evolutions  for  the  cases  in  Figures  5-7. 


amount  of  kinetic  energy  (JC)  in  the  bulk  material.  The  effect  of  TiBj  reinforcements  on  the 
energies  is  clearly  significant  regardless  of  the  interfacial  strength.  Without  the  reinforcements 
or  grains,  all  energies  are  only  small  fractions  of  the  values  for  the  cases  with  reinforcements. 
For  example,  the  cohesive  surface  energy  (<P)  for  the  pure  matrix  material  is  only  4  of  those 
for  the  composites,  see  Figure  9. 

The  distributions  of  cr??  at  four  different  times  for  two  calculations  using  the  microstruc¬ 
tures  in  Figure  4(b,  c)  are  shown  in  Figures  10  and  11.  respectively.  Here,  the  same  strong 
interphase  bonding  strength  as  that  in  Figure^ 6  is  assumed.  It  can  be  seen  that  reducing  the 
grain  volume  fraction  and  reducing  the  grain  size  while  keeping  the  grain  volume  fraction 
constant  have  similar  effects  of  facilitating  the  propagation  of  the  pre-crack  and  expediting 
its  subsequent  branching  into  multiple  cracks.  The  smaller  number  of  grains  and  the  smaller 
size  of  grains  pose  less  impediment  to  the  propagation  of  cracks.  Indeed,  the  results  in  these 
figures  look  closer  to  the  results  in  Figure  5  than  to  the  results  in  Figure  7. 

The  crack  length  histories  for  Figures  6,  10  and  1 1  are  shown  in  Figure  12.  These  curves 
allow  the  effects  of  the  different  microstructures  in  Figure  4  to  be  compared  quantitatively. 
Consistent  with  what  is  seen  in  Figures  10  and  11.  the  lower  grain  volume  fraction  for  mic¬ 
rostructure  B  and  smaller  grain  size  for  microstructure  C  both  cause  the  crack  length  in 
the  matrix  to  increase.  The  result  is  higher  overall  crack  lengths  for  microstructures  B  and 
C.  These  variations  in  microstructnre  do  not  seem  to  show  a  clear  influence  on  the  amount 
of  interfacial  crack  or  microcrack  growth  under  the  conditions  analyzed.  The  corresponding 
energy  evolutions  for  the  three  cases  are  shown  in  Figure  13.  As  a  result  of  the  increases  in 
the  matrix  crack  lengths,  microstructures  B  and  C  have  slightly  higher  cohesive  energy  levels 
than  microstructure  A.  These  results  suggest  that  the  variation  of  cohesive  energy  with  grain 
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Figure  JO.  Distributions  of  o?  2  four  different  times  for  a  specimen  with  microsmicture  in  Figure  4(b)  and  a 

strong  interfacial  bonding  between  the  grains  and  the  matrix.  Ip  —  50  m/s. 


size  oi  grain  volume  fraction  may  not  be  monotonic.  Rather,  phase  distributions  and  phase 
morphologies  may  have  an  impact.  This  effect  has  been  analyzed  in  Zhai  and  Zhou  (1998). 


6.  Discussion  and  conclusions 

A  micromechanical  framework  of  analysis  has  been  developed  and  used  to  provide  explicit 
incorporation  of  arbitrary  material  microstructures  and  to  resolve  arbitrary,  unconstrained 
fracture  patterns  in  heterogeneous,  brittle  solids.  The  approach  combines  descriptions  of  bulk 
constituent  response  and  fracture  surface  cohesion,  as  originally  used  by  Xu  and  Needleman 
(1994).  This  approach  is  especially  appropriate  for  analyzing  microscopic  failure  and  fracture 
over  a  range  of  length  scales  because  material  separation  is  a  natural  outcome  of  constitutive 
behavior,  microstructure  and  loading  in  this  model.  The  formulation  is  free  from  failure  cri¬ 
teria  valid  over  only  certain  length  scales,  e.g.,  continuum  criteria  based  on  the  existence  of 
AT-fields. 

The  dynamic  failure  behavior  in  a  class  of  Al20:,ATiB2  composites  is  simulated  under  the 
context  of  a  centered-cracked  specimen  and  the  conditions  of  plane  strain.  In  this  analysis, 
both  the  bulk  and  interfacial  constitutive  laws  are  hvperelastic.  The  results  demonstrated  the 
effects  of  microstructure  on  dynamic  failure  in  the  materials.  The  bonding  strength  between 
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Figure  //.  Distributions  of  0^2  at  four  different  times  for  a  specimen  with  microstrucuire  in  Figure  4(c)  and  a 
strong  interfacia]  bonding  between  the  grains  and  the  matrix.  Vq  =  50  m/s. 


the  alumina  matrix  and  the  titanium  diboride  reinforcements,  which  is  a  function  of  conditions 
during  hot-pressing  of  the  composites,  is  found  to  significantly  influence  the  failure  modes  in 
the  composites.  When  strong  bonding  exists,  the  mode  of  failure  is  the  extension  of  cracks 
and  crack  branching.  When  weak  bonding  exists,  the  mode  of  failure  is  microcrack  growth 
and  coalescence.  The  framework  of  analysis  also  allows  the  evolution  of  crack  lengths  and 
the  evolution  of  energies  to  be  tracked.  Calculations  have  shown  that  under  the  same  overall 
remote  deformation  microscopic  TiB?  reinforcements  significantly  increase  the  total  cohesive 
energy  consumed  in  the  generation  of  crack  and  microcrack  surfaces. 
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Figure  J2.  A  comparison  of  crnck  length  histories  for  the  three  different  microsmictures  in  Fleure  4. 
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Abstract:  The  dynamic  compressive  strength  and  microscopic  failure  behavior  of  TiB2/Al203  ceramic 
composites  with  a  range  of  microstructural  morphologies  and  size  scales  are  analyzed.  A  split 
Hopkinson  pressure  bar  is  used  to  achieve  loading  rates  of  the  order  of  400  s'1.  The  dynamic 
compressive  strength  of  the  materials  is  found  to  be  between  4.3  and  5.3  GPa,  indicating  a  strong 
dependence  of  strength  on  microstructure.  Microstructures  with  finer  phases  as  measured  by  linear 
intercept  length  (LIL)  have  higher  strength  levels.  The  dynamic  strength  levels  are  approximately  27% 
higher  than  the  values  of  3-4  GPa  measured  at  quasi-static  loading  rates  for  these  materials.  These 
strength  levels  are  also  higher  than  the  strength  levels  of  monolithic  TiB2  and  A1203  under  similar 
dynamic  conditions.  Scanning  electron  microscopy  (SEM)  and  energy  dispersive  spectrometry  indicate 
that  failure  associated  with  the  A1203  phase  is  transgranular  cleavage  in  all  microstructures  and  failure 
associated  with  the  TiB2  phase  is  a  combination  of  transgranular  cleavage  and  intergranular  debonding 
and  varies  with  the  microstructures.  The  measured  compressive  strength  of  the  materials  directly 
correlates  with  the  fraction  of  TiB2-rich  areas  on  fracture  surfaces. 


INTRODUCTION 

One  method  to  reduce  brittleness  and  enhance 
failure  resistance  of  ceramics  is  the  development  of 
ceramic  composites  (Niihara  et  al.,  1991).  Significant 
influence  of  microstructural  effects  on  properties  has 
been  reported.  For  example,  Niihara  et  al.  (1991) 
reported  that  a  5%  population  of  SiC  nanoparticles 
increases  the  tensile  strength  of  Si3N4  from  350  MPa 
to  1  GPa  and  improves  its  fracture  toughness  from 
3.25  MPaVm  to  4.7  MPaVm.  Although 
microstructure-induced,  size-dependent  toughening 
mechanisms  at  the  micro  and  nano  levels  are 
demonstrated  approaches  for  property  enhancement, 
so  far  such  effects  have  not  been  well  quantified.  In 
order  to  develop  more  advanced  materials,  it  is 
necessary  to  characterize  the  influences  of  phase 
morphology,  phase  length  scale,  and  interfacial 
behavior  on  fracture  toughness.  Recently,  two-phase 
ceramics  of  titanium  diboride/alumina  (TiB2/Al203) 
with  a  range  of  phase  sizes  and  phase  morphologies 


have  been  developed  (Logan,  1992a,b).  These 
materials  have  shown  a  wide  range  of  fracture 
toughness  values  and  some  of  the  values  are  higher 
than  those  of  both  constituents  produced  separately  in 
bulk.  These  materials  provide  an  opportunity  to  study 
the  correlation  between  microstructure  and 
mechanical  behavior. 

In  many  applications,  intensive  dynamic  loading 
occurs  under  normal  operating  conditions.  Examples 
include  particle  impact  on  ceramic  turbine  blades 
(Bilek  &  Helesic,  1990),  contact  of  a  high  speed 
cutting  tools  with  workpiece  (Komanduri,  1989),  and 
impact  of  ballistic  projectiles  with  ceramic  armor 
(Anderson  &  Morris,  1991).  The  dynamic  behavior 
of  materials  can  be  dramatically  influenced  by 
microstructural  characteristics,  including  phase  size, 
phase  morphology,  composition  and  texture 
(Viechnicki  et  al,  1991).  Thus,  it  is  crucial  to 
develop  intrinsic  relationships  between  fabrication, 
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resulting  microstructure,  and  dynamic  properties  in 
order  to  tailor  ceramic  systems  to  the  needs  of 
applications  (Bilek  &  Helesic,  1990).  In  order  to 
establish  the  relationships,  varying  processing 
approaches  have  been  used  to  produce  two-phase 
TiB2/Al203  ceramics  with  four  different  biased 
microstructures.  These  microstructures  are  analyzed 
in  this  research. 


Fig.  1  Microstructures  of  Materials  Analyzed  (TiB2 
white,  A1203  dark) 


MATERIALS 

The  materials  are  produced  through  self- 
propagating  high  temperature  synthesis  and 
mechanical  mixing  of  powders  followed  by  hot 
pressing  (Logan  1992b;  Carney  1997).  Each  material 
has  a  nominal  composition  of  70%  A1203  and  30% 
TiB2  by  weight.  Figure  1  shows  the  micro  structures 
of  the  four  Al203/TiB2  composites  analyzed.  The 
light  areas  consist  of  TiB2  and  the  dark  areas  are 
A1203.  The  microstructures  vary  in  phase  distribution 
and  length  scales.  This  microstructure  shows 
intertwining  of  alumna  and  titanium  diboride.  The 
microstructure  is  unique  in  that  its  complicated 
morphology  precludes  the  identification  of  simple 
crack  paths  through  each  phase  or  along  phase 
interfaces.  The  microstructure  in  Fig.  1(b)  consists  of 
particles  of  TiB2  embedded  in  an  A1203  matrix.  The 
microstructure  in  Fig.  1(c)  consists  of  isolated  A1203 


matrix  islands  surrounded  by  a  network  of  TiB2 
reinforcement.  The  isolated  matrix  areas  are 
approximately  20x100  Jim  in  size.  A  dispersion  of 
small  A1203  particles  nanometers  in  size  are  scattered 
in  the  TiB2  phase.  Clearly,  two  length  scales  are 
operative  in  this  microstructure,  influencing  its 
strength  and  fracture  toughness.  The  microstructure 
in  Fig.  1(d)  is  similar  to  that  in  Fig.  1(b)  except  for 
significantly  larger  TiB2  particle  sizes.  The  average 
TiB2  particle  size  is  5-10  |im  in  Fig.  1(d)  and  at  least 
one  or  two  orders  of  magnitude  smaller  in  Fig.  1(b). 
The  smallest  particles  in  Fig.  l-l(b)  are  of  submicron 
or  nanometer  sizes.  To  quantify  the  different  phase 
morphologies  and  sizes,  the  linear  intercept  length 
(LIL)  for  the  two  phases  are  measured  on  digital 
images  of  the  microstructures.  In  the  analysis  here, 
the  average  LIL  is  used  as  a  measure  of  phase  size. 

EXPERIMENTS 

In  order  to  characterize  the  response  of  the 
materials  to  dynamic  compressive  loading,  a  split 
Hopkinson  pressure  bar  apparatus  is  used.  Soft 
recovery  of  specimens  is  achieved  through  the  use  of 
a  momentum  trapping  technique  developed  by 
Nemat-Nasser  et  al.  (1991).  This  mechanism  allows 
specimens  to  be  recovered  after  loading  under  a 
single  stress  pulse,  eliminating  any  unintended 
reloading.  The  specimen  fragments  were  collected  for 
postmortem  characterization  of  microscopic  failure. 


A  B  C  D 


Microstructure 

Fig.  2  Compressive  Strength  and  Density 

COMPRESSIVE  STRENGTH 
The  measured  strengths  of  the  materials  and  their 
densities  are  summarized  in  Fig.  2.  The  correlation 
between  the  strength  and  the  linear  intercept  length  is 
shown  in  Fig.  3.  Microstructures  A  and  D  show 
significantly  higher  strengths  than  microstructures  B 
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and  C.  Since  material  B  has  a  much  higher  level  of 
porosity  than  the  other  materials,  attention  is  given 
only  to  materials  A,  C,  and  D.  Higher  strengths 
coincide  with  smaller  LIL.  This  can  be 
phenomenologically  explained.  Microstructures  with 
smaller  LIL  values  present  a  more  heterogeneous 
media  for  propagating  cracks,  inducing  more  tortuous 
cracks  paths  and,  therefore,  increasing  energy 
dissipation  for  crack  growth.  The  micromechanical 
quantification  of  the  fracture  resistance  of  the  four 
microstructures  accounts  for  the  actual  phase 
distributions  and  arbitrary  fracture  patterns. 


Fig.  3  Correlation  between  Strength  and  Linear 
Intercept  Length 

FRACTOGRAPHY 

Coinciding  energy  dispersive  spectrometry  maps 
showing  the  distributions  of  Al203-rich  and  TiB2-rich 
areas  on  fracture  surfaces  are  taken  and  shown  in  Fig. 
4.  The  use  of  these  phase  maps  and  fractographs 
greatly  facilitates  the  identification  of  phases  on 
fracture  surfaces  and  the  intergranular  or 
intragranular  nature  of  the  fracture  process.  The 
failure  process  in  Microstructure  A  shows  large 
regions  of  A1203  surrounded  by  a  continuous  network 
of  TiB2.  The  A1203  areas  are  smooth,  indicating 
transgranular  cleavage  inside  A1203  grains.  Failure 
involving  the  TiB2  phase  appears  to  be  primarily 
intergranular  pullout.  The  overall  failure  is  a 
combination  of  transgranular  cleavage  inside  A1203 
and  intergranular  separation  along  Al203/TiB2 
boundaries.  Approximately  74%  of  the  fracture 
surfaces  is  cleavage  planes  within  the  A1203  and  26% 
are  associated  with  the  TiB2.  Note  that  the  volume 
fraction  of  A1203  in  this  material  is  72.6%.  The 
fracture  surfaces  for  microstructure  B  are  dominated 
by  small  regions  of  transgranular  cleavage  of  the 


A1203  phase.  Failure  involving  the  TiB2  phase 
occurred  primarily  by  intergranular  pullout.  TiB2 
particles  are  relatively  evenly  distributed  on  the 
fracture  surfaces,  suggest  a  fracture  process  which  is 
not  significantly  affected  by  the  presence  of  the 
particles.  Also,  numerous  pores  were  observed  in  the 
A1203  phase,  consistent  with  the  measured  low 
density  (95%  of  theoretical  value)  of  this  material. 
Clearly,  excessive  porosity  is  detrimental  to  the 
resistance  to  crack  propagation  in  this  material. 
Approximately  84%  of  the  fracture  surfaces  is  A1203- 
rich,  significantly  higher  than  the  72.6%  volume 
fraction  of  A1203  in  this  material.  Thus,  the 
occurrence  of  fracture  is  predominantly  in  the  A1203 
phase.  Since  A1203  has  a  lower  strength  and  lower 
fracture  energy  compared  with  TiB2,  this  higher  ratio 
of  fracture  surfaces  in  A1203  appears  to  be  the 
primary  reason  for  the  observed  low  strength  reported 
earlier.  A  primarily  transgranular  mode  of  fracture 
inside  A1203  is  also  observed  for  microstructure  C. 
Fracture  surfaces  are  dominated  by  large  areas  of 
cleavage  within  the  A1203  phase.  However,  failure 
involving  the  TiB2  phase  occurred  through  a 
combination  of  transgranular  cleavage  and 
intergranular  pullout.  Most  parts  of  the  fracture 
surfaces  associated  with  the  TiB2  show  transgranular 
cleavage.  The  fracture  surface  characteristics  reflect 
the  morphology  of  this  microstructure-large  areas  of 
cleavage  in  A1203  surrounded  by  a  continuous 
network  of  TiB2.  The  fact  that  cleavage  in  A1203 
occurs  over  multiple  grains  contributes  to  its  lower 
failure  resistance,  due  to  the  fact  that  the 
effectiveness  of  the  TiB2  as  a  reinforcing  phase  is 
decreased.  Approximately  77%  of  the  fracture 
surfaces  is  Al203-rich,  higher  than  the  72.6%  volume 
fraction  of  A1203  in  the  material.  Thus,  there  appears 
to  be  a  slight  preference  for  failure  to  occur  in  the 
A1203  phase.  The  failure  in  microstructure  D  is 
unique  in  that  transgranular  cleavage  is  the  primary 
fracture  mechanism  for  both  A1203  and  TiB2. 
Although  failure  involving  the  TiB2  phase  occurred 
through  a  combination  of  transgranular  cleavage  and 
intergranular  pullout,  the  dominant  mechanism  is 
cleavage.  This  is  in  sharp  contrast  to  what  is  observed 
for  the  other  micro  structures.  Evenly  distributed  TiB2 
particles  across  the  fracture  surfaces  echo  the 
microstructural  morphology  of  fine  particles 
embedded  in  an  A1203  matrix.  It  also  appears  that  the 
homogeneous  distribution  of  TiB2  reinforcement 
inhibits  intergranular  separation,  thus  forcing  cracks 
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to  go  through  the  stronger  TiB2  phase  and  enhancing 
failure  resistance.  The  rough  and  ragged  nature  of  the 
surface  is  in  contrast  to  the  relatively  smooth 
appearance  of  the  surfaces  for  microstructures  A,  B, 
and,  most  prominently,  C.  Approximately  71%  of  the 
fracture  surfaces  is  Al203-rich,  lower  than  the  72.6% 
volume  fraction  of  A1203  in  this  material.  Note  that 
this  is  the  only  microstructure  that  shows  a  TiB2-rich 
fracture  surface  fraction  higher  than  its  corresponding 
volume  fraction  in  the  material.  This  bias  toward  the 
TiB2  points  to  an  unusual  shift  of  failure  into  the 
stronger  phase  and  can  be  associated  with 
significantly  higher  strength. 


Fig.  4  Fracture  Surfaces 


DISCUSSION 

The  analysis  has  yielded  dynamic  compressive 
strength  between  4.3  and  5.3  GPa,  indicating  a  strong 
dependence  of  strength  on  micro  structure. 
Microstructures  with  finer  phases  as  indicated  by 
smaller  linear  intercept  length  values  have  higher 
strength  levels.  The  dynamic  strength  levels  are 


approximately  27%  higher  than  the  values  of  3-4  GPa 
measured  at  quasi-static  loading  rates.  Scanning 
electron  microscopy  and  energy  dispersive 
spectrometry  indicate  that  failure  associated  with  the 
A1203  phase  is  transgranular  cleavage  in  all 
microstructures.  Failure  associated  with  the  TiB2 
phase  is  a  combination  of  transgranular  cleavage  and 
intergranular  debonding  and  varies  with  the 
microstructures.  The  fraction  of  TiB2-rich  fracture 
surface  areas  directly  correlates  with  the  measured 
compressive  strength  of  the  materials.  The  addition 
of  titanium  diboride  reinforcement  clearly  improves 
the  failure  resistance  of  the  materials. 
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ABSTRACT:  Dynamic  crack  propagation  in  ceramic  composites  is  analv?cd  numerically  The 
simulations  concern  the  effects  of  micros! met ura I  morphologies  on  fracture.  The  analysis  con- 
r  Ain  Ph‘,sc.thMri^unons  in  actual  microstructures  of  alumina/titanitim'dihorkle 

1 1  y  c<>mposites.  The  microstructures  analyzed  have  different  phase  morpholoeies  and 
ciitterem  phase  sizes  over  an  order  of  magnitude  in  length  (from  1-2  to  10-20  am).  A  micro- 
mechanical  model  that  provides  explicit  account  for  arbitrary  microstructures  and  arbitrary 
tract  lire  patterns  is  developed  and  used.  The  approach  uses  hoih  a  constitutive  law  for  the  bulk 
so  id  constituents  and  a  constitutive  law  for  fracture  surfaces.  The  model  is  based  on  the 
cohesive  surface  formulation  of  Xu  and  Needleman  and  represents  a  phenomenological  char¬ 
acterization  for  atomic  forces  on  potential  crack /microcrack  surfaces  This  framework  of  anal- 
ysts  dc>es  not  require  the  use  of  any  fracture  criteria.  Instead,  fracture  evolves  as  an  outcome 
of  bulk  material  response,  mterfacinl  behavior,  and  applied  loading.  This  approach  provides  a 
umfied  and  self-consistent  treatment  of  mixed  mode  fracture.  The  evolutions  of  crack  lengths 
in  different  phases  and  along  interphnse  interfaces  are  calculated  to  track  crack  orouth  The 
overall  local  crack  speed,  defined  as  the  time  rare  of  change  of  arc  leneth  a  Ion  if  ziczacoino 
crack  paths,  is  found  to  reach  the  intersonic  range,  i.e..  greater  than  the  shear  wave  snccdTnnd 
smaller  than  the  longitudinal  wave  speeds  in  the  constituent  phases.  The  model  also  allows  the 
energy  release  rate  to  be  evaluated  easily.  For  the  same  amount  of  crack  surfaces  ccnerated. 
the  average  energy  release  rates  for  fracture  panerns  in  four  microstructures  anah  zed  differ  by 
up  to  20 /r.  The  results  demonstrate  that  larger  TiB,  reinforcements  sienificantK  impede  crack 
propagation  and  increase  the  fracture  resistance  of  the  compositions  indicated  bv  hmher 
average  energy  release  rate  values.  c 

KEYWORDS:  micromechanical  modeling,  cohesive  force,  fraemre.  inhomnecneonx  materials, 
elasticity.  mixed-mode  fracture,  fracture  mode,  numerical  simulation,  crack  propn«ation  ce- 
ramie  composites  r  c 


Nomenclature 

Longitudinal  wave  velocity 
Shear  wave  velocity 
Rayleigh  wave  velocity 
Tmax  Maximum  shear  stress  for  a  cohesive  surface  pair 

°niax  Maximum  tensile  stress  for  a  cohesive  surface  pair 

Characteristic  length  for  normal  separation 
^  Characteristic  length  for  tangential  separation 


,  'Ih£  of  N*echnnical  Engineering.  Georgia  Institute  of  Technology.  At- 

t<  ita.  OA  30332-0405.  (M.  Zliou  is  author  to  whom  correspondence  should  be  addressed.) 
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<i>0 

ad>cr3ck/^/ 


Surface  energy  per  unit  area 
Fracture  energy  associated  with  crack  growth 
Energy  release  rate  per  apparent  crack  growth 
Energy  release  rate  per  local  crack  growth 


Microscopically  inhomogeneous  materials  derive  significantly  higher  strength  and  tough¬ 
ness  from  microscopic  reinforcements  such  as  fibers  and  particles.  The  enhancement  is  due 
to  the  higher  toughness  and  strength  of  the  additional  phases  and  deformation  mechanisms 
that  do  not  come  into  play  lor  monolithic  materials.  In  addition  to  different  length  scales 
associated  with  material  inhomogeneities,  different  time  scales  are  also  introduced  by  com¬ 
posite  microstructures  under  transient  loading  due  to  different  stress  wave  speeds  or  the 
inertia  effect.  These  spatial  and  temporal  effects  provide  opportunities  for  material  property 
enhancement  through  multiple  deformation  and  failure  mechanisms.  For  example,  material 
inhomogeneities  and  reinforcements  force  cracks  to  follow  tortuous  paths  therefore  increasing 
dissipation  and  fracture  resistance.  Micromechanical  modeling  and  simulation  represents  a 
unique  and  attractive  means  for  delineating  the  effects  of  micro  and  meso  failure  mecha¬ 
nisms.  The  dynamic  failure  of  brittle  materials  has  been  extensively  analyzed,  for  example, 
[1-20].  Most  available  models  for  the  failure  are.  for  the  most  part,  continuum  damage 
theories  in  which  the  net  effect  of  fracture  is  idealized  as  a  degradation  of  the  elasticity 
modulus,  for  example,  see  Refs  21-29.  While  capturing  the  macroscopic  or  effective  re¬ 
sponse.  these  models  do  not  explicitly  consider  the  discrete  nature  of  fracture  through  crack 
growth  and  coalescence.  Thus,  the  models  lack  the  ability  to  account  for  the  interaction 
between  cracks  and  resolve  specific  failure  modes  and  failure  patterns.  In  addition,  the  effects 
of  microstuctural  inhomogeneities  of  different  size  scales  such  as  inclusions,  fibers,  and 
grains  on  crack  path  and  fracture  toughness  cannot  be  explicitly  analyzed. 

The  complex  morphologies  of  material  microstructures  precludes  the  application  of  ana¬ 
lytical  methods.  Explicit  micromechanical  modeling  and  simulation  represent  a  unique  means 
for  analyzing  micro  and  meso  failure  and  for  elucidating  scaling  laws.  Through  the  consid¬ 
eration  of  representative  samples  of  actual  microstructures,  the  effects  of  various  fracture 
mechanisms  can  be  delineated.  The  required  features  of  this  framework  should:  (1)  allow 
explicit  account  of  real,  arbitrary  material  microstructures.  (2)  permit  explicit  modeling  of 
fracture  in  a  nonconstrained  (arbitrary  crack  paths  or  microcrack  patterns)  manner,  and  (3) 
be  free  of  limitations  ot  fracture  criteria  applicable  only  over  certain  length  scales  (such  as 
continuum  criteria  that  assume  the  existence  of  A'- fields).  Combined  use  of  the  cohesive 
surface  approach  of  Xu  and  Needleman  [30  Jj ]  and  bulk  constitutive  laws  is  a  good  can¬ 
didate  for  providing  such  a  framework.  Explicit  simulation  of  fracture  and  fragmentation 
have  also  been  carried  out  by  Camacho  and  Ortiz  [32]  and  Ortiz  [33]  using  a  cohesive  model 
in  the  analysis  of  dynamic  failure  of  materials. 

Alumina/titanium  diboride  (AECL/TiB.)  composites  with  a  wide  range  of  microstructurai 
morphologies  demonstrate  a  range  of  failure  resistance  and  a  strong  dependence  of  fracture 
toughness  on  microstructure  in  experiments,  see  Logan  [34].  These  materials  are  composed 
of  titanium  diboride  (TiB-,)  reinforcements  embedded  in  a  matrix  of  alumina  (A1203).  The 
principal  modes  of  failure  observed  in  experiments  are  transgranular  and  intergranular  frac¬ 
ture  of  A1203  matrix,  reflecting  the  effects  of  different  phase  distributions,  different  size 
scales,  and  different  interfacial  bonding  strengths  that  result  from  various  processing  con¬ 
ditions.  For  example,  the  principal  mode  of  failure  is  interfacial  fracture  through  microcrack 
formation  when  a  weak  interfacial  bonding  between  the  constituents  is  found.  Clearly,  in¬ 
terfacial  strength  significantly  influences  the  overall  behavior  of  the  composites.  A  numerical 
analysis  on  the  effect  of  interfacial  strength  on  failure  modes  has  been  conducted  in  Zhai 
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and  Zhou  f.?5).  The  analyses  here  concern  the  effects  of  microstructure  on  the  fracture 
resistance  of  AI,0.,/TiB,  composites.  The  focus  is  crack  path,  crack  speed,  and  energy  re¬ 
lease  rate.  Explicit  account  is  taken  of  arbitrary  crack  and  microcrack  fracture  patterns  as 
well  as  of  arbitrary  microstructures.  As  in  Xu  et  al.  [30).  the  crack  surfaces  are  regarded  as 
cohesive  surfaces  exhibiting  traction  forces  that  are  functions  of  interfacial  separations.  These 
relations  represent  phenomenological  characterizations  of  the  atomistic  attraction  and  repul¬ 
sion  forces  that  vary  with  inter  atomic  displacements.  All  finite  element  boundaries  between 
elements  are  cohesive  surfaces  serving  as  potential  crack  paths.  Arbitrary  crack  and  micro- 
crack  patterns  can  be  resolved  through  the  use  of  fine  meshes.  This  approach  does  not  require 
any  crack  initiation  and  propagation  criterion.  Instead,  the  cohesive  surface  constitutive  law 
allows  fracture  to  evolve  as  a  natural  outcome  of  the  combined  effects  of  bulk  constituent 
response,  interfacinl  behavior,  and  applied  loading. 


Problem  Formulation 

To  account  for  the  finite  strains  involved  in  crack  tip  regions,  a  Lagrangian  finite  defor¬ 
mation  formulation  is  used.  The  independent  variables  are  the  position  of  a  material  point 
in  the  reference  configuration  *.  and  time  /.  Relative  to  a  fixed  Cartesian  frame  (£'}.  a 
material  point  initially  at  at  occupies  position*  in  the  current  configuration.  The  displacement 
vector  and  the  deformation  gradient  are  defined  as  u  =  x  -  x  and  F  =  tix/dx.  respectively. 

The  principle  of  virtual  work  includes  a  contribution  from  the  cohesive  surfaces  and  is  written 
as 


j  s.8/*r/l  j  1  •  8 \dS  -  J  T  •  8 udS  -  J  p  •  8 udV  (I) 

v  $.<ii  v.  i 

where  5:  8 F  =  j"8F/f.,  5  is  the  nonsymmctric  first  Piola-Kirchhoff  stress;  A  is  the  displacement 
jump  across  a  pair  of  cohesive  surfaces:  V.  Scxr  and  Sin(  are  the  volume,  external  surface  area, 
and  internal  cohesive  surface  area,  respectively,  of  the  body  in  the  reference  configuration' 
The  density  of  the  material  in  the  reference  configuration  is  f>.  Also.  8 F.  8A  and  8«  denote 
admissible  variations  F.  A  and  11  respectively.  The  traction  vector  T  and  the  surface  normal 
in  the  reference  configuration  n  are  related  through  T  -  n  •  s.  The  volumetric  constitutive 
law  is  hyperelastic  so  that 


S  = 


dW 

c)E 


(2) 


where  S  -  s  •  F  T  is  the  second  Piola-Kirchhoff  stress.  The  strain  energy  density  IV  is  taken 
to  be  v 


W = ~ ELE 


with 


(3) 


L  = 


E 

I  +  v 


rV®') 


(4) 
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being  the  tensor  of  isotropic  elastic  moduli.  E  and  v  are  the  Young’s  modulus  and  Poisson’s 
ratio,  respectively.  E  is  the  Lagrangian  strain  given  by 

E  =  ~{FT  •  F  -  1)  (5) 

Also  in  the  above  formulas,  11  is  the  fourth  order  identity  tensor,  I  is  the  second  order 
identity  tensor,  /®/  denotes  the  tensor  product  of  two  second  order  tensors,  and  (  )r  and 
X  )"T  denote  inverse  and  inverse  transpose,  respectively. 

The  constitutive  law  for  cohesive  surfaces  relates  the  traction  and  displacement  jumps 
across  crack  surfaces  and  is  also  taken  to  be  hyperelastic  so  that  any  dissipation  associated 
with  separation  is  neglected.  Assuming  the  surface  potential  energy  is  <J>,  the  traction  on  the 
cohesive  surfaces  can  be  derived  through 


(6) 


In  two  dimensions,  the  specific  form  of  <j>  used  here  is 


<!>(A)  <f>0  <j>0 


(7) 


where  A„  -  n  *  A  and  A,  =  t  ■  A  are  the  normal  and  tangential  displacement  jumps,  with  n 
and  t  denoting  unit  vectors  along  the  surface  normal  and  tangent  in  the  reference  configu¬ 
ration,  respectively.  Potential  <J>  is  written  such  that  <f>(A)|w„0  =  0  and  therefore  <|>0  = 
4>(A)|uh.  is  the -work  of  separation.  Two  special  paths  of  decohesion  are  considered  to 
illustrate  the  variations  of  the  cohesive  traction  components.  The  first  path  is  pure  normal 
separation  with  tangential  traction  component  Tt  =  t  •  T  =  0  and  tangential  separation  A,  = 
0.  The  second  path  is  pure  tangential  separation  with  normal  traction  component  Tn  =  t  •  T 
—  0  and  normal  separation  A„  -  0.  The  works  of  normal  and  tangential  separations  along 
these  two  paths  are  equal  to  each  other  and  arejaken  to  be  <J>0  =  comax  8„  =  vW2titi.1x  5., 
with  e  —  exp  (I)  =  2.718281829.  and  5f/V2  are  normal  and  shear  displacement  jump 
values  at  which  maximum  normal  stress  crm.,x  and  maximum  shear  stress  Tmax  occur, 
respectively. 

The  cohesive  surface  taction  components  obtained  from  Eqs  6  and  7  are 


7„(A)  = 


(8) 


and 


U  A) 


HA 

8? 


(9) 


These  relations  are  illustrated  in  Fig.  1.  Clearly,  there  are  cross  dependencies  of  normal 
traction  on  shear  displacement  jumps  and  shear  traction  on  normal  displacement  jumps. 
When  displacement  jumps  exceed  certain  values,  both  traction  components  decrease  w'ith 
further  increase  in  separations  and  approach  zero  as  |A|  —  oc.  Under  negative  normal  sepa- 
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rations  (interpenetration),  the  normal  traction  component  increases  rapidly,  simulating  the 
strong  atomistic  repulsion  during  compression. 

The  balance  of  energy  requires  that  the  total  external  work  done  to  a  body  to  be  equal 
to  the  sum  of  the  strain  energy  u*  stored  in  the  bulk  material,  the  kinetic  energy  X  in  the 
body  and  the  cohesive  surface  energy  <I>.  meaning 


!l>  =  | 7 

’  ■  it  dS  = 

*K  +  (r 

+  <I> 

(10) 

5fCX( 

where 

.  *'  =  J 

|  1  du 

2PT' 

du  ,  ^ 
—  dV 

dt 

V=j 

\ 

HWV 

(11) 

and 

<F  = 

j  4>dS 

The  cohesive  energy  is  the  amount  of  energy  spent  on  generating  new  crack  surfaces, 
therefore  represents  a  measure  of  the  energy  consumption  on  fracture.  Under  the  same  load¬ 
ing  conditions,  the  higher  the  <J>  required  for  generating  the  same  amount  of  crack  surfaces, 
the  higher  the  fracture  resistance  of  the  material,  in  the  analyses  carried  out  there,  this  energy 
is  used  to  compare  the  fracture  resistance  of  specimens  with  different  microstructures  under 
the  same  overall  remote  loading. 


Problem  Analyzed 

Computations  are  carried  out  for  a  center-cracked  specimen,  as  illustrated  in  Fig.  2.  The 
specimen  has  an  initial  height  of  2H  =  1.6  mm  and  an  initial  width  of  2 L  =  1.6  mm.  An 
initial  crack  of  length  2 =  0.4  mm  exists  along  the  C  axis.  Only  one  half  of  the  specimen 
(£'  >  9)  is  discretized  and  modeled  in  the  simulations  because  of  symmetry  with  respect  to 
the  C  axis.  Conditions  of  plain  strain  are  assumed  to  prevail.  The  small  region  in  front  of 
the  crack  tip  contains  microstructures  digitized  from  micrographs  of  actual  composite  ma¬ 
terials,  see  Fig.  3.  Inside  this  region,  the  material  inhomogeneities  and  material  phase  dis¬ 
tribution  are  explicitly  modeled.  The  dimensions  for  this  region  are  limited  by  the  memory 
sizes  of  the  Cray  T90  and  J90  computers  used  in  this  work.  As  seen  below,  the  particular 
dimensions  chosen  here  (100  by  200  pm)  are  much  larger  than  most  size  scales  in  the 
microstructures  analyzed.  Digitized  microstructures  of  actual  Al203/TiB2  composites  shown 
in  Fig.  4  are  used  in  the  analyses.  In  these  micrographs,  the  Ai20,  matrix  appears  lighter 
and  the  harder  TiB2  reinforcements  appear  darker.  These  microstructures  have  TiB2  volume 
fractions  between  15  to  25%.  In  microstructures  (tf-c),  the  TiB,  particles  are  embedded  in 
the  alumina  matrix.  The  average  grain  size  is  approximately  10-20  pm  for  ( a ),  8-10  pm 
for  ( b ),  and  1-2  pm  for  (c).  The  composite  microstructure  in  Fig.  4(d)  consists  of  a  popu¬ 
lation  of  isolated  Al,03  matrix  islands  surrounded  by  a  network  of  TiB?  reinforcements.  The 
isolated  matrix  areas  are  approximately  20  by  100  pm  in  size.  A  dispersion  of  small  A120^ 
particles  several  microns  in  size  are  scattered  in  the  TiB2  phase.  Clearly,  two  length  scales 
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PIG.  2 — Problem  analyzed:  a  rente  re d  c  racked  specimen  under  tensile  loading 


are  operative  in  this  microstructurc.  influencing  its  deformation  and  failure,  and  consequently 
determining  its  strength  and  fracture  toughness.  The  different  microstructural  phase  mor¬ 
phologies  are  the  result  of  different  conditions  during  hot  pressing,  see  Logan  |J4]  Partly 
because  'of  the  differences  in  their  microstructures,  these  materials  have  shown  a  range  of 
fracture  toughness  and  failure  resistance.  The  current  analyses  focus  on  the  effects  of  these 
different  phase  morphologies  on  the  fracture  behavior  of  the  specimen  shown  in  Fig.  2.  The 
effects  of  different  interfacial  bonding  strengths  that  may  also  result  from  the  different  proc¬ 
essing  conditions  are  not  specifically  analyzed  here.  Some  results  concerning  the  influence 
of  interphase  bonding  strength  on  failure  modes  in  these  materials  are  given  in  Zhai  et  al. 
[55].  In  the  current  analyses,  each  of  the  phases  behaves  hvperelastically.  Depending  on 
material  properties,  boundary  conditions,  and  loading,  fracture  is  possible  inside  each  of  the 
phases  and  along  the  interfaces  between  the  phases.  The  bulk  properties  of  each  finite  element 
are  those  for  either  the  grains  or  those  for  the  matrix.  The  properties  ol  each  segment  of 
potential  fracture  surface  are  specified  according  to  location  as  those  belonging  to  the  matrix, 
the  reinforcements  or  the  matrix /reinforcement  interfaces.  The  material  and  model  param¬ 
eters  are  listed  in  Table  1  (bulk  properties)  and  Table  2  (cohesive  surface  constitutive  prop¬ 
erties).  For  comparison  and  analysis  purposes,  the  speeds  for  the  longitudinal  stress  waves 
(ct),  the  shear  stress  waves  (c,)  and  the  Rayleigh  surface  waves  (c„)  are  also  listed  in  Table 
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TABLE  b — Bulk  material  properties. 


Compound 

Density 

kg/m3 

KIC, _ 

MPaVm 

E . 

GPa 

V 

m/s 

m/s 

m/s 

ALO, 

3990 

4.0 

340 

0.23 

9  939 

5885 

5379 

TiB2 

4520 

7.2 

500 

0.12 

10  694 

7027 

6267 

Homogenized 

ALO-,/Ti2B  Composite 

4120 

3.6 

415 

0.15 

10  300 

6609 

5936 

1.  The  choice  of  the  cohesive  law  parameters  assumes  that  crmax  =  EJ 10  for  each  constituent 
with  E  being  the  Young’s  modulus  and  (J>0  =  (1  -  u2)K]cfE  with  K}C  being  the  Mode  I 
fracture  toughness  of  the  materials  in  question.  The  parameters  for  the  Al203/TiB2  interfaces 
are  chosen  such  that  the  maximum  stress  amax  is  lower  than  that  for  either  the  matrix  or  the 
reinforcements.  Specifically,  both  aniax  of  interfacial  work  of  separation  <f>0  are  1-3  of  those 
for  the  matrix,  consistent  with  what  is  typically  reported  for  ceramic  composites. 

Materials  outside  the  crack-tip  region  are  assumed  to  be  homogeneous  and  are  assigned 
a  set  of  effective  parameters  representative  of  those  for  the  Al203/TiB2  ceramic  composite. 
Both  regions  are  discretized  in  the  same  manner,  using  both  the  bulk  and  the  cohesive  surface 
constitutive  descriptions.  The  specimen  is  stress-free  and  at  rest  initially.  Tensile  loading  is 
applied  by  imposing 'symmetric  velocity  boundary  conditions  along  the  upper  and  lower 
edges  of  the  specimen.  For  the  results  discussed  here,  the  imposed  boundary  velocity  is  V0  = 
10  m/s  for  each  edge  with  a  linear  ramp  from  zero  to  this  maximum  velocity  in  the  first 
0.0 1  of  loading.  All  other  specimen  surfaces  have  traction-free  boundary  conditions. 

Specifically,  the  loading  conditions  are 


(±— —  Va,  t  <  0.01  \rs 

0001  ,  -L<£<L  (12) 

±  V0  t  >  0.001  (jur 

Tl(V,  ±H,  0  =  0,  -L<^<L  (13) 

T\±L.  £%  t)  =  T2(±L,  t)  =  0,  -H  <  ?  <  H  (14) 

This  set  of  conditions  represents  the  loading  of  the  pre-crack  by  a  tensile  wave  with  a  stress 
amplitude  of  424  MPa  [(pcL)composite  V0]  and  a  linear  ramp  from  zero  to  that  value  in  0.01 

px 


TABLE  2 — Constitutive  parameters  for  cohesive  surfaces. 


Cohesive  Surface  Pair 

GPa 

T 

max  ’ 

GPa 

8„,  S„ 

nm 

<1>0> 

i/m2 

ai2o3 

34.0 

78.2 

0.5 

46.2 

TiB2 

50.0 

11.5 

1.0 

135.9 

Homogenized  Al203/Ti2B  Composite 

41.5 

95.5 

0.3 

33.8 

AI203/TiB2  Interface 

11.3 

26.1 

0.5 

15.4 
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Finite  Element  Method 

Finite  element  discretization  is  based  on  linear-displacement  triangular  elements  arranged 
in  a  “crossed-triangle”  quadrilateral  pattern.  Neighboring  elements  are  connected  through 
cohesive  surfaces.  Hence,  for  the  uniform  mesh  region  in  front  of  the  crack  tip  in  Fig.  3, 
the  cohesive  surfaces  are  initially  oriented  along  four  directions,  horizontal  (0°),  vertical  (90°), 
positive  and  negative  45°  (±45°).  Since  a  very  fine  mesh  is  used  (the  element  size  is  1  pm), 
arbitrary  fracture  paths  or  patterns  can  be  resolved.  When  the  finite  element  discretization 
of  the  displacement  field  is  substituted  into  the  principle  of  virtual  work  (Eq  1).  the  discre¬ 
tized  equations  of  motion  take  the  form 


M 


dW 

fit2 


=  R 


(15) 


where  V  is  the  vector  of  nodal  displacements,  M  is  the  nodal  mass  matrix,  and  R  is  the 
nodal  force  vector  consisting  of  contributions  from  the  bulk  elements  and  the  cohesive  sur¬ 
faces.  A  lumped  mass  matrix  is  used  in  Eq  15  instead  of  the  consistent  mass  matrix  for 
reasons  of  efficiency  and  accuracy  during  explicit  time-integration  Krieg  and  Key  [36).  The 
explicit  time-integration  scheme  based  on  the  Newmark  P-method  with  p  =  0  and  y  =  0.5 
Belytschko  et  al.  [37]  is  employed  to  integrate  Eq  15. 


Results 

To  illustrate  the  progression  of  fracture  in  the  four  microstructures,  the  distributions  of  g22 
at  several  different  times  for  four  calculations  using  these  microstructures  are  shown  on  the 
corresponding  deformed  configurations  in  Figs.  5-8,  respectively.  The  microstructura!  phase 
boundaries  are  outlined  by  solid  dark  lines  for  visualization  of  the  phase  morphologies.  In 
order  to  quantify  the  progress  of  crack  propagation,  the  new  crack  surface  areas  (or  crack 
lengths  in  two  dimensions)  generated  in  the  two  phases  and  along  matrix /reinforcement 
interfaces  are  calculated  separately  and  analyzed.  The  crack  length  histories  for  these  cases  , 
are  shown  in  Figs.  9 a-d,  respectively.  The  crack  lengths  in  the  matrix  and  along  the  interfaces 
as  well  as  the  total  crack  length  f(0  in  the  composite  are  shown.  Note  that  f(0>  =  0,  therefore, 
crack  lengths  referred  to  here  concern  new  crack  surfaces  generated  and  are  independent  of 
the  initial  pre-crack.  Since  crack  propagation  occurs  mainly  in  the  matrix  and  along  the 
phase  interfaces,  the  crack  length  in  the  TiB,  reinforcements  is  relatively  small  and  therefore 
not  shown.  The  total  crack  length  histories  f(/)  shown  in  Fig.  9  are  differentiated  with  respect 
to  time  to  yield  the  rate  at  which  new  crack  surfaces  are  generated.  Note  that  f  is  the  arc 
length  along  crack  paths.  The  time  derivative  di'fdt  may  not  always  represent  crack  speed 
due  to  possible  occurrences  of  multiple  crack  fronts  or  simultaneous  microcrack  formation. 
However,  it  can  be  used  as  a  relative  measure  of  the  extent  of  damage  in  the  materials.  On 
the  other  hand,  when  fracture  occurs  in  the  form  of  a  single  crack,  di'fdt  represents  the  lo¬ 
cal  crack  speed  as  a  function  of  time.  For  comparison  purposes,  the  apparent  crack  speed, 
dafdt,  where  a  is  the  projection  of  l’  in  the  (horizontal)  direction,  is  also  calculated.  The 
apparent  crack  length  a  is  taken  to  be  nondecreasing,  therefore,  dafdt  >  0.  The  apparent 
crack  speed  is  what  is  most  often  measured  in  experiments.  The  profiles  of  both  di'fdt  and 
dafdt  for  the  four  cases  are  shown  in  Fig.  10.  The  oscillations  in  these  profiles  show  the 
dynamic  and  unsteady  nature  of  crack  growth  in  the  inhomogeneous  material  microstructure, 
although  some  of  the  noise  may  be  attributed  to  errors  in  the  numerical  differentiation  of 
(’(/)  and  a(/). 
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Crack  Propagation  in  a  Heterogeneous  Al203/TiB2  Microstructure 

First,  the  results  for  the  microstructure  in  Fig.  4 a  are  discussed.  Figure  5  and  Fig.  9 a  show 
that  crack  initiation  occurs  at  approximately  0.14  jxj.  The  crack  first  zigzags  in  the  matrix 
between  small  TiB2  particles  and  then  arrests  when  the  crack  tip  impinges  on  a  large  grain 
at  point  A,  Figs.  5a-b .  This  grain  represents  a  stronger  obstacle  to  crack  growth  due  to  its 
higher  bulk  and  cohesive  strengths.  As  a  result,  the  crack  is  essentially  stationary  between 
0.15  and  0.20  \xs,  see  Figs.  9 a  and  10a.  There  is  an  increase  in  stress  intensity  in  front  of 
the  crack  tip  during  this  period  as  indicated  by  the  size  of  the  area  with  high  stresses.  The 
attainment  of  higher  stresses  allows  the  crack  to  circumvent  the  impeding  grain  through  the 
separation  of  grain /matrix  interface,  Fig.  5c-d.  Part  of  the  debonding  process  involves  pri¬ 
marily  local  tangential  (shear)  displacement  along  the  interface.  After  the  grain  at  A,  sub¬ 
sequent  propagation  is  primarily  through  the  matrix  and  interfaces  between  small  particles 
and  the  matrix.  The  crack  eventually  branches  into  two  propagating  tips  at  approximately 
0.204  |xj  and  their  overall  characteristics  are  the  same  as  those  for  the  parent  crack.  Fig.  5 d. 
It  can  be  seen  in  Fig.  9 a  that  more  crack  surfaces  are  created  along  the  phase  boundaries 
than  in  the  matrix,  reflecting  the  effect  of  the  weaker  bonding  between  the  phases  compares 
with  the  cohesive  strengths  of  the  constituents. 

In  Fig.  10a,  the  average  apparent  crack  speed  is  approximately  3300  m/s  or  61%  of  the 
Rayleigh  wave  speed  for  the  matrix  which  is  the  more  compliant  material  in  the  two  con¬ 
stituents  in  the  composite.  The  local  crack  speed,  on  the  other  hand,  is  higher.  Before  crack 
branching  at  approximately  0.204  | ±s,  the  propagation  of  a  single  crack  is  observed  (Fig. 
5a-b)  and  dMt  represents  the  local  crack  tip  speed.  The  average  value  of  this  local  speed  is 
approximately  6000  m/s.  The  plot  shows  that  this  speed  approaches  transonic  levels,  that 
is,  above  the  Rayleigh  wave  and  shear  wave  speeds  of  the  matrix  phase.  Note  that  after 
branching  two  crack  fronts  are  involved  and  the  local  crack  tip  speed  for  each  is  one  half 
of  what  is  seen  in  the  profile.  Therefore,  the  mean  local  crack  tip  speed  is  approximately 
5500  m/s  for  each  branch,  well  below  the  longitudinal  wave  speeds  for  the  constituents 
(9939  and  10  694  m/s,  respectively)  and  slightly  lower  than  the  mean  local  tip  speed  for 
the  crack  before  branching.  These  results  are  consistent  w'ith  the  findings  reported  in  Xu  and 
Needleman  [3839]  and  the  experimental  measurements  reported  in  Lambros  et  al.  [40], 
Huang  et  al.  [4J]  and  Rosakis  et  al.  [42]  concerning  crack  propagation  along  bimaterial 
interfaces.  Experimentally,  apparent  crack  speeds  in  the  intersonic  range  have  been  measured 
in  carbon  fiber/polymer  matrix  composites  under  impact  loading,  Rosakis  et  al.  [43]. 


Effects  of  Different  Microstructural  Morphologies 

The  results  in  Figs.  5-8  show  that  crack  path  is  significantly  influenced  by  the  micro- 
structural  phase  morphologies.  Fracture  occurs  primarily  along  phase  boundaries  and  inside 
the  matrix.  Continuous  and  favorably  oriented  interfaces  (parallel  to  the  direction  of  crack 
propagation)  facilitate  crack  growth,  see  Fig.  M.  Crack  branching  occurs  at  different  times 
and  different  distances  of  crack  extension  from  the  initial  crack  tip.  Also  note  that  no  branch¬ 
ing  is  seen  in  Fig.  6  for  microstructure  ( b ).  Two  other  calculations  are  carried  out  under  the 
same  loading  conditions  using  uniform  matrix  or  grain  material  properties  for  the  micro- 
structural  region  in  front  of  the  crack  tip.  The  results  show  that  in  both  cases  crack  branching 
occurs  after  very  short  distances  of  crack  propagation  in  the  direction  from  the  pre-crack 
tip.  Both  distances  are  smaller  than  those  observed  in  the  composites  here.  These  results 
demonstrate  that  material  heterogeneities  retard  or  prevent  crack  branching  in  composite 
microstructures.  This  observation  is  consistent  with  the  same  finding  in  Xu  et  al.  [31]  for 
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materials  with  inhomogeneous  cohesive  strengths  but  otherwise  homogeneous  bulk 
properties. 

The  time  histories  of  the  local  crack  speed  deleft  and  apparent  crack  speed  dafdt  for  all 
four  cases  shown  in  Fig.  10  indicate  a  range  of  crack  initiation  times  for  the  microstructures, 
from  0.12  to  0.19  p5.  The  microstructures  have  a  clear  impact  on  the  fracture  initiation  time. 
Temporary  arrest  of  crack  propagation  due  to  impediments  of  larger  particles  is  also  ob¬ 
served.  see  Fig.  9a, d  and  Fig.  10 a,d.  The  apparent  crack  speed  seems  to  oscillate  around  a 
mean  value  of  approximately  3300  m/s  or  61%  of  the  Rayleigh  wave  speed  for  the  matrix 
material  in  all  microstructures.  This  average  seems  to  be  independent  of  the  differences  in 
the  microstmctural  morphologies  and  phase  sizes  considered.  In  contrast,  the  local  crack 
speed  seems  to  be  more  significantly  affected  by  the  microstructural  phase  distributions.  In 
Fig.  6  for  microstructure  ( b ),  no  branching  occurs  and  the  local  speed  reaches  well  into  the 
intersonic  range,  that  is  the  mean  crack  speed  is  well  above  the  shear  wave  speeds  of  the 
constituents  and  approaches  9000  m/s  or  approximately  90%  of  the  longitudinal  wave  speed 
of  the  matrix  material,  see  Fig.  10/;.  Recall  that  this  local  speed  represents  the  rate  of  change 
of  arc  length  along  crack  paths  and  is  not  subject  to  the  limit  of  physical  wave  speeds  such 
as  cL  and  cs  of  the  bulk  constituents.  In  Figs.  5,  7,  and  8,  crack  branching  occurs  and  the 
mean  local  crack  speed  for  each  tip  Figs.  10^,  c,  and  d  is  approximately  6000,  5800  and 
7000  m/s  for  microstructures  ( a ),  (c),  and  ( d ),  respectively.  These  values  are  lower  than  that 
in  Fig.  10/;.  It  appears  that  branching  reduces  local  crack  tip  speeds  while  allowing  overall 
rate  of  surface  generation  to  increase. 


Energy  Release  Rate 

The  evolutions  of  the  kinetic  energy  :k,  strain  energy  W,  and  cohesive  surface  energy  cf> 
in  the  whole  specimen  for  the  four  cases  are  shown  in  Fig.  11.  The  times  of  crack  initiation 
are  clearly  signified  by  upturns  in  the  curves  for  in  ( a ),  (c),  and  (d).  In  Fig.  1  lc,  the  crack 
arrests  shortly  after  initiation.  Earlier  crack  initiation  allows  more  crack  surfaces  to  be  gen¬ 
erated,  resulting  in  more  energy  to  be  expended  in  the  process  (higher  <$>)  and  less  strain 
energy  to  be  stored  in  the  specimen  (lower  V).  This  is  the  case  in  Figs.  1 1 a-c  in  which  the 
crack  initiation  times  are  between  0.124  and  0.144  p s .  On  the  other  hand,  delayed  crack 
initiation  (Figs.  1  Id,  approximately  0.192  ps)  corresponds  to  less  new  crack  surfaces  and 
less  energy  consumption  (lower  <!>),  leading  to  more  strain  energy  stored  in  the  specimen 
(higher  tf).  Although  crack  branching  has  a  clear  signature  in  the  profiles  for  and  X  in 
Fig.  I  la  for  microstructure  (a),  such  an  effect  can  not  be  easily  identified  in  Figs.  II b-d 
for  other  microstructures. 

To  analyze  the  energy  expended  during  fracture,  the  total  crack  surface  energy  4>crart  is 
calculated  for  each  case.  This  energy  is 


‘I’crack  =  /  4>0 (MS  (16) 

Av'raok 

where  Scr3ck  is  defined  as  the  part  of  the  internal  surface  Sinl  that  satisfies  <J>(A)  ^  0.99<|)o.  In 
the  current  analysis,  <I>crack  is  calculated  for  the  first  350  pm  of  the  total  crack  length  t’(/),  or 
for  0  ^  f(/)  —  350  pm.  Note  that  it  takes  different  amounts  of  time  to  reach  this  crack  length 
in  different  microstructures.  The  calculated  values  of  0>cnck,  the  apparent  crack  length  a,  and 
the  time  t  at  a  total  crack  length  of  (’  =  350  pm  for  the  four  microstructures  are  listed  in 
Table  3.  Figure  12  compares  the  histories  of  <I>crack  as  function  of  f  and  the  apparent  crack 
length  a  for  the  four  cases  considered.  This  figure  and  the  table  show  that  there  is  a  significant 


of  energy  evolutions  for  the  four  micro  structures. 
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TABLE  3 — Average  energy  release  rates  for  different  microstructures. 


Microstmcture 

Average  size 
of  TiB2,  p.m 

^cfacfc’ 

x  10-3  y 

Time, 

a(t)  at 

C  =  350  |xm 

d'tWk 
df  1 
J/m 

^cac., 

da 

J/m 

(a) 

10-20 

12.5 

0.23 

122 

35.7 

102.5 

(b) 

8-10 

1 1.0 

0.19 

142 

31.4 

77.5 

( c ) 

1-2 

10.3 

0.22 

102 

28.6 

98.0 

id) 

20  X  100 

12.1 

0.24 

110 

34.3 

11 0.0 

variation  of  fracture  energy  among  the  four  microstructures.  The  waviness  in  the  curves 
reflects  the  different  values  of  work  of  separation  <j>()  in  different  constituents  and  along 
matrix /reinforcement  interfaces.  The  energy  release  rate  or  energy  consumption  per  unit 
length  of  crack  generated  is 


„  _  dfrcncVc 
dt 


(17) 


CRACK  LENGTH  (^m) 

FiG.  12 —Fracture  energy  4>cwi  as  functions  of  local  and  apparent  lengths  for  the  four 
microstructures. 


ZHAI  AND  ZHOU  ON  MODELING  CRACK  GROWTH  IN  CERAMIC  195 


When  fracture  occurs  in  the  form  of  a  single  crack  propagating  solely  in  the  horizontal 
direction,  f  =  a,  «  =  d<bcrack/aC  =  d&enc J da  represents  the  energy  release  rate  for  crack 
propagation  in  the  regular  fracture  mechanics  sense.  Figure  13  compares  histories  of  the 
instantaneous  energy  release  rate  defined  in  Eq  17  as  a  function  of  C  for  the  four  cases 
analyzed.  In  these  profiles,  there  are  three  characteristic  levels  which  correspond  to  the  works 
of  separation  inside  the  grains  (4>o  =  135.9  J/m2),  inside  the  matrix  (<|>0  =  46.2  J/m2)  and 
along  the  matrix /reinforcement  interfaces  (4>0  —  15.4  J/m2).  Instantaneous  8  values  between 
the  above  three  levels  occur  when  simultaneous  fracture  happens  at  multiple  sites  with  dif¬ 
ferent  values  of  4>0.  The  profiles  in  the  plot  show  that  such  incidents  occur  quite  frequently 
for  all  four  microstructures.  Most  of  these  are  microcrack  sites  along  the  main  crack  path  or 
paths  seen  in  Figs.  5—8,  since  all  four  cases  show  only  one  main  crack  (Fig.  6)  or  one  crack 
followed  by  two  branches  and  no  large  scale  microcracking  is  observed.  It  can  be  seen  that 
q  oscillates  mainly  between  46.2  J/m2  and  15.4  J/m2,  providing  a  quantitative  measure  for 
the  observation  made  from  Figs.  5—9  that  crack  propagation  is  primarily  through  the  matrix 
and  the  phase  boundaries.  To  facilitate  comparison,  the  average  values  of  3<Fcrack/dC  and 
d(Fcrack/dfl  corresponding  to  the  four  profiles  in  Fig.  13  are  listed  in  Table  3*  Note  that  it 
takes" different  amounts  of  time  for  fto  reach  350  p.m  in  different  microstructures.  For  higher 
resistance  to  fracture,  it  is  desirable  to  induce  cracks  to  go  through  the  hard  TiB2  reinforce¬ 
ments,  as  seen  in  Fig.  5 a  and  indicated  by  the  spikes  of  q  values  up  to  135.9  J/nrr  in  Fig. 
13a.  It  can  be  seen  That  as  the  average  size  of  TiB2  reinforcements  is  decreased  from  mi¬ 
crostructures  (a)  to  (c)  the  average  energy  release  rate  d$cnckf  dt  decreases,  reflecting  de¬ 
creasing  amounts  of  crack  growth  in  higher  strength  TiB2  or  matrix  in  Figs.  13 b  and  13c. 
Indeed,  Fig.  13c  shows  that  no  fracture  occurs  inside  TiB2  in  microstructure  (c)  for  the 
duration  analyzed.  Also,  there  seems  to  be  a  higher  fraction  of  fracture  along  matrix/grain 
interfaces  that  has  the  lowest  work  of  separation  among  the  three  types  of  fracture  sites.  The 
larger  sizes  of  the  TiB2  reinforcements  in  microstructure  (d)  cause  its  average  value  of 
d<j>crac Jdilo  be  higher.  Clearly,  there  is  significant  fracture  inside  the  TiB2  phase,  as  indicated 
in  Fig.  13d.  However,  the  existence  of  continuous  and  favorably  oriented  TiB2/Al203  bound¬ 
aries  inherent  in  the  morphology  of  a  network  of  TiB2  filled  with  A1203  colonies  also  allow 
cracks  to  propagate  along  the  interfaces  with  relatively  low  hindrance. 

Discussions 

The  results  discussed  above  show  that  it  is  possible  to  change  the  fracture  resistance  of 
brittle  composites  by  altering  microstructural  morphologies  through  favorable  material  proc¬ 
essing  conditions.  It  is  conceivable  that  energy  release  rates  higher  than  those  calculated 
here  may  be  achieved  by  other  microstructural  morphologies  not  considered  since  the  four 
variations  analyzed  represent  only  a  small  sample  of  microstructural  morphologies  that  can 
be  produced  by  various  methods  of  material  synthesis.  This  analysis  did  not  consider  the 
effects  of  varying  the  bulk  properties  for  each  phase  on  the  failure  behavior  of  the  compos¬ 
ites.  There  are  ample  opportunities  for  performance  enhancement  through  the  use  of  con¬ 
stituents  with  dissipation  mechanisms.  The  effects  of  various  elasto-plastic  constituents  on 
fracture  should  be  characterized  in  future  analyses.  Interfacial  bonding  strength  between  the 
phases  is  another  factor  that  can  significantly  influence  the  behavior  of  the  materials.  In  a 
related  analysis,  it  is  found  that  under  certain  loading  rates  this  interphase  bonding  plays  a 
role  in  determining  the  characteristics  of  fracture  or  the  mode  of  fracture  in  the  TiB2/Al203 
composites,  Zhai  et  al.  [35].  Specifically,  when  weak  interfacial  strength  exists,  microcrack 
initiation  and  growth  are  the  principal  mode  of  failure;  whereas  when  strong  interfacial 
strength  is  derived  from  material  processing,  the  advancement  of  a  dominant  crack  and  crack 
branching  are  observed.  It  is  expected  that  the  fracture  behavior  of  heterogeneous  materials 
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is  also  loading  rate  sensitive  due  to  the  rate-dependence  of  bulk  constituent  behavior  and 
inertia.  These  issues  shall  be  addressed  in  future  publications. 

A  micromechanical  model  that  provides  explicit  account  for  arbitrary  microstructures  and 
arbitrary  is  used  here.  The  approach  incorporates  both  a  constitutive  law  for  the  bulk  solid 
constituents  and  a  constitutive  law  for  fracture  surfaces.  The  interfacial  separation  law  can 
be  regarded  as  a  phenomenological  characterization  of  atomic  interaction  forces.  This  ap¬ 
proach  is  especially  appropriate  for  analyzing  microscopic  fracture  over  a  range  of  length 
scales  because  (i)  material  separation  is  a  natural  outcome  of  bulk  constitutive  response, 
interfacial  behavior,  microstructure,  and  loading;  and  (ii)  the  formulation  is  free  from  any 
fracture  initiation,  growth,  and  coalescence  criteria  that  are  necessary  in  most  other  contin¬ 
uum  analyses,  therefore  circumventing  the  limitations  of  failure  criteria  that  are  valid  only 
over  certain  length  scales,  such  as  continuum  criteria  based  on  the  existence  of  /Gfields.  This 
framework  allows  for  a  unified  treatment  of  mix-mode  fracture  without  resorting  to  the 
differentiation  of  fracture  into  separate  modes.  The  elimination  of  the  need  for  ad  hoc  criteria 
is  accompanied  by  the  introduction  of  the  need  for  interfacial  separation  laws  at  various 
microscopic  settings  of  interest.  Therefore,  an  important  task  is  the  development  of  interfacial 
constitutive  laws  appropriate  for  specific  materials  and  length  scales.  For  specific  forms  of 
cohesive  laws,  the  determination  of  cohesive  parameters  for  constituents  is  another  relevant 
issue.  For  example,  over  a  certain  size  scale  the  calibration  of  model  parameters  by  assuming 
4>0  =  (1  -  i r)K2lc!E  may  be  justifiable  for  brittle  constituents.  However,  the  same  assumption 
may  become  indefensible  if  the  objective  is  to  quantify  nano-scale  fracture  mechanisms. 
Indeed,  there  is  a  serious  lack  of  experimental  characterization  as  well  as  understanding  for 
interfacial  constitutive  behaviors  on  the  sub-micron  length  scales.  Issues  include,  among 
others,  dissipation  in  the  separation  process  and  the  gradient-dependence  of  traction  along 
interfaces.  Kim  [44]  proposed  a  relation  for  interface  decohesion  that  accounts  for  nonlocal 
effects  on  the  nano-scale.  Significant  challenges  and  opportunities  exist  for  experimental  and 
analytical  research  on  nano-scale  and  micro-scale  interfacial  separation  mechanisms.  As  such 
separation  laws  are  made  available  and  as  long  as  appropriate  bulk  response  laws  are  used 
including  the  quasi-continuum  formulation  of  Tadmor  et  al.  [45],  the  unique  framework  off 
analysis  employed  here  represents  a  systematic  and  consistent  approach  for  explicit  modeling 
of  fracture  over  a  wide  range  off  length  scales. 


Conclusions 

The  dynamic  failure  behavior  in  a  class  of  Al203/TiB2  composites  is  analyzed  in  the 
context  of  fracture  in  a  centered-cracked  specimen  under  remote  tensile  loading.  Numerical 
simulations  are  carried  out  under  the  conditions  of  plane  strain.  A  micromechanical  frame¬ 
work  of  analysis  has  been  developed  and  used,  providing  an  explicit  account  for  arbitrary 
material  microstructural  morphologies  and  resolution  for  arbitrary,  unconstrained  fracture 
patterns.  The  approach  combines  constitutive  description  of  bulk  material  responses  and 
characterization  of  fracture  surface  cohesion.  The  framework  of  analysis  allows  mixed-mode 
fracture  behavior  of  materials  with  composite  microstructures  to  be  examined  without  re¬ 
sorting  to  any  mode-dependent  fracture  initiation  or  propagation  criteria.  In  this  analysis, 
both  the  bulk  and  the  interfacial  constitutive  laws  are  assumed  to  be  hyperelastic.  The  results 
of  numerical  simulations  demonstrated  the  microscopic  phase  distribution  and  phase  size 
scale  can  significantly  influence  the  fracture  behavior.  Due  to  their  higher  cohesive  strength 
and  therefore  higher  energy  release  rate,  TiB2  reinforcements  represent  stronger  obstacles  to 
the  growth  of  cracks  and  increase  the  energy  required  in  generating  the  same  amount  of 
crack  surfaces.  The  increase  in  fracture  resistance  is  higher  in  microstructures  with  larger 
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TiB,  particles.  A  25%  variation  in  average  energy  release  rate  is  observed  for  the  four 
microstructures  of  actual  A1203/TiB2  composites  prepared  by  high-temperature  synthesis. 
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